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An  equation  for  rcr  of  sand  beds  in  the  presence  of  interstitial  clay  particles  was 
developed  based  on  the  balance  of  forces  on  a  sand  particle  at  the  surface  of  the  bed.  The 
equation  was  applied  to  data  collected  from  a  series  of  laboratory  flume  experiments  on 
medium  sand  beds  in  order  to  examine  the  variation  of  xcr  with  clay  weight  fraction  \|/ .  A 
quasi-hydrodynamic  model  of  interstitial  clay  flow,  possessing  some  of  the  features  of 
hydrodynamic  lubrication,  was  developed  to  explain  the  experimental  results  through 
changes  in  the  lubricating  effect  of  the  clay  suspension  between  the  sand  grains  as  the 
clay  fraction  increased. 

Based  on  the  application  of  the  equation  for  rcr  to  experimental  data,  insight  was 
gained  on  the  behavior  of  ij/min ,  the  weight  fraction  at  which  the  minimum  xcr  occurs  at 
the  point  of  maximum  lubrication.  It  was  shown  that  \\imin  was  affected  by  the  frictional 
resistance  offered  by  the  clay  suspension  and  by  pore  water  salinity,  but  not  sand  size. 


For  a  bed  of  0.83  mm  sand  to  which  kaolinite  in  fresh  water  was  added  in  varying 
amounts,  measurements  yielded  i|/min  =  0.04.  This  result  was  explained  through  y^ 
being  the  fine  weight  fraction  at  which  the  thickness  of  the  clay  layer  is  close  to  the  size 
of  sand  grain  asperities  (~  20  |xm)  at  which  a  transition  from  thin-film  to  boundary  layer 
lubrication  occurs.  It  was  thus  shown  that  y^,,  is  not  the  space-filling  clay  weight 
fraction  as  postulated  by  previous  investigators. 

The  space-filling  clay  weight  fraction  was  found  to  be  \|/r,  the  weight  fraction  at 
which  the  shear  stress  once  again  equals  tco,  and  the  point  at  which  the  sediment  starts  to 
be  dominated  by  the  cohesive  component.  For  a  bed  of  0.83  mm  sand  to  which  kaolinite 
in  fresh  water  was  added  in  varying  amounts,  measurements  yielded  \|/r=  0.13  which  was 
found  to  be  the  point  at  which  the  pore  volume  not  occupied  by  clay  becomes  negligible 
(<  1%  of  the  total  volume). 

Through  the  application  of  a  shear  resistance  model  to  the  experimental  data 
collected  in  a  Schulze  ring  shear  tester,  it  was  concluded  that  the  strongest  lubricating 
effect  occurred  when  the  over-all  confining  stress  was  carried  by  both  the  sand  matrix 
and  the  thin  interstitial  clay  suspension. 

The  effect  of  salinity  on  zcr  and  \\/min  was  explained  through  a  change  in 
flocculation  and  a  consequent  change  in  the  number  of  edge-to-face,  particle-particle 
contacts. 

The  two  models  for  clay  lubrication  presented  in  this  study  gave  differing  values 
for  ij/min,  and  this  discrepancy  warrants  further  study.  The  quasi-hydrodynamic  model  for 
clay  lubrication  predicted  a  value  for  y^,,  for  kaolinite  of  0.04,  and  the  shear  resistance 
model  for  clay  lubrication  predicted  a  value  for  y^  for  kaolinite  of  0.01.    A  possible 
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cause  of  this  discrepancy  is  that  while  the  dimensional  analysis  that  leads  to  the  Stribeck 
curve  is  reasonable  for  the  "thick-film"  lubrication  regime  where  the  grains  are  not  in 
contact,  it  is  quite  possible  that,  in  the  lower  region  of  the  curve  where  the  grain  asperities 
become  significant,  other  variables  such  as  grain  asperity  size  and  the  load  carried  by  an 
individual  grain  might  become  contributing  factors.  Consequently,  the  application  of 
normal  loads  over  grains  with  differing  asperities  needs  to  be  conducted  to  determine  the 
effect  of  normal  load  on  v|/min . 
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CHAPTER  1 
INTRODUCTION 


1.1  Problem  Statement 

For  a  proper  management  of  coastal  and  estuarine  waters  it  is  essential  for 
engineers  and  scientists  to  accurately  predict  the  movement  of  mixtures  of  cohesionless 
and  cohesive  material.  There  are  reliable  methods  for  predicting  the  erosion  of  sand,  but 
not  for  mixed  sediments.  A  practical  problem  one  faces  with  sandy  beds  is  that  whereas 
engineers  and  scientists  typically  use  the  Shields'  diagram  to  predict  the  critical  shear 
stress  for  erosion,  recent  studies  have  shown  that  the  presence  of  even  small  fractions  of 
clay  in  the  sandy  bed  can  measurably  alter  the  critical  shear  stress.  For  example, 
Mitchener  and  Torfs  (1996)  attempted  to  characterize  the  erosional  behavior  of  mixed 
sediments  in  a  way  that  could  be  used  for  predictive  models  and  in  so  doing  made  some 
important  observations: 

(1)  The  critical  shear  stress  for  erosion  increases  when  mud  is  added  to  sand,  and  also 
when  sand  is  added  to  mud  (mud  is  <  62.5  [im).  The  addition  of  up  to  50%  sand  to  a  mud 
bed  can  typically  increase  critical  erosion  shear  stress  by  a  factor  of  2.  Conversely,  the 
addition  of  30%  mud  to  a  sand  bed  can  increase  the  critical  shear  stress  by  as  much  as  a 
factor  of  10. 

(2)  The  addition  of  mud  to  sand  significantly  increases  the  critical  erosion  shear  stress 
with  a  maximum  value  occurring  at  a  sand  content  between  50%  and  70%  by  weight 
(corresponding  to  a  30%  to  50%  mud  content). 


(3)  The  addition  of  mud  to  sand  not  only  delays  the  onset  of  erosion  by  increasing  the 
threshold  conditions  required  for  sediment  motion,  it  also  diminishes  the  erosion  rates 
once  this  occurs.  Similarly,  the  erosion  rate  of  mud  decreases  with  increasing  sand 
content.  The  erosion  rate  of  a  mixed  bed  can  be  reduced  by  a  factor  of  10  if  more  than 
20%  sand  is  added  to  a  mud  bed. 

The  above  findings  were  primarily  descriptive  in  nature,  but  provided  a  very 
thorough  characterization  of  the  erosional  behavior  of  sediment  mixtures. 

Panagiotopolous  et  al.  (1997)  looked  at  the  erosion  threshold  (i.e.,  critical  shear 
stress)  of  mixed  sediment  deposits  under  the  action  of  unidirectional  currents  and 
attempted  to  establish  the  critical  erosion  conditions,  in  terms  of  near-bed  current  speed 
and  bed  shear  stress  necessary  to  erode  each  sediment  mixture,  and  to  examine  the 
influence  of  clay  content  on  the  established  erosion  threshold  conditions.  The  results 
suggested  that  the  critical  stress  for  erosion  of  a  mixed  bed  increases  with  increased  mud 
content.  The  authors  argued  that  two  regimes  could  be  identified:  one  regime  for  mud 
content  less  than  30%  and  another  for  contents  greater  than  30%.  Under  steady  flow 
conditions,  a  linear  increase  in  the  threshold  shear  stress  occurred  for  low  mud  contents 
(<  30%).  These  observations  agreed  with  the  observation  that  a  clay  mineral  content  of 
5-10%  (by  weight)  is  high  enough  to  dominate  sediment  erodibility  characteristics. 

The  authors  applied  the  analysis  of  Wiberg  and  Smith  (1997)  to  mixed  sediment 
beds  and  used  this  application  to  suggest  that  at  low  mud  contents  (<  30%)  the  sand 
particles  continue  to  be  in  contact  with  each  other  so  that  pivoting  is  the  main  mechanism 
of  initiation  of  motion.  However,  they  noted  that  mud  fills  the  voids  as  the  fine  fraction 
is  increased,  and  in  particular  the  spaces  on  either  side  of  the  contact  point  of  the  two 


sand  particles  and  that  this  results  in  slightly  altering  the  pivoting  characteristics  and 
consequently  the  value  of  the  internal  friction  angle.  At  higher  mud  content  (>  30%),  a 
matrix  of  clay  particles  is  formed  which  incorporates  the  sand  and  silt  grains  at  which 
point  the  individual  sand  particles  are  no  longer  in  contact  with  each  other  and  pivoting 
stops  to  be  the  main  mechanism  for  initiation  of  motion.  As  a  result,  erosion  is  then 
controlled  mainly  by  the  resistance  of  the  clay  fraction. 

Torfs  et  al.  (2000)  developed  a  formula  for  the  critical  stress  for  mixture  erosion, 
Tcm,  and  the  rate  of  erosion  of  the  fine  fraction  through  a  diagnostic  evaluation  using  data 
on  the  erosion  of  selected  clay/sand  and  mud/sand  mixtures,  and  in  so  doing  made  some 
additional  important  observations: 

(1)  The  critical  shear  stress  of  the  mixture,  Tcm,  was  found  to  vary  non-monotonically 
with  the  fine  sediment  weight  fraction,  \|/.  Starting  with  a  sandy  bed,  as  \|/  increased,  at 
first  Tcm  seemingly  decreased  from  its  value  for  sand  and  reached  a  minimum  at  a  weight 
fraction  tj/mjn,  then  increased. 

(2)  It  appeared  that  \|/min,  the  fine  grain  weight  fraction  at  which  the  minimum  value  of 
critical  shear  stress  occurred,  may  be  a  measure  of  the  space-filling  concentration  of  fine 
particles  within  the  pores  of  the  sandy  bed  matrix. 

(3)  The  logarithm  of  the  rate  of  erosion  was  shown  to  vary  linearly  with  the  logarithm  of 
the  excess  bed  shear  stress,  and  this  variation  depended  on  \|/. 

(4)  It  was  their  opinion  that  the  highly  non-linear  dependence  of  the  critical  shear  stress 
on  the  fine-grained  weight  fraction  was  indicative  of  the  interactive  role  of  the  mixture 
components  in  controlling  erosion,  and  suggested  a  potential  route  for  improvement  in 
modeling  erosion  by  accounting  for  these  interactions.  In  that  context,  they  considered 


additional  experimental  observations  on  the  erosion  of  mixtures  in  the  vicinity  of  space- 
filling weight  (or  volume)  fraction  to  be  essential. 

1.2  Problem  Approach 

It  was  the  goal  of  this  study  to  develop  a  physics-based  explanation  for  the 

behavior  of  the  critical  shear  stress  of  sand  beds,  rcr,  in  the  presence  of  interstitial  clay 
particles.  As  suggested  by  the  previous  works  on  this  subject,  there  exists  a  definite  trend 
in  the  behavior  of  sand  beds  containing  pore  water  clay  particles  with  fine  fractions  less 
than  -0.15.  Starting  with  pure  sand,  as  the  weight  fraction  of  clay  in  the  pore  water 
increases  the  critical  shear  stress  decreases  and  reaches  a  minimum  value  at  a 
characteristic  weight  fraction,  then  increases.  This  work  has  attempted  to  both  further 
quantify  and  explain  this  behavior,  so  that  the  behavior  of  rcr  could  be  predicted.  This 
goal  was  accomplished  by  experimental  observations  on  the  erosion  of  sand  beds  in  the 
vicinity  of  the  aforementioned  range  of  clay  weight  fraction. 

To  accomplish  this  goal,  several  tasks  were  completed.  First,  an  expression, 
based  on  a  force-balance  for  a  sand  grain  at  the  surface  of  a  bed,  was  developed  for  the 
critical  shear  stress  of  sand  beds  in  the  presence  of  interstitial  clay  particles,  Tcr,  in  the 
range  of  fine  fraction,  \\i,  from  0.0  to  0.15. 

Second,  a  narrow,  unidirectional  flow  flume  with  a  false  bottom  designed  for 
emplacing  sediments  beds,  was  used  to  determine  the  relationship  of  the  critical  shear 
stress  of  the  sand  beds  as  a  function  of  clay  weight  fraction.  This  flume  was  calibrated  by 
comparing  clay  erosion  data  obtained  in  it  with  analogous  data  obtained  in  a  cylindrical 
apparatus  in  which  the  shear  stress  on  the  bed  could  be  measured  directly. 


Third,  the  expression  for  the  critical  shear  stress  of  sand  beds  in  the  presence  of 
interstitial  clay  particles  was  applied  to  the  results  from  the  flume  tests  to  determine  the 
effect  of  the  fine  grain  fraction  on  the  effective  particle  packing  angle,  Oeff. 

Fourth,  a  quasi-hydrodynamic  model  of  clay  lubrication,  a  mechanism  of 
interstitial  clay  flow  which  possesses  some  features  of  hydrodynamic  lubrication,  was 
developed  to  explain  the  interaction  between  sand  grains  and  the  fine  fraction. 

Fifth,  a  Schulze  ring  shear  tester  was  modified  for  use  with  clays  and  then  used  to 
explain  the  behavior  of  xcr  of  sand  at  low  values  of  \|/. 

Finally,  a  shear  resistance  model,  which  relates  to  the  nature  of  shear  resistance 
between  particles,  was  developed  to  explain  the  dependence  of  xcr  on  the  fine  grain 
fraction. 

1.3  Outline  of  Chapters 

This  study  is  presented  in  the  following  order. 

Chapter  2  is  devoted  to  a  discussion  of  an  analytical  and  experimental  background 
of  subjects  related  to  critical  shear  stress  of  non-cohesive  sediments,  cohesive  sediments, 
and  mixtures  of  the  two. 

Chapter  3  first  defines  incipient  motion  for  both  cohesionless  and  cohesive 
sediments,  then  develops  an  expression  for  incipient  motion  of  cohesionless  sediment, 
and  finally  presents  in  detail  the  derivation  of  an  expression  describing  incipient  motion 
of  mixtures  of  both  cohesive  and  non-cohesive  sediment. 

In  Chapter  4,  the  experimental  equipment,  material,  and  procedures  are  discussed 
in  detail. 


In  Chapter  5,  the  results  for  each  of  the  three  types  of  equipment  used  in  this 
study  are  presented.  This  chapter  is  organized  into  a  section  for  each  device,  with  each 
section  giving  first  a  table  of  the  results,  plots  if  pertinent,  and  then  an  analysis  of  the 
results  from  that  device. 

In  Chapter  6,  the  results  from  each  of  the  three  types  of  equipment  used  in  this 
study  are  analyzed.  This  chapter  is  organized  into  three  sections.  The  first  two  describe 
two  different  experimental  models  used,  and  the  final  section  gives  a  comparison  of  these 
models. 

In  Chapter  7,  conclusions  from  this  study  and  recommendations  for  future 
research  are  presented. 


CHAPTER  1 
LITERATURE  REVIEW 

2.1  Introduction 


Selected  studies  on  the  critical  shear  stress  and  erosion  of  cohesionless,  cohesive 
and  mixed  cohesionless/cohesive  soil  beds  are  reviewed  here  in  chronological  order  to 
provide  the  background  for  the  mechanisms  that  underlie  the  physical  processes  germane 
to  the  present  study. 

2.2  Einstein  and  El  Samni  (1949) 

The  purpose  of  this  study  was  to  measure  the  drag  and  lift  forces  acting  on 
particles  in  the  surface  layer  of  a  bed.  Spherical  plastic  balls,  6.9  cm  in  diameter,  were 
used  as  the  roughness  elements  in  a  recirculating,  unidirectional  flume.  In  order  to 
prevent  the  seepage  of  water  between  these  balls  and  the  bottom  of  the  flume,  only  the 
upper  halves  of  the  spheres  were  used  by  gluing  the  hemispheres  in  a  hexagonal  pattern 
to  the  bottom  of  the  flume.  The  drag  and  lift  components  of  the  average  dynamic  forces 
on  this  rough  surface  were  then  measured. 

Drag  is  the  resultant  force  in  the  direction  of  velocity,  and  lift  is  the  resultant  force 
in  the  direction  of  the  normal  to  velocity.  These  resultants  of  shear  stress  and  pressure 
distributions  can  be  obtained  by  integrating  the  effect  of  these  two  quantities  on  the  body 
surface: 


FD  =  [Pcos8dA+[rw  sinBdA 

J  .  J  .  (2.1) 

FL  =  -  J  Psin  G  dA  +  J  xw  cos0  dA 

where  FD  is  the  drag  force,  Fl  is  the  lift  force,  Tw  is  the  shear  stress  along  the  bed,  P  = 
pressure  at  any  point  on  the  surface  of  the  body,  dA  =  differential  element  of  surface  area 
for  the  body,  and  6  is  the  angle  to  the  tangent  of  the  body  surface.  The  lift  and  drag 
coefficients  are  defined  as: 

CL  =  -=—£ CD  =  — -^ (2.2) 

In  these,  the  two  effects  are  considered  together  «  net  drag  is  produced  by  both 
pressure  and  shear  effects,  and  so  is  net  lift.  pw  =  mass  density  of  the  fluid,  U  =  average 
fluid  velocity  over  the  cross  sectional  area  of  the  particle,  and  Ax  =  average  cross 
sectional  area  of  the  particle.  To  determine  Q,  the  measured  pressure  difference  between 
flume  side  wall  above  the  hemispheres  and  the  flume  floor  between  the  hemispheres  was 
described  in  the  following  form: 

( P  \     i 

*P  =  CL       V       "035  (2-3) 

where  wo.35  was  measured  0.35  hemisphere  diameters  from  the  theoretical  bed,  assumed 
to  be  a  distance  of  0.2^  below  the  plane  tangent  to  the  top  of  the  spheres  whose  diameter 
was  taken  as  ks.  The  coefficient,  CL,  was  found  to  have  a  constant  value  of  0.178  when 
and  only  when  the  flow  velocity  was  measured  at  a  distance  of  0.35  diameters  from  the 
theoretical  bed. 


To  determine  FD,  hence  CD,  the  investigators  used  a  logarithmic  velocity  profile: 


—  =  C2    +    5.75   log 


v*7 


(2.4) 


where  u  is  the  average  velocity  over  the  hydraulic  radius  R  and  Ci  is  a  constant.  The 
friction  velocity  w.  was  obtained  from  the  above  profile  and  then  Tb  =  pwu* .  This  shear 
stress  was  then  integrated  over  the  entire  bed.  The  results  for  the  average  shear  stress 
matched  previous  results.  The  coefficient  C2  was  found  to  have  the  value,  8.50,  as 
previously  determined  by  Keulegan  (1938).  Tests  were  also  conducted  of  the 
applicability  of  Equation  2.4  to  turbulent  flow  over  a  natural  gravel  bed.  In  these 
experiments,  very  close  agreement  was  achieved  if  the  roughness  ks  was  chosen  to  equal 
6.1  cm. 

2.3  Benedict  and  Christensen  (1967) 

The  objective  of  this  work  was  to  develop  a  physical  understanding  of  the  lift 
force  in  the  hydrodynamically  rough  range  of  flow.  The  investigators  developed  an 
analytical  approach  for  the  lift  force  on  idealized  beds  composed  of  hemispheres.  A 
potential  flow  solution  was  obtained  for  the  bed  and  linked  to  the  actual  velocity  profile 
over  the  bed  to  predict  pressure  and  hence  lift.  The  integration  of  the  theoretical  solution 
for  the  closely  packed  hemispheres  yielded  a  value  of  d  =  0.259. 

2.4  Christensen  (1975) 

This  work  was  meant  to  assign  a  probability  of  erosion  to  the  critical  shear  stress 
of  a  cohesionless  sediment  bed.  Bed  scour  is  usually  considered  to  begin  when  the  local 
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time-mean  bed  shear  stress  exceeds  the  critical  shear  stress  of  the  bed  material.  This 
leads  to  the  definition  of  the  critical  bed  shear  stress  as  the  time-mean  value  of  the  shear 
stress  when  erosion  is  initiated.  However,  the  forces  that  cause  a  grain  to  move  are  of  a 
stochastic  nature  due  to  the  turbulent  flow.  Therefore,  it  is  possible  to  observe  local  bed 
failure  due  to  an  instantaneous  velocity  fluctuation  causing  these  forces  to  exceed  the 
failure  condition  when  the  time-mean  bed  shear  stress  is  less  than  the  critical  shear  stress. 
It  is  therefore  appropriate  to  assign  a  probability  of  erosion  to  the  critical  shear  stress. 
The  investigator  started  with  the  criterion  for  incipient  motion  of  a  grain, 

=     tanq>  (2.5) 


Wb-FL 
where 

FD=xhaxd]  (2.6) 

FL=laxd)  (2.7) 

and 

Wb=(ys-Y)a2d3e  (2.8) 

in  which  Fd  is  the  hydrodynamic  drag  force,  Fl  is  the  hydrodynamic  lift  force,  Wb  is  the 
buoyant  weight,  ys  is  the  specific  weight  of  the  sediment  grain,/ is  the  specific  weight  of 
the  fluid,  (p  is  the  angle  of  repose,  r0  is  the  instantaneous  value  of  the  fluctuating  bed 
shear  stress,  and  A  is  the  instantaneous  value  of  the  fluctuating  lift  per  unit  bed  area;  both 
T0  and  A  correspond  to  incipient  motion.  a,\  and  da  are  shape  factors;  the  first  is  the  bed 
area  associated  with  one  grain  divided  by  the  grain  size  squared  and  the  other  is  the  grain 
volume  divided  by  the  grain  size  cubed.  It  is  important  to  note  that  drag  force  in 
Equation  2.5  is  the  absolute  value  because  it  is  possible  for  the  instantaneous  flow  to  be 
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from  right  to  left  even  though  the  time-mean  flow  is  from  left  to  right.    Consequently, 
two  failure  conditions  need  to  be  considered. 

Equations  2.6,  2.7,  and  2.8  were  substituted  into  Equation  2.5  resulting  in 


a'     A 


*J  =  — •     ,    "   2    '  (2-9) 


k2 


I r  +  COt<)0 


Fm 


The  instantaneous  lift  to  shear  stress  ratio  was  then  replaced  by  the  ratio  of  the 
corresponding  time-mean  values,  since  A  and  Tt  both  are  proportional  to  the  same 
instantaneous  velocity  squared.  The  absolute  value  of  the  instantaneous  bed  shear  stress 
corresponding  to  incipient  motion  may  be  written 

W  =  CBPwy  (2-10) 

where  Co  is  assumed  constant.  The  instantaneous  velocity  w  is  measured  at  the  top  of  the 
grain,  which,  according  to  Einstein  and  El  Samni  (1949)  is  located  at  approximately  0.2 de 
height  above  the  theoretical  bed.  When  the  velocity  u  is  expressed  by  its  time-mean  and 
fluctuating  component,  u  and  u  ,  the  corresponding  absolute  value  of  the  bed  shear 
stress  becomes 

W  =  7lV(1  +  5.»)2  (2-11) 

with  the  time-mean  value 


*t  -  CD  Pw 


-if      /■  -  \2\ 


1  + 


V  u    ) 


(2.12) 


where  Su  is  the  dimensionless  value  of  the  standard  deviation  of  the  velocity  fluctuation 
au,  defined  as  aju  ,  and  n  is  the  normalized  velocity  fluctuation  u'lau. 
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An  analogous  treatment  was  done  for  the  time-mean  value  of  lift  and  substituted 
into  Equation  2.9  along  with  Equation  2.1 1  to  obtain 

T"    4  +  nSj  =%-%'-&-  (2.13) 


1+5.^ 


a,    X 


— +  cot<p 


From  Equation  2.13  it  is  seen  that  the  critical  shear  stress  of  a  horizontal  bed  may  be 


written 


F„  = 


cc2{Ys-Y)de      i+s2u 


a,    A 


—  +  cot(p 


(1  +  nSj 


(2.14) 


The  investigator  then  introduced  an  expression  for  the  ratio  of  the  time-mean  values  of 
lift  and  shear  based  on  Einstein's  (1943)  observations  to  obtain  an  expression  for  critical 
shear  stress 


*cr=0Cr(Ys  ~  Y)de 


(2.15) 


where 


e  =^ 


i  +  s7. 


a, 


cot  cp  +  0.556 


In 


1   + 


10.4 


/j 


(1  +  nSuf 


(2.16) 


The  symbol  r  in  Equation  2.16  is  the  roughness/grain-size  ratio  defined  as  the  equivalent 
sand  roughness  ks  divided  by  the  effective  grain-size  de. 

Equations  2.15  and  2.16  were  found  to  be  in  good  agreement  with  most 
observations,  with  Su  =  0.164  and  n  =  3.09  corresponding  to  the  probability  or  erosion  Pr 
=  0.001  related  to  the  Gaussian  distribution  of  «'. 
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2.5  Wiberg  and  Smith  (1988) 


Expressions  were  derived  for  the  critical  shear  stress  for  beds  of  both  uniform  and 
heterogeneous  sediment  under  unidirectional  flow.  The  investigators  start  by  developing 
a  relationship  for  the  critical  shear  stress  for  a  sediment  particle  at  the  surface  of  a 
uniform  bed: 

2  1  (tan<I>cos06-sin06) 


T„  =  ■ 


CD«3     lf2\ 


1  + 


V 


tanO 


(2.17) 


where  a?  is  a  parameter  representing  the  geometry  of  the  grain,  /  2(z/z0)  comes  from  the 
expressions  for  FD  and  Fl  and  represents  the  cross-sectional  averaged  velocity  profile 
function,  z  =  height  above  the  bed,  za  =  bottom  roughness  parameter,  <J>  is  the  particle 
packing  angle  (or  bed  pocket  geometry  angle),  a  measure  of  the  difficulty  of  moving  a 
sediment  grain  up  and  out  of  the  pocket  in  which  it  is  resting,  and  &b  is  bed  slope.  The 
investigators  use  the  following  expression  of  Reichardt  (1951)  for  the  velocity  profile, 
which  provides  a  transition  between  the  viscous  sub-layer  and  the  logarithmic  flow 
above,  and  is  in  agreement  with  measurements  (Schlichting,  1979): 


u(z)=ut 


— ln(l  +  Kz+)-c 

K 


l_^-".6_J_e-0.33Z+ 

11.6 


yj 


(2.18) 


where  z+  =  u*  zjv  =  R,z/ks  ,  z+0  =  u.zjv  =  R,zjks  and  /?,is  the  roughness  Reynolds 
number  Rt  =utks/v ,  where  ks  is  the  scale  length  for  the  bed  roughness.  It  should  be 
noted  that,  with  this  expression,  all  of  the  terms  in  Equation  2.17  can  be  either  measured 
or  calculated  with  the  exception  of  the  particle  packing  angle,  O. 
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The  investigators  then  develop  an  expression  for  the  critical  shear  stress  for 
heterogeneous  beds.  They  argue  that  there  are  two  lengths  scales  for  such  beds  while 
there  is  only  one  (the  grain  diameter)  for  homogeneous  beds.  The  diameter  of  the  particle 
size  fraction  of  interest  is  the  length  scale  for  lift,  drag,  and  gravity,  specifically  as  it 
governs  how  far  into  the  flow  the  grain  protrudes.  The  flow  velocity,  however,  also 
depends  on  the  roughness  scale  ks  of  the  surrounding  bed.  The  presence  of  these  two 
length  scales  is  primarily  responsible  for  the  difference  between  the  computed  critical 
shear  stress  for  sediment  on  uniform  versus  heterogeneous  beds.  The  investigators  bridge 
this  gap  using  an  expression  for  O.  They  argue  that  the  particle  angle  of  repose  changes 
for  a  bed  composed  of  mixed  sizes.  For  example,  it  is  easier  for  a  large  grain  to  roll  over 
a  bed  composed  of  smaller  particles  than  for  a  small  grain  to  roll  over  a  bed  of  larger 
particles.  This  is  indicated  by  an  increase  in  the  value  of  <t>  as  the  ratio  of  grain  diameter 
to  bed  roughness  length,  d/ks,  decreases.  The  investigators  use  an  expression  derived  from 
the  work  of  Miller  and  Byrne  (1966): 


O  =  cos" 


d/k+l 


(2.19) 


where  z*  is  the  average  level  of  the  bottom  of  the  almost  moving  grain. 

2.6  Dade  and  Nowell  (1991) 

It  was  the  purpose  of  this  work  to  develop  an  expression  for  critical  bed  shear 
stress  using  the  force  balance  of  Wiberg  and  Smith  (1988),  but  substituting  a  term  for  the 
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inter-particle  cohesive  force  to  obtain  an  equation  for  entrainment  from  flat  cohesive  mud 
beds: 


T„  = 


40pwv' 


tan<$>d' 


'j3  »„„2/^/ 


gpsd'  tan'O 


pwv2     240 


l+3(l-cos4>) 


A 


giPs-PJd 


+  1     -1 


(2.20) 


where  %  is  the  interfacial  mud  shear  stress,  ps  is  the  particle  density,  v  is  the  fluid 

kinematic  viscosity,  pw  is  the  fluid  density,  g  is  the  acceleration  due  to  gravity,  and  the 
average  particle  packing  angle,  <J>  is  65°.  Relative  contributions  of  either  particle  cohesion 
or  submerged  weight  to  mud  erosion  resistance  are  embodied  in  the  dimensionless  term 
?bl{gpsd). 


2.7  Dade  et  al.  (1992) 

It  was  the  purpose  of  this  work  to  define  the  critical  bed  shear  stress  required  to 
initiate  particle-by-particle  entrainment  of  cohesive  sediment.  The  investigators  start  their 
development  with  a  force  balance  similar  to  that  of  Wiberg  and  Smith  (1988).  The  only 
difference  is  the  inclusion  of  a  net  inter-particle  force,  FA,  which  includes  the  effect  of 
both  cohesion  and  adhesion  acting  in  the  same  direction  as  gravity.  Their  balance  is 

Wb  +  FA=FD  cotO  +  FL  (2.21) 

where  Wb  is  the  immersed  weight  of  the  grain,  and  FD  and  FL  are  the  drag  and  lift  forces 
due  to  the  overlying  flow.  This  equation  can  be  manipulated  to  take  the  form 
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D    _ 


tanO 


W  F, 

b     1  +  -ManO 


(       F  ^ 


(2.22) 


It  is  important  to  note  that  the  particle  packing  angle  <I>,  first  proposed  by  Wiberg  and 
Smith  (1988),  is  analogous  to  a  coefficient  of  static  friction  when  the  incipient  force 
balance  is  expressed  as  it  is  in  Equation  2.22.  The  modifications  then  proposed  by  the 
authors  were  an  explicit  inclusion  of  the  cohesion  term  and  appropriate  expressions  for 
FD  and  FL. 

The  investigators  assume  the  mud  bed  component  particles  to  be  spheres  that  are 
small  relative  to  the  thickness  of  the  near-bed  viscous  dominated  region  of  overlying 
turbulent  flows.  Consequently,  under  these  conditions  flat  beds  of  fine  sediment  are  said 
to  be  hydrodynamically  smooth  and  particle  Reynolds  number  is  around  unity.  Turbulent 
stress  then  makes  up  only  a  fraction  of  total  near-bed  stress.  The  resulting  near-viscous 
flow  exhibits  a  relatively  constant,  linear  velocity  profile.  Under  these  conditions  it  is 
appropriate  then  to  use  existing  expressions  for  FD  and  FL  acting  on  spheroidal  particles. 
These  expressions  are  substituted  into  Equation  2.22.  Ultimately  the  following  relation 
for  6cr,  the  Shields  parameter,  is  obtained: 


6,  =- 


-|2 


5X   VX 


-3  -1 


tanO       1  -  X 


(2.23) 


where 


~_g(ps-Pw)d3n 


Pwv 


(2.24) 


and 
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X  =  H 


1  +  -^- 
WbJ 


71  .3 


48006, 


tan'O  (2.25) 


where  b\  a  shape  factor  relating  viscous  forces  acting  on  actual  spherical  particles  to  that 
acting  on  nominal  spheres,  dn  is  the  nominal  spherical  particle  diameter,  and  S  is  referred 
to  as  the  Yalin  (1972)  parameter  which  is  a  dimensionless  submerged-particle  weight. 
Thus  the  critical  shear  stress  is  a  function  of  submerged  weight  and  shape  of  bed- 
component  particles,  relative  particle  protrusion  (embodied  in  tan<I>),  and  cohesion 
relative  to  submerged  particle  weight  (FA/Wt).  This  ratio  is  obtained  from 

-*-  =  3b2  (l  -  cos<I>)-^ — Ty     x  (2.26) 

where  fo  is  a  shape  factor  that  adjusts  the  ratio  of  particle  surface  area-to- volume  as  a 
function  of  platy-ness. 

The  strength  of  this  approach  lies  in  its  attempt  to  extend  non-cohesive  sediment 
transport  theory  to  mud  beds  by  including  an  additional  parameter  for  adhesion/cohesion 
force  FA  in  a  standard  force  balance.  Calculations  were  found  to  be  in  good  agreement 
with  the  observations  of  Migniot  (1968)  and  Ohtsubo  and  Muraoka  (1986). 

2.8  Nalluri  and  Alvarez  (1992) 

The  objective  of  this  work  was  to  identify  the  influence  of  cohesive  additives  to 
the  erosion  threshold  of  non-cohesive  sediments.  The  investigators  conducted 
experiments  using  both  a  pipe-flow  and  a  unidirectional  flume. 

They  started  with  a  determination  of  the  frictional  properties  of  both  flow 
systems.    Once  both  systems  were  calibrated,  they  prepared  sand  beds  in  both  systems 
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and  compared  the  mean  friction  coefficients  (fm)  with  the  bed  friction  coefficients  (/J,) 

computed  by  the  Einstein- Vanoni's  separation  technique  (Vanoni  and  Brooks,  1957). 

Tests  were  then  run  to  determine  the  influence  of  additives  such  as  china  clay,  oil,  and 

petroleum  jelly. 

Based  on  a  rheological  study  of  sewer  sediments,  the  investigators  used  the 

following  test  ranges: 

Concentration  of  laponite  clay  gel:  18  -  40  g/1 
Sand/clay  gel  ratio:  0  -  100%  by  weight 
Sand  sizes  (dso):  0.12  -  2.9  mm 

The  study  led  to  the  following  main  conclusions: 

1 .  Even  a  low  level  of  cohesion  can  increase  the  critical  shear  stress  significantly 
when  compared  with  equivalent  non-cohesive  sediments. 

2.  The  size  of  sand  has  no  significant  effect  on  the  critical  shear  stress  of 
cohesive  sediment. 

3.  For  a  given  sand  mixed  with  a  given  laponite  concentration,  there  was  an 
optimum  sand  to  clay  gel  ratio  for  which  a  maximum  resistance  to  erosion 
was  obtained. 

2.9  Mehta  and  Lee  (1994) 

In  this  study,  the  threshold  condition  for  the  incipient  transport  of  cohesionless 
grains  was  linked  to  cohesive  particles.  The  investigators  start  by  stating  the  well-known 
condition  for  incipient  transport 

where  FD  represents  the  hydrodynamic  drag  force  acting  on  a  grain,  Wb  is  the  buoyant 
weight  of  the  grain,  FL  is  the  hydrodynamic  lift  force  acting  on  a  grain,  and  Fc  is  the 
active  cohesive  force.  Therefore,  given  Tcr  is  equal  to  the  threshold  or  critical  bed  shear, 
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d  is  a  representative  grain  size,  a\  is  an  area  shape  factor,  cii  is  equal  to  «iC//Cd,  where 
d  and  Co  are  the  lift  and  drag  coefficients,  we  have  Fd  =  a\tCsd2,  and  Fl  =  ocitad2.  The 
buoyant  weight,  W  =  a3g(pSfpw)d  ,  where  a?,  is  a  volumetric  shape  factor,  g  is  the 
acceleration  due  to  gravity,  pSf  is  the  grain  (or  floe)  density,  and  pw  is  the  fluid  (water) 
density.  Substituting  for  these  quantities  in  Equation  2.27,  the  following  relationship  is 
obtained: 

rcr  a3tan(?       [  Fc  tan i<p (ax+a2  tan  g) 

g(psf-Pw)d     (a,+a2tanp)  g(psf  - pw)d3 

where  the  term  on  the  left  hand  side  is  Shields'  parameter,  \[/,  representing  the  ratio  of 
drag  force  to  buoyancy.  The  first  term  on  the  right  hand  side  is  seen  to  be  a  sediment- 
dependent  parameter,  and  the  second  term  is  essentially  the  ratio  of  cohesive  force  to 
buoyancy.  When  cohesion  is  unimportant  then,  this  term  is  essentially  zero  and  Equation 
2.28  becomes  a  representation  of  Shields'  entrainment  relationship  for  coarse  sediment. 

Further,  it  is  evident  from  Equation  2.28  that  Shields'  relationship  is  not 
applicable  to  situations  involving  flocculated  sediments,  and  there  is  further  implication 
that  extrapolation  of  this  relationship  into  the  cohesive  size  range  can  not  be  achieved 
unequivocally  without  including  some  measure  of  cohesion  as  a  behavior-governing 
parameter. 

The  investigators  point  out  that  an  evaluation  of  several  previous  laboratory 
studies  suggests  that  Ts  reasonably  correlates  with  the  solids  volume  fraction  0a.  In 
general, 
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*.»«4  (#.-#«/  (2-29) 

where  0flC  is  the  critical  solids  volume  fraction  such  that  for  all  <j>a  <  <j)ac  the  soil  does  not 
possess  a  measurable  structural  integrity  (even  though  it  does  have  cohesion),  and  is  thus 
not  a  bed  but  is  essentially  a  fluid-supported  slurry.  The  coefficients  0C4  and  /?  depend  on 
sediment-specific  properties. 

The  investigators  also  point  out  that  it  should  be  recognized:  (1)  for  cohesionless 
soils  the  threshold  condition  for  transport  is  specified  by  a  particular  value  of  the  applied 
stress  that  causes  the  grain  on  the  bed  surface  to  be  entrained,  and  (2)  the  threshold 
condition  for  a  cohesive  soil  bed  occurs  when  the  applied  stress  equals  the  shear  strength, 
associated  with  the  bed  surface  floes.  In  the  latter  case  the  floe  is  detached  from  its 
neighbors,  to  which  it  is  otherwise  linked  by  cohesion,  and  is  entrained. 

2.10  Mitchener  and  Torfs  (1996) 

The  paper  describes  the  results  on  the  erosion  of  sediment  mixtures  of  sand  and  mud 

based  on  laboratory  and  field  tests.   The  investigators  used  the  following  expression  to 

calculate  the  erosion  rate  for  each  shear  stress  step: 

dm  ,  x 

—  =  mc(tb-Tcr)  (2.30) 

where  dm/dt  is  the  rate  of  change  of  suspended  sediment  mass,  mc  is  the  erosion  constant, 
%  is  the  applied  bottom  shear  stress,  and  xcr  is  the  critical  shear  stress  for  erosion.  The 
following  are  the  main  conclusions  of  this  study: 

1.    The  critical  shear  stress  for  erosion  increases  when  mud  is  added  to  sand,  and 
consequently  the  onset  of  erosion  is  delayed.  The  authors  found  that  the  addition 
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of  30%  mud  to  a  sand  bed  can  increase  the  critical  shear  stress  by  as  much  as  a 
factor  of  10. 

2.  The  addition  of  mud  to  sand  diminishes  the  erosion  rates  once  sediment 
movement  does  occur. 

3.  The  results  in  this  work  agree  with  the  theory  of  Dade  and  No  well  (1991)  and 
indicate  that  the  transition  from  cohesionless  to  cohesive  erosion  occurs  between 
3%  to  15%  mud  by  weight.  Consequently,  cohesive  type  sediment  transport 
equations  should  be  used  where  the  sediment  bed  contains  more  than  15%  mud  by 
weight. 

2.11  Huygens  and  Verhoeven  (1996) 

It  was  the  aim  of  these  experiments  to  examine  similarities  and  differences  in  the 
behavior  of  mixtures  of  cohesionless  and  cohesive  sediment,  in  relation  to  the  well- 
established  cohesionless  sediment  transport  approach. 

Using  well-known  non-cohesive  sediment  transport  experiments  as  a  starting 
point  for  research  of  partly  cohesive  material,  the  investigators  conducted  a  series  of 
experiments  in  a  straight  semi-cylindrical  flume  which  was  1 1.2  m  long  and  had  an  inner 
diameter  of  390  mm.  The  inner  wall  was  covered  with  a  sand  paper  whose  wall 
roughness  was  evaluated  to  be  ks  =  0.54  mm. 

In  all  the  tests  sediment  mixtures  of  four  basic  materials  were  used:  two  non- 
cohesive  materials  (a  fine  white  sand  and  a  medium  beach  sand)  and  two  cohesive 
materials  (kaolinite  and  montmorillonite).  To  calculate  the  bed  shear  stress,  the 
investigators  used  the  Vanoni-Brooks  technique  for  removing  the  effect  of  the  side  walls 
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of  the  flume  (Huygens  and  Tito,  1994).  Using  this  technique,  it  was  found  that  the 
critical  shear  stresses  corresponding  to  the  incipient  motion  of  the  pure  sand  beds  was 
compatible  with  the  traditional  Shields  values.  The  critical  bed  shear  stress  of  the  partly 
cohesive  sediment  mixtures  varied  linearly  with  the  clay  content  up  to  a  clay  content  of 
30%.  For  sediment  mixtures  with  clay  content  greater  than  30%,  a  slight  decrease  of  the 
critical  bed  shear  was  noticed.  The  investigators  noted  that  the  cohesive  clay  particles 
cause  a  significant  increase  in  the  resistance  to  erosion. 

In  an  effort  to  quantify  the  bed  load  transport  of  mixtures,  the  investigators 
developed  an  expression  for  the  bed  load  transport  of  partly  cohesive  mixtures  that  was 
analogous  to  well-known  formula  for  cohesionless  sediment  previously  found  to  be 
successful  (Einstein,  1943;  Ackers-White,  1993;  Van  Rijn,  1984).  They  determined  the 
bed  load  transport  of  non-cohesive  and  partly  cohesive  sediment  to  be  related  to  the 
dimensionless  bed  shear  stress  rbt,  the  dimensionless  Chezy  coefficient  C6„,  and  the 

dimensionless  unit  stream  power  USf  I  w» ,  where   u„   is  the  critical  shear  velocity 

obtained  from  the  Shields  diagram  for  pure  sand  but  varies  with  clay  content  for 
mixtures,  U  is  the  mean  flow  velocity  and  S/  is  the  energy  slope.  They  found  a  general 
bed  load  transport  calculation  that  covered  all  mixtures  to  be 

/USA 


0,=O.OO2rJfcr 


V        *     J 


(2.31) 


where  the  dimensionless  bed  shear  stress  is  given  by 


*»  =    /     Xb    x,  (2.32) 
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and  the  dimensionless  Chezy  coefficient  is  given  by 
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Cb,  =  &  (2.33) 


2.12  Panagiotopolous  et  al.  (1997) 

The  objective  of  this  study  was  to  establish  the  critical  erosion  condition 
necessary  to  erode  a  sediment  mixture,  and  to  examine  the  influence  of  clay  content  on 
the  erosion  threshold. 

The  sediments  were  tested  in  a  recirculating  flume  (5.00  m  long,  0.30  m  wide  and 
0.45  m  deep)  with  an  open  top  and  glass-sided  walls.  A  video  camera  system  was  used  to 
observe  and  record  critical  erosion  condition.  Time-averaged  (2  min)  current  speeds  at 
threshold  were  obtained  0.4  cm  above  the  bed  using  a  laser-Doppler  velocimeter  operated 
at  a  sampling  frequency  of  20  Hz  and  measuring  velocities  to  within  plus  or  minus  0.2 
cm/s. 

The  sand  component  of  the  examined  soils  was  angular,  fine-grained  with  dso 
equal  to  153  and  215  Dm.  The  cohesive  component  was  Combwich  mud,  a  natural 
estuarine  silt-clay.  The  two  sands  were  mixed  separately  with  various  proportions  of 
Combwich  mud. 

Two  series  of  experiments  were  conducted  in  unidirectional  flow,  referred  to  as 
set  1  and  set  2.  In  set  1  the  erodibility  of  14  sediment  deposits  was  examined.  The 
threshold  condition  was  reached  within  2  min.  of  increasing  the  flow  from  zero. 
Following  each  test  the  water  content  of  the  mixture  was  determined.  In  set  2  the 
sediment  samples  were  first  subjected  to  various  pre-threshold  flow  speeds  (70%,  80%, 
and  90%  of  the  critical  values  defined  by  set  1),  over  different  times  (5  min,  10,  and  20 
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min)  and  then  tested  like  set  1.  These  experiments  focused  on  examining  the  erosion  of 
the  sediment  beds,  after  the  samples  were  exposed  to  current-induced  shear  stresses  at 
values  lower  than  the  critical  one. 

The  results  suggest  that  the  critical  stress  for  erosion  of  a  mixed  bed  increases 
with  increasing  mud  content.  Overall,  two  regimes  can  be  identified:  one  for  mud 
content  less  than  30%  (clay  mineral  content  10%),  and  another  for  greater  than  30%.  For 
mud  contents  <  30%,  under  a  steady  flow  a  slight  linear  increase  in  the  critical  stress 
occurs.  For  mud  contents  >  30%,  the  resistance  to  erodibility  increases  significantly, 
once  the  mud  content  exceeds  30%.  These  observations  agree  with  earlier  findings 
(Dyer,  1986;  Raudkivi,  1990)  that  a  clay  mineral  content  of  5-10%  is  high  enough  to 
dominate  sediment  erodibility. 

In  tests  on  the  erosion  of  the  mixtures  after  the  samples  were  exposed  to  current- 
induced  shear  stresses  at  values  lower  than  the  critical  one,  it  was  found  that  the  critical 
stress  increased  with  increased  pre-threshold  speed.  This  is  to  be  expected  because  the 
bed  becomes  more  compact  and  resistant  to  erosion  as  the  sub-threshold  velocity  is 
increased.  Grains  or  aggregates  that  do  not  satisfy  the  threshold  criterion  adjust  to  more 
stable  positions  within  the  bed. 

2.13  Torfs  et  al.  (2000) 

The  objective  of  this  work  was  to  a  develop  formula  for  the  critical  stress  for 
mixture  erosion,  rcm.  Following  the  work  of  Mehta  and  Lee  (1994)  reviewed  in  Section 
2.9,  the  investigators  noted  that  Equation  2.28  is  strictly  applicable  to  erosion  of  a  bed  of 
single-sized  grains,  even  as  it  allows  this  size  to  vary  from  fine  to  coarse.  When  fine  grain 


25 

aggregates  occur  in  a  mixture  with  coarse  grains,  the  interpretation  of  this  equation 
becomes  complex,  especially  because  the  two  particle  populations,  having  different  sizes, 
densities  and  cohesion  (equal  to  zero  for  coarse  grains),  do  not  necessarily  entrain  at  the 
same  rate. 

Therefore,  the  investigators  assume  that  the  entrainment  of  this  fine/coarse 
mixture  can  be  characterized  by  size,  dm,  and  granular  density,  psm,  the  latter  being  an 
indirect  measure  and  approximate  substitute  for  fine  sediment  aggregate  density.  For  the 
composite  (mixture)  bed,  the  corresponding  angle  of  repose  will  be<pm,  and  the  critical 

stress,  rcm.  Further,  for  fine  grains,  Fc  must  be  replaced  by  a  composite  force,  F„,  which 
depends  not  only  on  cohesion,  but  also  on  the  influence  of  the  inter-bedded  coarse  grains. 
Such  an  influence  can  either  enhance  the  effect  of  cohesion,  e.g.,  by  a  "hiding,"  or 
sheltering  effect,  or  reduce  cohesion  by  decreasing  inter-particle  bonding.  Thus,  the 
investigators  consider  Fm  =  KFC,  where  K  depends  on  bed  composition.  Hence,  Equation 
2.28  becomes 

rcm  _      a3  tan<pm        |  KF  tan ?„/(«,  +  a2  tan«pm ) 

g(Psm-Pw)dm     (a,+a2tan<pj  g{psm-pM 

in  which  cxi,  a2  and  a3  now  relate  to  dm.  For  a  uniform,  cohesive  bed  K=\,  and  the  first 

term  on  the  right  hand  side  loses  meaning  because,  as  noted,  q>m  is  not  definable  for  such 

a  bed.  Accordingly,  equating  zcm  with  the  bed  shear  strength  against  erosion,  r,  (Parchure 
and  Mehta,  1985),  from  Equation  2.34  we  obtain 

xs=-k-  (2.35) 

aag 
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where  K\  is  a  lumped,  sediment-specific  parameter  and  dag  is  the  aggregate  diameter.  An 
analysis  of  cohesive  bed  erosion  experiments  has  yielded  the  following  relation  (Mehta 
and  Parchure,  2000): 

t.  =  C<t  Hk  J6  (2-36) 

in  which  §vc  is  the  fine  solids  volume  fraction  obtained  by  dividing  the  dry  density  of 
solids,  i.e.,  sediment  dry  weight  per  unit  volume  of  fine  sediment- water  mixture,  by  the 
granular  density  of  the  fine  particles.  The  quantity  <t>v  is  the  "threshold"  value  of  <t>v  below 
which  the  bed  matrix  is  not  fully  particle-supported  and  is  therefore  fluid-like,  and  £  and 
%  are  coefficients  which  depend  on  bed  composition  and  the  degree  of  consolidation. 
From  Equations  2.35  and  2.36  we  obtain 

F=^k*>v-^  (2-37) 

Hence,  (2.34)  can  be  restated  as 

Tcm  a3tan<t>m 


g(psm-P)dm     (a1+a2tan0m) 


(2.38) 


ffi  gipsm-p)dV%    %c> 


Equation  2.38  can  be  simplified  by  recognizing  that  the  first  term  on  the  right  hand  side  is 
coarse  grain  dominated,  while  the  second  term  is  dominated  by  the  fine  constituent. 
Therefore,  in  the  first  term  the  subscript,  eg,  implies  coarse  grain  dominance.  As  for  the 
second  term  the  ratio,  Ktan^m/(a1  +a2tan<])m)#'1,  will  be  considered  to  be  a  bed 

dependent  characteristic  coefficient,  K' .  Then,  (2.38)  can  be  restated  as 
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T       = 

cm 


<*3cgtan(j)cg         |  K'Cfo0-*J 


g^sm-?)dm  (2.39) 


(alcg  +  a2cg  tan^)     g(psm-p)dr, 

where  the  authors  have  also  invoked  the  assumption,  dag  ~  dm.  In  (2.39),  K'=Q  implies 
the  absence  of  fine  material,  and  increasing  K'  signifies  increasing  modulation  of  the 
critical  stress  for  erosion  of  coarse  material  by  fines. 

Equation  2.39  is  significant  because  it  points  to  the  interactive  role  of  the  mixture 
components  in  controlling  erosion,  and  through  the  parameter  FC  suggests  a  potential 
approach  for  both  improving  the  modeling  of  erosion  of  sediment  mixtures  and 
explaining  the  observed  behavior.  Consequently,  the  authors  suggest  additional 
experimentation  on  the  erosion  of  mixtures  in  the  vicinity  of  space-filling  volume 
fraction. 

2.14  Sharif  (2002) 

The  objective  of  this  work  was  to  develop  a  model  for  the  critical  shear  stress  for 
the  erosion  of  cohesive  soils.  The  author's  model  included  a  relation  for  critical  shear 
stress  of  pure  cohesive  soils,  a  relation  for  critical  shear  stress  of  mixed  cohesive/non- 
cohesive  soils  with  low  clay  content,  and  a  relation  for  critical  shear  stress  of  mixed 
cohesive/non-cohesive  soils  with  high  clay  content. 

The  relation  for  critical  shear  stress  of  mixed  cohesive/non-cohesive  soils  with 
low  clay  content  was  based  on  a  force  balance  on  a  single  non-cohesive  particle  similar  to 
Equation  2.27  and  resulted  in  the  following  expression  for  Tcn,  the  critical  shear  stress  for 
a  non-cohesive  particle  in  a  cohesive  bed 


Ta,=8(pm-P)dn0  +  - 
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^ccjnd        (    fcPdPs 


.   yjr  sn         J  nr^d  /  . 


(2.  40) 


a  a  nd  d2p2p 

vc     vn         n     c  f     f  s 

where  ps„  is  the  granular  density  of  the  non-cohesive  sediment  grain,  psc  is  the  granular 
density  of  the  fine  cohesive  particles,  dnc  is  the  diameter  of  the  non-cohesive  particle,  fc  is 
the  weight  fraction  of  fine  cohesive  particles,  /„  is  the  weight  fraction  of  non-cohesive 
particles,  Pd  is  the  bulk  density  of  the  bed,  £„  is  the  magnitude  of  the  force  associated 
with  the  bonds  between  a  fine  cohesive  particle  and  a  non-cohesive  particle,  6  is  the 
Shields  parameter,  dc  is  the  average  diameter  of  fine  cohesive  particles,  cCvC  and  a„„  are 
volumetric  shape  factors  for  fine  cohesive  particles  and  non-cohesive  particles 
respectively,  and  ccs„  is  as  area  shape  factor  for  non-cohesive  particles. 

The  relation  for  the  critical  shear  stress  for  erosion  of  pure  cohesive  soils  was 
based  on  a  force  balance  for  the  separation  of  the  aggregates  at  the  surface  of  the  bed 
from  the  bed  itself.  This  force  balance  ultimately  resulted  in  the  following  expression  for 
Tea,  the  critical  shear  stress  for  a  pure  cohesive  bed,  as  a  function  of  aggregate  density 

*i*0.         ,  18P„% 


7 '-^y^T+         V  (2-41) 


where  pag  is  the  granular  density  of  the  aggregates,  §,  is  the  magnitude  of  the  average 
force  associated  with  the  bonds  between  two  fine  cohesive  particles,  6a  is  the  Shields 
parameter  for  cohesive  aggregates  estimated  from  the  Shields  diagram,  and  k\  and  k2  are 
coefficients  from  the  following  relation  for  floe  diameter,  da,  as  a  function  of  floe  density 
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The  relation  for  critical  shear  stress  for  erosion  of  mixed  cohesive/non-cohesive 
soils  with  high  clay  content  was  based  on  the  presumption  that  the  critical  shear  is 
governed  by  the  critical  shear  stress  of  the  non-cohesive  particles  and  the  critical  shear 
stress  of  the  aggregates,  and  the  degree  by  which  these  two  different  critical  shear 
stresses  influence  the  bed  shear  stress  overall.  The  author  made  the  assumption  that  the 
degree  of  impact  of  each  critical  stress  type  depended  on  the  average  area  at  the  plane  of 
failure  covered  by  each  type  of  material  (cohesive  or  non-cohesive).  Using  this 
assumption,  the  author  arrive  at  the  following  expression  for  zCh,  the  overall  critical  shear 
stress  of  the  mixed  cohesive/non-cohesive  beds  with  high  clay  content 

*c=*cnWn+TcaWa  (2.43) 

where  W„  and  Wa  are  weighting  factors  representing  the  contribution  of  the  critical  shear 
stress  of  non-cohesive  particles  and  cohesive  aggregates,  respectively,  to  the  overall 
critical  shear  stress  of  the  bed.    Wn  and  Wa  are  defined  as  follows 

W„  =  Abni  (2.44) 

Wa  =  Abai  (2.45) 

where  A„b  is  the  specific  area  of  the  non-cohesive  particles,  Aab  is  the  specific  area  of  the 
aggregates,  and  b„  and  ba  are  experimentally  determined  coefficients. 

After  applying  the  model  for  low  clay  concentrations  and  analyzing  the  results, 
the  author  made  the  following  important  conclusion:  the  critical  shear  stress  for  erosion 
of  mixed  cohesive/non-cohesive  soils  with  low  clay  content  is  influenced  solely  by  the 
forces  acting  on  the  non-cohesive  particles  and  that  this  was  the  case  up  to  the  clay 
content  associated  with  the  maximum  critical  shear  stress. 


CHAPTER  3 
CRITICAL  STRESS  FOR  SAND  BED  WITH  INTERSTITIAL  CLAY  PARTICLES 


3.1  Introduction 

An  expression  for  the  critical  stress  of  sand  in  the  presence  of  clay  particles  in 
pore  water  is  developed  in  this  chapter.  To  begin  with,  the  incipient  motion  for 
cohesionless  grains  is  defined  and  an  expression  for  the  critical  stress  of  cohesionless 
sediment  is  developed. 

3.2  Incipient  Motion  of  Cohesionless  Particles 

Shields  (1936)  defined  the  critical  condition  as  the  point  at  which  zero  sediment 
transport  occurs.  In  other  words,  the  critical  condition  is  the  point  of  maximum  stress  at 
which  no  sediment  transport  occurs,  therefore  allowing  for  weak  sediment  movement  but 
no  mass  transport  of  particles.  Kramer's  (1935)  definition  of  weak  movement  seems  to 
correspond  to  Shields'  critical  condition,  where  it  is  understood  that  with  weak 
movement  a  few  of  the  smallest  sand  particles  are  in  motion  in  isolated  spots.  Vanoni 
(1946)  advanced  a  criterion  for  incipient  motion  by  counting  the  number  of  bursts  per 
second;  critical  is  defined  with  the  burst  frequency  of  1/3  to  1  bursts  per  second. 
Sutherland  (1967)  found  that  near  the  threshold  on  a  flat  bed,  motion  was  restricted  to 
isolated  spots  with  a  small  number  of  grains  moving  briefly  a  few  centimeters.  Bursts  of 
motion  occurred  from  positions  all  over  the  bed  with  frequencies  that  appeared  to  be 
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constant  from  one  location  to  another.   Sutherland  connected  the  intermittent  movement 
with  penetration  of  eddies  into  the  viscous  sublayer. 

Given  the  above  observations  it  is  necessary  to  explain  the  formal  use  of 
"incipient  motion"  in  relation  to  "incipient  transport."  From  the  work  on  incipient 
motion  done  to  date,  the  first  term  is  most  frequently  used  to  describe  the  instantaneous 
state  at  which  a  solid  moves.  This  definition,  while  concise,  may  still  lead  to  ambiguities 
between  the  work  of  one  researcher  and  the  work  of  another.  For  example,  it  does  not 
follow  that  once  surface  solids  have  moved  in  a  particular  area  of  the  bed,  that  they  will 
continue  to  do  so,  since  the  condition  of  "armoring"  may  occur  (Mantz,  1977).  It  also 
does  not  follow  that  when  an  individual  grain  or  small  grouping  of  surface  solids  is 
moving  at  one  location,  that  others  will  be  moving  at  a  different  location,  due  to  the 
stochastic  nature  of  the  properties  of  both  the  flow  and  the  bed.  A  time  and  spatially 
averaged  critical  stress  which  determines  the  fluid  condition  for  incipient  motion  need 
not,  therefore,  be  that  which  describes  the  initial  continuous  motion  of  solids,  or  incipient 
transport. 

Defining  incipient  motion  using  average  flow  properties  tends  to  over-simplify 
the  entrainment  process  in  light  of  the  stochastic  nature  of  turbulent  flow.  In  that  regard, 
Christensen  (1972)  and  others  have  shown  that,  due  to  the  randomness  of  the  velocity 
fluctuations,  the  instantaneous  bed  shear  stress,  tb,  may  be  related  to  a  time-mean  value 
of  the  shear  stress,  fcr,  when  the  probability  of  erosion,  i.e.,  the  probability  that  the 
instantaneous  shear  stress  exceeds  fcr,  is  prescribed.  Christensen  (1972)  proposed  the 
following  equation: 


7cr 


u  +u 
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v     "      y 


(3.2.1) 


1  +  ^-n 


in  which  w '  is  the  instantaneous  value  of  the  velocity  fluctuation,  u  is  the  average  value 
of  the  velocity,  au  is  the  standard  deviation  of  u ',  and  n  =  u  7ctu.  Using  Einstein  and  El 
Samni's  (1949)  observations  of  the  fluctuating  lift,  Christensen  (1965)  demonstrated  that 
the  velocity  fluctuations  are  normally  distributed  near  the  bed  and  that  ou/u  must  be 
constant  and  equal  to  0.16,  i.e.,  a  value  that  is  in  good  agreement  with  direct  observations 
(e.g.,  Reichardt,  1938).  Consequently,  n  may  be  found  from  the  table  of  the  Gaussian 
probability  integral  when  the  probability  of  exceeding  %,  and  thereby  risking  erosion,  is 
decided.  Christensen  and  Bush  (1971)  showed  that  good  agreement  with  observations  of 
the  critical  shear  stress,  fcr ,  may  be  obtained  if  this  probability  is  0.001  corresponding  to 
an  n  value  of  3.09.  Hence,  since  sediment  erosion  rate  tends  to  increase  rapidly  with  the 
average  bed  shear  stress  in  slight-to-moderate  excess  of  calculated  critical  values,  the 
definition  of  a  critical  shear  stress  using  time-averaged  flow  properties  (fcr)  remains 
useful  for  sediment  transport  modeling. 

3.3  Critical  Stress  of  Cohesionless  Sediment  with  Clay  Particles  in  the  Pore  Water 

Consider  a  planar  stationary  bed  of  cohesionless  particles  subjected  to 
unidirectional  turbulent  flow.  It  is  assumed  the  flow  and  bed  are  statistically  in  a  steady 
state.  As  soon  as  the  fluid  starts  to  flow,  hydrodynamic  forces  are  exerted  on  the 
particles,  and  with  an  increase  of  the  flow  velocity  a  condition  is  eventually  reached 
where  the  particles  at  the  surface  of  the  bed  are  unable  to  resist  the  forces  and,  as  a  result, 
the  first  grain  get  dislodged  and  moves.  This  condition  of  incipient  motion  is  defined  by 
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tanp  =  -5-  (3.3.1) 

where  F,  and  Fn  are  the  instantaneous  forces  parallel  and  normal  to  the  angle  of  repose,  j. 
F,  and  Fn  are  resultants  of  the  hydrodynamic  drag  FD,  the  lift  force  Fl,  and  the  submerged 
weight  Wb.  Substituting  in  these  forces,  Equation  3.3.1  becomes 

tan<p  =  -^ L_L_£l  (3.3.2) 

Wbcosdb-FL 

where  the  angle  Ob  is  the  inclination  of  the  bed  from  the  horizontal  at  which  incipient 

sediment  motion  takes  place.  Consequently, 

,     .^(taiKDcosg.-sineJ 

\rD\- 7; {5.5.5) 

1  +  — tan<I> 

FD 

where  the  particle  packing  angle  (or  bed  pocket  geometry  angle),  O,  is  substituted  for  the 
angle  of  repose,  q>,  following  Wiberg  and  Smith  (1987).  Figure  3.1  illustrates  the  force 
balance.  The  bed  pocket  geometry  angle  is  a  measure  of  the  difficulty  of  moving  a 
sediment  grain  up  and  out  of  the  pocket  in  which  it  is  resting.  The  instantaneous  drag 
force,  the  instantaneous  lift  force,  and  the  buoyant  weight  are  all  expressed  as  follows: 

FL=Xaxd2  (3.3.4) 

FD=xbaxd2  (3.3.5) 

Wb=a2(Ys-Y)d3  (3.3.6) 

where  Oil  is  a  shape  factor  representing  the  cross-sectional  area  of  the  grain,  0C2  is  a  shape 
factor  representing  the  volume  of  the  grain,  A  is  the  value  of  the  instantaneous  lift  force 
per  unit  area,  %  is  the  value  of  the  instantaneous  shear  stress,  ys  is  the  unit  weight  of  the 
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Figure  3.1  Force  balance  on  a  particle  at  the  surface  of  a  bed. 


sediment,  y  is  the  unit  weight  of  water,  and  d  is  the  grain  size.  Substituting  into  Equation 
3.3.3  gives 


I     i_  a2  (Ys  -y)d(tanOcos0A  -s'mdb) 


a, 


1  +  — tanO 

T, 


(3.3.7) 


This  is  the  instantaneous  bed  shear  stress  at  which  failure  is  about  to  occur  for  surface 
particles  on  a  bed  with  a  sloping  angle  of  $,.  Following  Christensen  (1972)  we  will 
define  Tcr  as 


K\  =  cdPy  =  cd  p^u  +  U^  =  c°  P\ 


1  +  ZjlJL 


u    a \. 


(3.3.8) 


where  u  is  the  time-averaged  velocity,  w'is  the  instantaneous  fluctuation  of  velocity, 
and  au  is  the  standard  deviation  of  the  instantaneous  velocity  fluctuations.  Consequently, 
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(3.3.10) 


where  5W  =cth/m  and  n  =  u'/cru .  Therefore, 


T      =T 

cr  cr 


,     (1  +  ^n)2 


1  +  5,: 


(3.3.11) 


Equivalently, 


1  /_ 


\l\  =  CLp^  =  CLp±(u+u'f=CLp^ 
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(3.3.12) 
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=  CLPy(l  +  5„2)  (3.3.13) 
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(3.3.14) 


and  therefore, 


-C^"H^.!) 


Substitution  back  into  Equation  3.3.7,  results  in 

_    (l  +  Sun)2  _a2  (ys  -y)d(tanOcos0fc  -sin0fc) 


1  +  5,: 


a, 


1  + 


CL 


f77\2 


ku.j 


(l  +  S2)tanO 


(3.3.15) 


(3.3.16) 


which  simplifies  to 


36 


-    _«2  (y,  -  y)d (tanO cos0 b  -  sin 6b  jl  +  S2U ) 


T„  = 


a, 


1  + 
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(  tt\ 


\U*J 


(l  +  S„2)tanO 


(3.3.17) 


(l  +  Sunf 


At  this  point,  an  expression  for  the  turbulent  velocity  profile,  w/w, .  For  a 
turbulent  flow  the  vertical  profile  of  velocity  for  the  near-bed  flow  is  a  function  of  the 
roughness  Reynolds  Number,  /?,  =  utks  /v ,  where  ks  is  the  length  scale  for  the  bed 

roughness.    When  the  flow  is  in  the  hydraulically  rough  range,/?,  MOO,  the  flow  is 
found  to  increase  logarithmically  away  from  the  bed  for  z  >  3^: 


w. 


u(z)  =  ^\n 
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(3.3.18) 


where  the  von  Karman's  constant  is  k  =  0.407,  and  za  =  ks  /30.  In  fully  rough  flow,  the 
viscous  sublayer,  &,  is  compressed  compared  to  hydraulically  smooth  or  transitional 
flow,  and  the  individual  grains  are  found  to  extend  up  through  the  sublayer  and  into  the 
turbulent  flow,  resulting  in  flow  separation  around  the  individual  grains  and  a  consequent 
independence  of  critical  shear  (Shields  parameter)  from  Reynolds  number.  When  the 
flow  is  in  the  hydraulically  smooth  range,  R.  >  3 ,  the  grains  at  the  surface  of  the  bed  are 
considered  to  be  totally  submerged  in  the  viscous  sublayer,  which  has  an  average  height 
of  6V  =11.6v/m,  .     Within  the  viscous  sublayer,  the  velocity  profile  is  linear,  and 

for/?,  > 3,  u(z)=  u]  zjv  in  this  region.  Above  the  viscous  sublayer,  the  velocity  profile 

is  logarithmic  and  obeys  Equation  3.3.18  with  z0  =v/(9«.)  for  hydraulically  smooth 

flow. 

Wiberg  and  Smith  (1987)  use  the  equation  of  Reichardt  (1951)  for  the  velocity 
profile  for  transitional  flow  which,  they  argue,  provides  a  nice  transition  between  the 
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viscous  sub-layer  and  the  logarithmic  flow  above  and  is  also  in  excellent  agreement  with 
measurements,  and  use  Equation  3.3.18  for  hydraulically  rough  flow.  One  of  the  ideas 
advanced  in  their  work  is  that  in  a  heterogeneous  mixture  of  cohesionless  grains,  the 
grains  that  are  most  likely  to  move  are  those  that  sit  relatively  high  on  the  bed  (e.g.,  those 
with  their  bottoms  near  the  zero  velocity  level).  They  argue  that  these  particles  are 
minimally  affected  by  the  wakes  of  upstream  particles  when  the  flow  is  in  the  fully  rough 
region  of  flow,  and  that  it  was  therefore  acceptable  to  extend  the  logarithmic  velocity 
profile  down  to  z  =  z0-  In  the  hydraulically  smooth  region  of  flow,  they  agree  that  the 
aforementioned  expression  for  flow  inside  the  viscous  sublayer  was  only  valid  for  z  >  ks, 
but  argue  that  since  ks  «  &,  it  was  acceptable  to  approximate  the  near-bed  flow  by 
extending  the  linear  velocity  profile  to  the  bed.  They  also  argue  that  the  transition 
between  the  viscous  sublayer  and  the  logarithmic  flow  above  is  not  abrupt  because  the 
turbulent  eddies  depress  and  expand  the  viscous  sublayer  as  they  pass,  and  therefore  the 
velocity  at  z  =  d  in  reality  lies  between  the  linear  and  logarithmic  values.  Consequently, 
they  sought  an  expression  for  velocity  that  yields  a  linear  profile  for  z  «  &,  a 
logarithmic  profile  for  z  »  &,  and  provides  a  smooth  transition  in  the  region  in  between. 
For  this  reason  they  use  the  expression  of  Reichardt  (1951) 


u(z)=u. 


-ln(l  +  *z+)-, 


(  ,+  \ 

l_e-zt-M-6  Z       c-033z+ 


/J 


(3.3.19) 


11.6 

where  z+  =  u,  z/v  =  R,  z/ks  and  z+0  =  «.  zjv  =  /?,  zjks ,  to  determine  the  velocity  in  the 
transitional  region  of  flow  (3  <  Rt  <  100).  In  the  transitional  region,  the  bed  roughness 
parameter,  za,  is  a  function  of  roughness  Reynolds  number  as  given  by  Nikuradse  (1933). 
Figure  3.2  shows  the  vertical  velocity  profile  calculated  from  Equation  3.3.19  for/?,  =  1 . 
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Figure  3.2  Velocity  as  a  function  of  non-dimensional  height  above  the  bed.  Equation  of 
Reichardt(1951). 


It  should  be  noted  that  with  Equation  3.3.19,  all  of  the  terms  in  Equation  3.3.17 
can  be  either  measured  or  calculated  with  the  exception  of  the  particle  packing  angle,  O. 
It  should  also  be  mentioned  that  to  use  these  two  Equations  (3.3.18  and  3.3.19)  for  non- 
dimensional  velocity  with  the  expression  for  critical  shear  stress,  Equation  3.3.17,  it  is 
necessary  to  use  Equation  3.3.19  in  the  transitional  region  of  flow  and  to  use  Equation 
3.3.18  in  the  hydraulically  rough  region  of  flow.  It  is  also  necessary  to  select  the 
appropriate  elevation  above  the  bed  at  which  to  determine  the  non-dimensional  velocity 
(this  elevation  is  discussed  subsequently). 
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Christensen  (1972)  notes  that  the  classical  expression  based  on  Prandtl's  well- 
known  mixing  length  hypothesis  gives  very  good  results  at  moderate  to  great  distances 
from  the  bed,  but  that  it  fails  near  the  bed  where  z  approaching  zero  results  in  a  time- 
mean  velocity  approaching  the  highly  unrealistic  value  of  minus  infinity.  Consequently, 
he  used 


u         .  (29.6z 


—  =  -ln 
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(3.3.20) 


This  velocity  profile  is  identical  to  the  classical  expression  except  for  the  +1  term  which 
makes  the  time-mean  velocity  equal  to  0  at  the  bed,  and  is  based  on  the  work  of  Rossby 
(1932).  This  profile  is  plotted  in  Figure  3.3.  Since  it  is  only  valid  for  the  hydraulically 
rough  range,  it  can  be  modified  it  so  that  it  may  be  used  in  both  the  transitional  region  of 
flow  and  the  fully  rough  region: 


JL  -I 

U,        K 


(      d  1 

q—e^+l 


(3.3.21) 


where  q  is  the  elevation  above  the  bed  expressed  as  a  fraction  of  the  grain  size,  d,  and  B 
is  the  Nikuradse  velocity  coefficient.  It  should  be  noted  that  as  with  the  Reichardt  (1951) 
formula,  in  this  expression  all  of  the  terms  in  Equation  3.3.17  can  be  either  measured  or 
calculated  with  the  exception  of  the  particle  packing  angle,  O.  As  before,  it  is  necessary 
to  select  the  appropriate  elevation  above  the  bed  at  which  to  determine  the  non- 
dimensional  velocity. 

Guoren  (1981)  developed  an  expression  that  accurately  describes  the  velocity,  at  a 
given  point  above  the  bed,  throughout  all  three  regions  of  flow.  He  considered  the 
occurrence  of  the  three  regions  to  be  due  to  the  degree  of  separation  of  the  flow  passing 
over  the  roughness  elements.  He  argued  that  in  the  case  of  homogeneous  roughness,  the 
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roughness  element  can  be  considered  approximately  as  a  sphere  and  that,  consequently, 
the  degree  of  separation,  as,  may  be  expressed  as  follows: 
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Figure  3.3  Turbulent  time  mean  velocity  distribution  in  proximity  of  rough  bed,  rough 
range.  Profile  of  Christensen  (1972). 


a,=(l-cos0,J/2 


(3.3.22) 


where  6sep  is  the  separation  angle.    Based  on  this  concept  he  developed  the  following 
generalized  velocity  profile: 
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where 
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and 


P,=l-(l-P0Xa,+eg/n)/2 


6S  /n  =  (in  Re  -  In  Re,  )/(ln  Re2  -  In  Re, ) 


(3.3.25) 
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in  which  Rei  =  1.25,  Re2  =  100,  and  fi0  =  0.107.  A  comparison  of  this  profile  with  the 
other  three  and  the  data  of  Nikuradse  (1933)  is  shown  in  Figure  3.4  from  which  it  may  be 
seen  that  a  good  agreement  is  obtained. 

Coleman  (1981)  developed  an  expression  that  accurately  describes  the  velocity 
above  the  bed  of  a  sediment  laden  fluid.  He  accounted  for  the  effect  of  suspended 
sediment  through  a  modified  expression  for  kinematic  viscosity  and  a  wake  region 
augmentation  function.  His  expression  for  viscosity  is 


v  = 


_n(l  +  2.5(p  +  625(p2w+l5.62<pl) 
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(3.3.28) 


where  \i  is  the  dynamic  viscosity  of  the  fluid,  (pw  is  the  sediment  concentration  at  the 
point  in  question  expressed  as  volume  of  sediment  per  unit  volume  of  sediment-fluid 
mixture,  ps  is  the  sediment  density,  and  p  is  the  fluid  density.  His  velocity  profile  for 
bounded  shear  flows  is 
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where  A  is  an  integration  constant,  Aw/w,  is  the  channel  roughness  velocity  reduction 
function  (Hinze,  1951),  and  (jJ/K)m(z/d)  is  the  wake  region  velocity  augmentation 
function  containing  the  wake  strength  coefficient  I"I,  and  the  boundary  layer  thickness  <5. 
The  symbol  co  is  merely  a  functional  symbol.  The  augmentation  function  is  defined  as 
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where  um  is  the  maximum  velocity.  A  comparison  of  this  profile  with  the  other  three,  as 
they  vary  over  the  transitional  and  hydraulically  rough  ranges  of  flow,  and  with  the  data 
of  Nikuradse  (1933)  is  shown  in  Figure  3.4.  All  four  profiles  seem  to  compare  well  with 
each  other  and  with  the  data. 
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Comparison  of  the  Nondimensional  Velocity  Profiles 
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Figure  3.4  Comparison  of  four  velocity  profiles  with  the  data  of  Nikuradse  (1933). 
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The  above  four  profiles  were  then  substituted  into  the  expression  for  critical 
bottom  shear,  Equation  3.3.17,  and  were  found  to  compare  well  with  the  data  of  Shields 
(1936)  as  can  be  seen  in  Figure  3.5. 

Values  for  several  of  the  parameters  in  Equation  3.3.17  need  to  be  defined  before 
it  may  be  used: 

a, :  This  parameter  is  the  shape  factor  representing  the  bed  area  associated  with  one 
grain  divided  by  the  grain  size  squared.  If  the  grains  are  spherical  and  arranged  in  a 
pattern  that  may  be  described  between  a  square  and  a  hexagon,  it  can  be  shown  that  al  = 
0.933. 
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Figure  3.5  Comparison  of  profiles  for  critical  bottom  shear  with  the  data  of  Shields 
(1936). 
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a2 :  This  parameter  represents  the  grain  volume  divided  by  the  grain-size  cubed.  By 
definition,  this  parameter  is  equal  to  p/6  or  0.524. 

Su  and  n:  Using  Einstein  and  El  Samni's  (1949)  observations  of  the  fluctuating  lift, 
Christensen  (1965)  showed  that  the  velocity  fluctuations  are  normally  distributed  near  the 
bed  and  Su  is  constant  and  equal  to  0.164,  a  value  that  is  in  good  agreement  with 
observations.  Consequently,  n  may  be  found  from  the  table  of  the  Gaussian  probability 
when  the  probability  of  exceeding  tcr  is  selected.  As  noted  earlier,  Christensen  and  Bush 
(1971)  showed  that  good  agreement  with  observations  was  obtained  if  this  probability  is 
0.001  corresponding  to  n  =  3.09. 

CL:  Einstein  and  El  Samni  (1949)  and  Christensen  (1965)  found  a  value  of  CL  =  0.178 
produced  the  closest  agreement  with  measurements,  when  the  velocity  used  was 
measured  at  a  distance  of  0.35  dso  above  the  theoretical  bed. 

O:  Wiberg  and  Smith  (1987)  argue  that  the  particle  packing  angle  changes  for  a  bed 
composed  of  mixed  sizes  because  it  is  easier  for  a  large  grain  to  roll  over  a  bed  of  smaller 
particles  than  for  a  small  grain  to  roll  over  a  bed  of  larger  particles.  This  is  indicated  by 
an  increase  in  the  value  of  <I>  as  the  ratio  of  grain  diameter  to  bed  roughness  length,  d/ks, 
decreases  and  is  illustrated  in  Figure  3.6.  <J>  was  defined  using  an  expression  derived 
from  the  work  of  Miller  and  Byrne  (1966): 

~  ds 


<J>  =  cos" 


(3.3.32) 


45 

where  z*  is  the  average  level  of  the  bottom  of  the  almost  moving  grain.  The  best 
agreement  with  Miller  and  Byrne  (1966)  was  obtained  with  zt  =-0.02  for  natural  sand 
having  a  sphericity  of  0.7  and  roundness  of  0.5. 
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Figure  3.6  Relationship  between  the  particle  packing  angle  and  the  ratio  of  particle  size  to 
bed  roughness  scale,  D/ks,  where  D  is  the  grain  diameter.  In  this  representation  the 
bottoms  of  the  almost-moving  grains  are  at  the  zero  level  of  the  bed,  indicated  by  arrows 
at  z  =  0  (Wiberg  and  Smith,  1987). 


Dade  and  No  well  (1991)  suggested  a  threshold  value  of  the  solids  fraction  of  a 
clay  suspension  above  which  a  networked  structure  is  formed  at  which  point  the  sediment 
mixture  behaves  primarily  as  if  it  were  a  mud.  They  suggested  threshold  values  of  solids 
fraction  y  between  0.01  and  0.05  (1  and  5  percent)  where  vj/  is  defined  as 
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r1 clay 


(3.3.33) 


^clay  Pclay  +  ® sand  P  sand 

The  results  of  Torfs  (1995)  concerning  the  mode  of  erosion  agreed  with  this  concept  and 
indicated  a  transition  from  sandy  to  muddy-type  erosion  between  3%  and  15%.  The 
results  of  Nalluri  and  Alvarez  (1992)  indicated  a  transition  between  2%  and  17%. 
Mitchener  and  Torfs  (1996)  mention  a  transition  zone  in  the  mode  of  erosion  from  sand- 
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type  behavior  to  mud-type  behavior  between  7%  and  13%  from  their  results  with 
montmorillonite  mixtures.  The  results  of  Huygens  and  Verhoeven  (1996)  for  kaolinite 
mixtures  show  a  transition  between  2%  and  9%.  Panagiotopoulos  et  al.  (1997)  noted  two 
different  regimes  of  critical  conditions  for  erosion;  one  at  mud  content  less  than  30%,  and 
one  at  mud  content  greater  than  30%.  The  results  of  Sharif  (2002)  show  a  transition  zone 
between  5%  and  15%  for  kaolinite  and  sand  mixtures.  Examples  of  this  trend  are  given 
in  Figures  3.7  and  3.8. 

The  non-linear  dependence  of  the  critical  shear  stress  on  the  fine-grained  weight 
fraction  is  indicative  of  the  interactive  role  of  the  sand  grain  matrix  and  the  interstitial 
clay  particles  in  controlling  critical  shear  stress,  and  further  suggests  a  potential  route  for 
improving  this  relationship.  Incipient  motion  of  an  individual  sand  grain  with  interstitial 
clay  particles  can  be  examined  by  considering  the  forces  acting  on  the  grain  resting  at  the 
bed  surface  subjected  to  turbulent  flow. 

As  clay  is  added  to  sand  the  clay  particles  start  to  fill  in  the  pore  spaces  of  the 
sand  matrix  altering  the  pivoting  characteristics  and  consequently  the  value  of  the  internal 
friction  angle.  At  low  clay  content,  clay  apparently  acts  to  lubricate  the  sand  grain  at  its 
contact  points  reducing  the  internal  friction  angle  and  resulting  in  the  decrease  in  f  as 
shown  in  Figures  3.7  and  3.8.  Hence,  this  lubricating  effect  of  clay  can  be  shown  through 
the  adoption,  in  Equation  3.3.17,  of  an  effective  particle  packing  angle,  Oeff: 
_    _a2(yI-r)d  (tan Oeff  cos0fc  -  sin 6b  jl  +  S2U ) 
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Figure  3.7  Critical  shear  stress  for  erosion  of  kaolinit/sand  and  Scheldt  mud/sand  mixtures 
versus  fine  grain  weight  fraction:  data  (Torfs,  1995). 
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Figure  3.8  Critical  shear  stress  for  erosion  of  Montmorillonite/sand  mixture  versus  fine 
grain  weight  fraction:  data  (Torfs,  1995) 


CHAPTER  4 
EXPERIMENTAL  SETUP  AND  PROCEDURES 


4.1  Introduction 

To  test  the  validity  of  the  analysis  carried  out  in  Chapter  3,  which  derives  the 
critical  shear  stress  of  sand  in  the  presence  of  clay,  laboratory  experiments  were 
conducted.  The  experimental  setup,  procedures  and  tests  are  presented  in  this  chapter. 

4.2  Flume 

Tests  were  carried  out  in  a  small  flume  at  the  Coastal  Engineering  Laboratory  of 
the  University  of  Florida.  McAnally  (2000)  has  provided  a  brief  summary  of  the 
underlying  design-and-operation  concept  of  such  open  channel  flumes.  Typically,  a 
pump  is  used  to  deliver  water  from  a  tail  basin  to  a  head  basin.  The  water  then  flows 
through  the  flume  back  into  the  tail  basin  completing  the  circuit.  The  velocity,  and 
consequently  the  bottom  shear  stress,  are  controlled  by  regulating  the  pumps'  discharge. 

Figure  4.1,  which  shows  the  schematic  of  the  flume,  was  constructed  of  3 -cm- 
thick  Plexiglas  to  dimensions  of  4.3  m  length,  15  cm  width,  and  19  cm  depth.  It  was 
elevated  to  a  height  of  92  cm  above  the  floor  and  supported  by  a  wooden  table  the  base  of 
which  was  adjustable  and  used  to  level  the  flume. 

The  head  basin  was  constructed  of  3-cm-thick  Plexiglas  to  dimensions  of  89  cm 
length,  183  cm  width,  and  29  cm  depth.  It  was  fed  by  a  Worthington  model  D814 
centrifugal  pump  with  a  Baldor,  single-phase,  1  hp,  L5510F  motor.   The  basin  was  fed 
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via  a  7.6  cm  OD  PVC  pipe  into  which  several  holes  were  drilled  so  as  to  evenly  distribute 
the  incoming  water  and  minimize  any  currents  in  the  head  basin.  The  current  velocity  in 
the  input  pipe  was  measured  using  a  Signet  model  MK525  metalex  "paddle-wheel"  type 
current  meter. 
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Figure  4.1  Schematic  drawing  of  flume. 
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A  vertical  sluice  gate  was  installed  at  the  end  of  the  flume  to  regulate  the  water 
surface  height.  It  was  constructed  of  6.6  mm  Plexiglas  slightly  wider  than  the  width  of 
the  flume,  and  was  placed  at  the  end  of  the  flume  immediately  before  the  tail  basin.  The 
sluice  gate  could  be  raised  or  lowered,  which  would  in  turn  increase  or  decrease, 
respectively,  the  flow  area  underneath  the  gate.  Consequently,  for  a  given  discharge  the 
sluice  gate  could  be  used  to  control  the  surface  height  of  the  water  in  the  flume. 


50 


The  tail  basin  was  constructed  of  3-cm-thick  Plexiglas  to  dimensions  of  120  cm 
long,  120  cm  wide,  and  122  cm  deep.  This  basin  was  supported  by  a  frame  constructed 
of  5.1  cm  aluminum,  the  base  of  which  was  adjustable  and  was  used  to  level  the  basin. 
Water  flowed  into  the  top  of  the  basin  and  out  of  the  bottom  side  of  the  basin  through  a 
15.2  cm  OD  PVC  outflow  pipe.  This  pipe  connected  to  both  a  7.6  cm  OD  PVC  re- 
circulating pipe  and  the  pump.  A  bypass  pipe  was  added  to  control  the  pump's  discharge. 

Current  profiles  were  measured  using  a  Kent  miniflo  model  265  propeller-type 
current  meter.  This  current  meter  was  attached  to  a  vernier  scale,  so  that  the  depth  of  the 
sensor  could  be  accurately  measured.  The  vernier  scale  was  attached  to  a  sliding 
aluminum  mount,  so  that  the  current  meter  could  be  moved  along  the  length  of  the  flume, 
and  be  maintained  in  a  position  perpendicular  to  the  flow  so  only  the  horizontal 
component  of  velocity  was  measured.  The  current  meter  and  the  vernier  scale  may  be 
seen  in  Figure  4.2. 

Seven  0.3  cm  holes  were  drilled  into  the  side  of  the  flume  at  positions  115  cm, 
165  cm,  and  425  cm  from  the  upstream  end  of  the  flume,  so  that  water  samples  at 
different  depths  could  be  extracted  during  experimental  runs  to  obtain  sediment 
concentrations  profiles,  and  a  manometer  was  used  to  determine  the  slope  of  the  water 
surface  across  three  different  points  along  the  length  of  the  flume.  The  manometer  and 
two  of  the  three  sets  of  side-holes  may  be  seen  in  Figure  4.3. 

A  false  bottom  was  constructed  of  3-cm-thick  Plexiglas  and  secured  to  the  bottom 
of  the  flume  so  that  a  sediment  bed  could  be  placed  into  the  flume.  The  false  bottom, 
seen  in  Figures  4.3  and  4.4,  was  made  into  three  segments.  The  front  segment  was  198.5 
cm  long,  2.5  cm  high,  the  width  of  the  inside  of  the  flume,  and  secured  at  the  upstream 
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end  of  the  flume  with  silicon  gel  so  that  it  would  neither  move  nor  leak  but  could  be 
removed  later  if  needed.  The  end  segment  was  20  cm  long,  2.5  cm  high,  the  width  of  the 
inside  of  the  flume,  and  similarly  secured  at  a  distance  of  400  cm  from  the  upstream  end 
the  flume.  The  middle  segment  was  226  cm  long,  2.5  cm  high,  the  width  of  the  inside  of 
the  flume,  and  designed  to  be  removed  and  replaced  with  a  sediment  bed.  It  could  lay 
flush  with  and  in  between  the  front  and  back  segments.  In  this  way,  when  emplaced  it 
resulted  in  a  flat  bottom  2.5  cm  above  the  flume  bottom.  As  seen  in  Figure  4.4,  the 
middle  segment  began  at  a  distance  of  199  cm  from  the  upstream  end  of  the  flume  and 
ended  at  a  distance  of  425  cm. 


Figure  4.2  Kent  propeller  type  current  meter  and  its  mount. 
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Figure  4.3  Manometer  and  sampling  locations. 
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Figure  4.4  Schematic  drawing  of  flume  detailing  sampling  points  and  false  bottom. 


4.3  Simulator  of  Erosion  Rate  Function 

The  Simulator  of  Erosion  Rate  Function  (SERF)  was  similar  in  principle  to  those 
designed  by  earlier  investigators  (Moore  and  Masch,  1962;  Arulanandan  et  al.,  1975; 
Croad,  1981;  Chapuis  and  Gatien,  1986).  Lee  and  Mehta  (1994)  have  provided  a  brief 
summary  of  the  underlying  design-and-operation  concept.  Typically,  these  devices 
consist  of  a  rotating  outer  cylinder,  which  is  filled  with  the  chosen  eroding  fluid.  A 
remolded  or  intact  sample  in  the  shape  of  a  smaller  cylinder  is  placed  into  the  outer 
cylinder.  While  rotating,  the  outer  cylinder  imparts  rotation  to  the  fluid,  which  translates 
into  a  shear  stress  being  uniformly  imparted  to  the  sample.  The  torque  applied  to  the 
stationary  sample  is  measured  and  converted  to  the  corresponding  shear  stress.  Erosion  of 
the  sample  follows,  and  a  relation  is  obtained  between  erosion  rate  and  shear  stress.  In 
SERF  a  load  cell  allowed  the  eroded  mass  to  be  measured  without  removing  the  sample 
assembly  from  the  outer  cylinder. 

Figure  4.5  shows  the  schematic  of  the  SERF,  and  Figure  4.6  is  a  photographic 
view.  The  clay  sample  is  molded  into  a  cylinder  with  a  diameter  of  7.6  cm  and  a  height  of 
9.6  cm.  The  cylinder  is  impaled  on  a  shaft  which  secures  an  upper  and  a  lower  disk  of 
equal  diameters  to  the  two  ends  of  the  sample.  The  upper  end  of  the  shaft  affixes  to  an 
Omega  (Stamford,  CT)  model  TQ201-100Z  torque  cell,  which  is  in  turn  connected  to  an 
Omega  model  LC601-5  load  cell.  When  the  acrylic  outer  cylinder  is  rotated,  the  chosen 
eroding  fluid  imparts  a  torque  to  the  sample,  and  in  turn  the  shaft  and  torque  cell.  The  cell 
is  connected  to  an  Omega  model  DP25-S  LED  (light-emitting  diode)  display,  and  shows 
the  applied  torque  in  g-cm.  The  load  cell  is  attached  to  the  entire  shaft/torque  cell/sample 
assembly,  and  can  be  used  to  measure  mass  eroded  after  a  given  run.  This  unit  can  move 
vertically  out  of  the  outer  cylinder  for  sample  and  water  replacements. 
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The  outer  cylinder  is  capable  of  rotating  at  a  maximum  of  2,350  rpm,  and  a 
minimum  of  150  rpm.  The  torque  cell  has  a  1,800  g-cm  capability,  and  the  load  cell  can 
accommodate  a  maximum  of  2.3  kg. 
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Figure  4.5  Schematic  of  SERF. 


55 


Figure  4.6  Front- view  of  SERF  showing  the  rotating  cylinder  assembly,  speed  control 
panel  (lower  left)  and  digital  readouts  (lower  right)  for  the  torque  cell  and  the  load  cell. 


4.4  Schulze  Ring  Shear  Tester 

The  Schulze  Ring  Shear  Tester  shown  in  Figure  4.7  is  similar  in  principle  to  the 
common  Jenike  shear  tester.  It  consists  primarily  of  a  lower  rotating  annular  tray  that 
holds  the  sample,  and  an  annular  disk  that  is  placed  on  top  the  sample  and  just  barely 
inside  the  bottom  tray  as  can  be  seen  in  Figure  4.8. 

The  lower  rotating  tray  is  used  to  apply  a  shear  stress  to  the  sample,  while  the 
upper  annular  disk  is  used  to  apply  a  normal  force.  The  lower  tray  is  connected  to  an 
electric  motor,  while  the  upper  tray  is  connected  to  two  horizontal  load  cells  in  order  to 
measure  the  torque  being  applied  to  the  sample,  and  is  also  connected  to  a  single  vertical 
load  cell  used  to  measure  the  normal  force.  These  components  may  be  seen  in  Figure 
4.9. 
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Figure  4.7  Schulze  Ring  Shear  Tester. 
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Figure  4.8  Cross-sectional  view  of  sample  inside  upper  annular  disk  and  lower  annular 
tray. 
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Figure  4.10  shows  a  schematic  of  the  upper  disk,  and  Figure  4.11  contains  a 
schematic  of  the  fins  that  were  attached  to  the  bottom  of  the  upper  disk  to  prevent  sample 
material  from  sliding. 
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Figure  4. 10  Schematic  of  upper  annular  disk. 


■  '         '  ■ 


1 


—  mst  »» 

1 

("•              <"» 

\ 

f                V 

1 

D3T 

-        \ 

Drill  clearance  for  10-32  bolt 


Figure  4. 1 1  Fin  schematic. 
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The  lower  tray  was  machined  out  of  solid  aluminum  while  the  upper  disk  was 
constructed  out  of  0.16  cm  thick  perforated  aluminum  supported  by  1.3  cm  thick 
Plexiglas  designed  according  to  Figure  4.10.  The  perforated  aluminum  on  the  upper  disk 
allowed  water  to  seep  through  when  a  normal  force  was  applied  to  the  sample.  This 
controlled  the  pore  water  pressure  and  kept  the  load  on  the  sediment  matrix. 
Additionally,  the  fins  were  screwed  onto  the  perforated  aluminum  of  the  upper  disk. 

The  Schulze  tester  allowed  for  automated  load  application.  The  vertical  load,  N, 
as  shown  in  Figure  4.9,  acting  on  the  bulk  solid  sample  could  be  adjusted  by  the  operator. 
The  tester  also  had  a  digitally  controlled  motor  drive  for  shear  cell  rotation.  Thus,  cell 
rotation  w  (Figure  4.9)  was  computer-controlled.  The  shear  force,  F\  =  Fi  (Figure  4.9), 
was  automatically  measured  and  the  output  was  viewed  using  a  Linseis  plotter.  Testing 
at  very  low  pressures  was  accomplished  with  a  counterweight  system  to  reduce  the 
weight  of  the  lid  and  other  parts  connected  to  the  lid.  Thus,  tests  at  stresses  <  500  Pa  were 
possible. 

4.5  Materials 

The  soil  materials  used  for  these  experiments  were  clay  minerals,  commercially 
prepared  clays,  and  sand  particles.  All  three  were  obtained  commercially,  and  were  well 
graded.  All  three  materials  were  analyzed  for  particle  size  distribution,  and  the  clay  were 
additionally  analyzed  for  pertinent  physicochemical  properties  such  as  the  Cation 
Exchange  Capacity  (CEC). 

Three  types  of  commercially  available  clays  were  selected  which  together  cover  a 
wide  range  of  cohesive  soil  properties;  a  kaolinite,  a  bentonite,  and  an  attapulgite. 
Kaolinite  (pulverized  kaolin),  a  light  beige-colored  powder,  was  obtained  from  the  EPK 
Division  of  Feldspar  Corporation  in  Edgar,  Florida.  Its  Cation  Exchange  Capacity  (CEC) 
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as  given  by  the  supplier,  was  5.2  -  6.5  milliequivalents  per  100  grams.  Its  granular 
density  was  2,603  kg/m3.  Bentonite  was  obtained  from  the  American  Colloid  Company 
in  Arlington  Heights,  Illinois.  It  was  a  sodium  montmorillonite  (commercial  name 
Volclay)  of  a  light  gray  color.  Its  CEC  was  105  milliequivalents  per  100  grams,  and  its 
granular  density  was  2,760  kg/m3.  Attapulgite,  of  a  greenish-white  color,  was  obtained 
from  Floridin  Company  in  Quincy,  Florida.  Also  known  as  palygorskite,  its  CEC  was  28 
milliequivalents  per  100  grams,  and  its  granular  density  was  2,300  kg/m3. 

The  median  particle  sizes  of  attapulgite,  bentonite,  and  kaolinite  were  0.86  |im, 
1.01  urn,  and  1.10  |im,  respectively.  These  three  clays  were  used  to  create  two  different 
mixtures.  Mixture  #1  was  composed  of  50%  kaolinite,  45%  attapulgite,  and  5% 
bentonite.  Mixture  #2  was  composed  of  50%  kaolinite,  35%  attapulgite,  and  15% 
bentonite. 

Three  commercially  prepared  ceramic  clays  were  obtained  from  Bennett  Pottery 
Supply  (Orlando,  Florida)  #10,  #20,  and  #75.  Clay  #10  had  a  d50  value  of  0.45  [im,  #20 
had  a  d50  of  0.36  ^m  and  #75  had  a  d5o  of  0.46  p.m.  The  cumulative  size  distribution 
plots  may  be  seen  in  Figures  4.12  through  4.14.  These  three  clays  were  also  analyzed 
using  a  computer-controlled  x-ray  diffractometer.  Each  of  the  clays  was  scanned  from  2° 
to  60°  of  20  using  copper  ka  radiation.  The  results  of  this  analysis  can  be  seen  in  Table 
4.1. 

Five  angular  quartz  sands  from  the  Feldspar  Corporation  in  Edgar,  Florida  were 
used  to  provide  a  range  of  d50  values.  The  sand,  which  had  the  largest  d50  value  of  1.40 
mm,  was  Feldspar  F  8-20.  The  cumulative  size  distribution  plot  for  this  sand  may  be 
seen  in  Figure  4.15  and  a  picture  of  the  grains  in  Figure  4.20.   The  second,  with  a  d50 
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value  of  0.94  mm,  was  Feldspar  sand  sieved  in  through  a  #20  screen  (Figures  4.16).  The 
third,  with  a  d5o  value  of  0.83  mm,  was  the  Feldspar  20-30  (Figures  4.17  and  4.21).  The 
fourth,  with  a  d50  value  of  0.47  mm,  was  Feldspar  sand  sieved  through  a  #30  screen 
(Figure  4.18).  The  smallest,  with  a  d50  value  of  0.41  mm,  was  Feldspar  sand  sieved 
through  a  #40  screen  (Figures  4.19  and  4.22). 


Table  4. 1  Composition  of  commercially  pre 

sared  clays  in  order  of  predominance 

Material 

Bennett  Pottery  #10 

Bennett  Pottery  #20 

Bennett  Pottery  #75 

Kaolinite 

Yes 

Yes 

Yes 

Quartz 

Yes 

Yes 

Yes 

Mica 

Yes 

Yes 

No 

Vermiculite 

Yes 

Yes 

No 
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Figure  4.12  Cumulative  size  distribution  plot  of  Bennett  Pottery  Supply  Clay  #10. 
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Percent  Finer  for  Bennett  Pottery  Supply  Clay  #20 
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Figure  4.  13  Cumulative  size  distribution  plot  of  Bennett  Pottery  Supply  Clay  #20. 
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Figure  4.14  Cumulative  size  distribution  plot  of  Bennett  Pottery  Supply  Clay  #75. 
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Percent  Finer  for  Feldspar  8-20  Sand 
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Figure  4.15  Cumulative  size  distribution  plot  of  Feldspar  F  8-20  sand. 
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Figure  4.16  Cumulative  size  distribution  plot  of  Feldspar  sand  sieved  through  a  #20 
screen. 
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Figure  4.17  Cumulative  size  distribution  plot  of  Feldspar  20-30  sand. 
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Figure  4.18  Cumulative  size  distribution  plot  of  Feldspar  sand  sieved  through  a  #30 
screen. 
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Figure  4.19  Cumulative  size  distribution  plot  of  Feldspar  sand  sieved  through  a  #40 
screen. 


Figure  4.20  Picture  of  sand  grains  of  Feldspar  sand  sieved  through  a  #40  screen. 
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Figure  4.21  Picture  of  sand  grains  of  Feldspar  20-30  sand. 


Figure  4.22  Picture  of  sand  grains  of  Feldspar  F  8-20  sand. 
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4.6  Flume  Calibration  and  Experimental  Procedure 

The  flume  was  calibrated  to  obtain  the  relationship  between  the  sluice  gate  height 
and  discharge,  so  that  the  discharge  could  be  calculated  knowing  the  water  surface  height 
immediately  before  the  gate  and  the  gate  height,  and  to  obtain  the  Manning's  coefficient 
of  the  Plexiglas,  so  that  the  wall  friction  could  later  be  subtracted  out  to  obtain  the  bottom 
shear  values  for  a  given  flow  condition. 

The  discharge  Q  was  related  to  the  gate  height  and  the  water  surface 
elevation  immediately  upstream  of  the  gate  through 

Q  =  ygbCdJlgAH  (4.6.1) 

where  yg  is  the  gate  height  in  meters,  b  is  the  flume  width  in  meters  (=0.15  m),  Q  is  the 
discharge  coefficient  (=  0.45),  g  is  acceleration  due  to  gravity,  and  AH  =  yi-yg  where  y;  is 
surface  elevation  immediately  upstream  of  the  gate.  Figures  4.23  and  4.24  show  the 
results  for  two  gate  heights. 

Manning's  n  value  was  calculated  by  using  the  standard  step  method  to  calculate 
the  back-water  profile  in  the  flume  for  a  series  of  runs  each  with  a  different  discharge. 
For  each  run,  the  surface  elevation  in  the  flume  was  measured  91  cm  upstream  of  the  bed, 
in  15.2  cm  increments  along  the  length  of  the  bed  starting  at  the  front  of  the  bed  and 
ending  at  the  back,  and  61  cm  from  the  back  of  the  bed.  In  the  used  for  each  run, 
Manning's  n  was  varied  until  the  calculated  and  measured  surface  elevations  matched. 
Figures  4.25  and  4.26  provide  an  example  of  two  different  runs. 
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Figure  4.23  Relationship  for  discharge  calculation  for  gate  height  of  3.18  cm. 
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Figure  4.24  Relationship  for  discharge  calculation  for  gate  height  of  3.81  cm. 
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Figure  4.25  Comparison  of  calculated  and  measured  surface  elevations  for  Run  20-1. 


0.1280 
?  0.1270 
I,  0.1260 
I   0.1250 


0.1240 


,5   0.1230 


c/> 


0.1220 


-2.50 


68 


Standard  Step  Method  for  Flume  Run  20-2 
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Figure  4.26  Comparison  of  calculated  and  measured  surface  elevations  for  Run  20-2. 

Manning's  n  value  for  each  run  was  then  averaged  to  obtain  a  value  of  0.014  for 
the  bed  which  compares  well  with  the  reported  value  for  Plexiglas  of  0.010. 


4.7  SERF  Calibration  and  Experimental  Procedure 

The  torque  cell  was  calibrated  using  an  aluminum  ring  affixed  to  the  shaft  of  the 
cell.  A  two-point  calibration  was  done  using  no  load,  and  then  with  a  1  kg  weight 
hanging  from  the  ring,  which  had  a  radius  of  1  cm.  This  resulted  in  a  calibration  between 
0  and  1  kg-cm.  The  load  cell  was  calibrated  in  the  same  two-point  fashion,  using  a  1  kg 
weight  hung  from  the  cell. 

The  calibration  to  convert  torque  into  shear  stress  required  running  the  device 
with  dummy  samples.  SERF  was  run  with  three  aluminum  cylinders  of  varying  diameters 
but  identical  height.  The  varying  of  diameters  provided  a  relationship  between  the 
applied  torque  and  the  sample  radius,  considered  to  be  necessary  in  instances  when  the 
radius  of  a  sample  decreases  measurably  due  to  erosion.  For  a  given  cylinder,  SERF  was 
rotated  from  a  standstill  to  a  "reasonable"  speed  (up  to  1,600  rpm)  in  time-steps,  and 
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torque  readings  were  taken  at  each  rpm  value  after  steady-state  was  reached  (typically 
within  5  seconds).  Once  at  peak  rpm,  the  device  was  slowed  in  time-steps  and  torque 
readings  were  taken  again.  This  procedure  was  done  three  times  for  each  cylinder,  the 
shaft  and  cylinder  being  removed  and  re-installed  each  time  in  order  to  gauge  the 
precision  of  the  device  after  re-installment  of  samples.  The  bottom  disk  of  the  sample 
holder  was  also  tested  in  this  manner.  By  subtracting  the  bottom  disk  torque  from  the 
total  torque,  the  torque  applied  to  the  sample  area  can  be  obtained.  The  water  meniscus 
was  located  just  below  the  upper  disk  (Fig.  4.5)  for  the  fluid  to  act  on  the  cylinder  and  the 
bottom  disk.  It  was  at  times  necessary  to  add  or  remove  water  to  obtain  the  appropriate 
meniscus. 

A  trend-line  was  plotted  for  the  variation  of  torque  with  rpm  for  the  disk,  and  the 
corresponding  relation  (Eq.  4.7.1)  was  used  to  determine  the  torque  provided  by  the 
sample  area  by  subtracting  the  bottom  disk  torque  (Fig.  4.5)  from  the  overall  torque.  This 
adjusted  torque  value  was  converted  into  shear  stress  using  Eqs.  4.7.1  and  4.7.3: 

T  =  0.0086RPM-1  (4.7.1) 

FT=-  (4.7.2) 

t  =  ^-  (4.7.3) 

\ 

where  T  is  the  torque  reading  in  kg-m  (g-cm  for  Eq.  4.7.1),  rc  is  the  cylinder  radius  in  m, 
Flan  is  the  applied  tangential  force  in  kg,  r  is  the  shear  stress  in  Pa,  g  is  the  acceleration 
due  to  gravity  in  m/s2,  and  Ac  is  the  cylinder  surface  area  in  m2.  The  shear  stress  values 
were  then  plotted  versus  rpm  to  yield  calibration  curves  for  three  cylinders  (Figs.  4.27, 
4.28,  4.29).  The  spread  in  data  was  attributed  to  slight  but  inevitable  changes  in  the 
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alignment  of  the  cylinders  with  respect  to  the  central  axis,  due  to  the  removal  and  re- 
installment  of  the  shaft/cylinder  assembly  between  runs.  This  procedure  of  removal/re- 
installment  was  carried  out  in  order  to  identify  the  degree  of  precision  of  the  device. 
Overall,  the  small  and  large  cylinders  followed  a  power  fit,  and  the  medium  cylinder 
followed  a  linear  trend  Figs.  4.27,  4.28,  4.29).  Fitting  a  power  function  to  the  medium 
cylinder  resulted  in  overestimation  of  shear  stress  at  high  rpm  (>  1,000). 


1 


Figure  4.27  Calibration  curve  for  large  aluminum  cylinder  (dia.=7.6  cm). 


Figure  4.28  Calibration  curve  for  medium  aluminum  cylinder  (dia.=7.2  cm). 
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Figure  4.29  Calibration  curve  for  small  aluminum  cylinder  (dia.=6.8  cm). 

For  the  large,  medium,  and  small  cylinder,  respectively,  the  following  relations 
were  obtained: 

T  =  4.9xlO^RPM196  (4.7.4) 

T  =  4.7  xl(T3RPM- 1.32  (4.7.5) 

T  =  8.0;cl(r8RPM2-63  (4.7.6) 

As  an  illustration  of  the  differences  among  the  above  relations  note  that  a  10% 

decrease  in  surface  area  from  the  small  cylinder  to  the  large  cylinder  led  to  a  67% 

increase  in  shear  stress  at  an  rpm  of  1,000  for  the  small  cylinder.  This  resulted  from  the 

decreased  moment  arm  of  the  small  cylinder  leading  to  an  increased  tangential  force  and 

therefore  increased  shear  stress. 

The  three  ceramic  clays  from  Bennett  Pottery  Supply  (Orlando,  FL),  #10,  #20, 
and  #75  were  used.  Each  clay  was  molded  to  two  different  bulk  densities  (both  lower 
than  initial)  by  adding  tap  water  and  working  by  hand.  This  procedure  was  repeated 
twice,  resulting  in  four  densities  tested  for  each  clay.    A  third  set  was  prepared  by 
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molding  two  equal  amounts  of  clay  with  water,  one  with  tap  water  and  another  with  35 
ppt  salt  (sodium  chloride)  water.  Water  was  added  in  equal  amounts  to  yield  a  similar 
density.  Since  the  clays  arrived  with  existing  pore  water,  the  actual  salinity  of  the  pore 
water  was  less  than  35  ppt. 

The  first  set  of  clays  was  shaped  into  cylinders  of  the  same  diameter  as  the 
medium  cylinder,  hence  the  medium  cylinder  calibration  curve  was  used.  The  second  set 
of  samples  was  shaped  into  the  diameter  of  the  large  cylinder,  as  was  the  third  set.  The 
volume  and  surface  area  of  each  sample  were  known,  and  sample  mass  was  determined 
by  weighing  the  sample  with  the  bracket  and  lid  (Figs.  4.5,  4.6),  and  then  weighing  only 
the  bracket  and  lid.  This  provided  the  density  of  each  respective  sample. 

Once  a  sample  was  emplaced  within  into  the  outer  cylinder  the  load  cell  was 
zeroed.  After  5  to  15  min  of  rotation  at  a  given  rpm,  the  device  was  stopped  and  allowed 
to  rest  for  approximately  one  minute.  The  actual  time-step  (between  5  and  15  min) 
depended  on  sample  erodibility.  In  most  cases  the  maximum  speed  of  the  outer  cylinder 
was  1,500  rpm.  Above  this  speed,  the  possibility  of  mass  erosion  (i.e.,  large  clasts  being 
eroded)  was  present,  which  could  have  resulted  in  an  imbalance  and  excessive  torque, 
thereby  damaging  the  torque  cell.  After  reaching  steady-state,  the  device  was  stopped  and 
the  load  cell  displayed  the  loss  in  material  in  grams. 

Since  the  load  cell  measured  buoyant  mass,  only  the  mass  loss  of  solids  was  being 
displayed.  This  value,  along  with  the  time-step  and  the  surface  area  resulted  in  the 
erosion  rate: 

Am, 
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where  Ams  is  the  dry  sediment  mass  loss,  Acs  is  the  sample  surface  area  in  m2,  and  At  is 
the  time  interval  in  seconds. 

A  comment  is  in  order  concerning  the  operating  speed  range  of  SERF,  the 
generation  of  Taylor  vortices  in  the  annular  gap  between  the  outer  cylinder  and  the 
sample,  and  the  development  of  turbulence  in  this  gap.  When  a  certain  Reynolds  number 
has  been  exceeded  in  the  gap  between  two  rotating  concentric  cylinders,  vortices  appear. 
The  axes  of  these  vortices  are  along  the  circumferences  of  the  cylinders,  and  they  rotate 
in  alternately  opposite  directions  (Schlichting,  1979).  The  Taylor  number  Ta  is  used  to 
predict  the  formation  of  these  irregularities 


U  w     \w 
r.  =  -^-i    -*-  (4.7.8) 

v     \Rt 

where  Ua  is  the  annular  peripheral  velocity  of  the  cylinder,  wg  is  the  gap  width,  /?,  is  the 

inner  cylinder's  radius,  and  v  is  the  kinematic  viscosity  of  the  fluid.  Taylor  vortices  are 

formed  in  the  range  41.3<  Ta  <400,  within  a  laminar  flow  regime  (Schlichting,  1968). 

Even  at  the  lowest  operating  speed  of  SERF,  which  corresponds  to  a  peripheral  velocity 

of  0.8  m/s  (at  150  rpm),  using  a  gap  width  of  0.024  m,  an  inner  cylinder  radius  of  0.038 

m,  and  water  as  the  fluid  ((v*=10'6  m2/s),  the  Taylor  number  was  15,000,  well  over  the 

range  of  Taylor  vortex  formation. 

With  the  characteristic  Reynolds  number  defined  as 

Uaw. 
Re  =  -i-i  (4.7.9) 

v 

using  an  annular  peripheral  velocity  of  0.8  m/s  and  the  same  fluid  and  gap  width,  Re  = 

18,880  is  obtained.  The  SERF  was  typically  operated  above  this  speed,  and  in  these  tests 

the  lowest  Re  attained  was  100,000  at  a  speed  of  4.3  m/s  (at  800  rpm).  Operating  above 
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the  Taylor  vortex  regime  and  above  the  laminar  flow  regime,  one  can  expect  a 
comparatively  uniform  distribution  of  shear  stress  over  the  sample  surface  area.  A 
detailed  operator's  manual  may  be  found  in  Ganju  et  al.  (2000). 

4.8  Schulze  Ring  Shear  Tester  Calibration  and  Experimental  Procedure 

The  calibration  of  the  Schulze  Ring  Shear  Tester  was  accomplished  at  the  start  of 
each  run  to  assure  that  the  data  being  collected  were  accurate.  The  calibration  was 
carried  out  using  an  80g  aluminum  weight  holder  designed  to  connect  to  the  load  cells 
and  to  carry  a  series  of  weighted  plates.  A  measure  was  first  taken  on  the  Linseis  strip 
chart  with  no  load.  Next,  the  80g  weight  holder  was  connected  and  another  reading  was 
taken  using  the  strip  chart.  After  this,  a  1  kg  weight  was  placed  onto  the  weight  holder 
and  a  reading  was  taken  using  the  strip  chart.  Again,  an  additional  1  kg  weight  was 
loaded  onto  the  weight  holder  and  a  reading  taken  completing  the  calibration  procedure. 
These  readings  provided  a  linear  scale  between  the  load  placed  on  the  load  cell  and  the 
displacement  on  the  strip  chart  against  which  readings  taken  during  tests  could  be  scaled. 

Next,  the  sample  bulk  solid  was  placed  into  the  shear  cell  and  preconsolidated 
manually.  It  was  then  sheared  at  a  normal  stress  <7pre  which  was  adjusted  by  a  specific 
normal  force  N.  The  shear  force  S,  acting  over  area  A,  increased  with  time  t  as  indicated 
in  the  left  diagram  of  Figure  4.30.  Besides  the  shear  stress,  the  bulk  density  of  the  sample 
also  increased  as  the  sand  matrix  decreased  in  volume  as  sand  grains  moved  to  a  tighter, 
stronger  packing  structure.  After  a  certain  time  (on  the  order  of  10  minutes)  a  constant 
bulk  density  and  a  constant  shear  stress  (rpre)  prevailed. 
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preshear       shear 


Figure  4.30  Plot  of  the  shear  stresses  with  time  illustrating  the  data  collection  process, 
and  the  corresponding  yield  locus. 


Both  bulk  density  and  shear  stress  are  characteristic  for  the  normal  stress  crpre, 
which  means  that  each  normal  stress  <7pre  results  in  a  specific  constant  bulk  density  and 
shear  stress.  The  deformation  of  the  bulk  solid  at  constant  normal  and  shear  stresses  and 
constant  bulk  density  represents  a  steady-state  (or  stationary)  flow,  and  the  process  of 
attaining  steady-state  is  called  "preshear".  With  the  preshear  process  the  bulk  solid  is 
placed  into  a  defined  consolidation  state  (similar  to  the  application  of  the  consolidation 
stress  in  a  uniaxial  compression  test).  Preshear  is  arrested  once  steady  state  is  reached  by 
reversing  the  direction  of  rotation  of  the  sample  disk  until  S  =  0  (t  =  0).  The  values  of 
normal  stress  <7pre  and  the  corresponding  shear  stress  Tpre  form  the  point  of  steady  state 
flow  in  the  cr,T-diagram  (Figure  4.30,  right  diagram) 

Afterwards,  the  sample  is  sheared  again  with  a  reduced  normal  stress  cSh  <  <7pK. 
Since  the  sample  is  now  sheared  under  a  smaller  normal  load  than  at  preshear,  it  will  start 
to  flow  (break)  at  a  certain  shear  stress  (similar  to  the  breakage  of  the  sample  in  a 
uniaxial  test). 
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A  bulk  solid  starts  to  flow  when  the  Mohr  stress  circle  representing  the  actual 
stress  state  touches  the  yield  locus.  The  start  of  flow  is  connected  to  a  decrease  in  the 
bulk  density  and  a  corresponding  reduction  of  the  shear  stress  (see  Figure  4.30).  The 
maximum  in  the  plot  of  the  shear  stress  versus  time  indicates  the  start  of  flow.  The 
corresponding  normal  and  the  shear  stresses  describe  a  point  on  the  yield  locus  in  the  a,  t- 
diagram.  If  several  samples  of  bulk  solid  are  presheared  under  the  identical  normal  stress 
(Tpre  but  sheared  under  different  normal  stresses  <7Sh  <  opK,  the  path  of  the  yield  locus  in 
the  a,  T-diagram  can  be  determined. 

The  parameters  which  describe  the  flow  properties  can  be  determined  from  the 
yield  locus  (Figure  4.31).  The  relevant  consolidation  stress  0\  is  equal  to  the  major 
principal  stress  of  the  Mohr  stress  circle  which  is  tangential  to  the  yield  locus  and 
intersects  at  the  point  of  steady  state  flow  (crpre  ,  Tpre).  This  stress  circle  represents  the 
stresses  in  the  sample  at  the  end  of  the  consolidation  process  (stresses  at  steady  state 
flow).  The  unconfmed  yield  strength  ac  results  from  the  stress  circle  which  is  tangential 
to  the  yield  locus  and  which  runs  through  the  origin  (minor  principal  stress  o2=  0).  This 
stress  circle  represents  the  same  stress  state  as  the  one  which  prevails  in  the  second  step 
of  a  uniaxial  compression  test.  In  contrast  to  a  uniaxial  compression  test,  the  unconfmed 
yield  strength  ac  has  to  be  determined  on  basis  of  the  yield  locus  and  does  not  follow 
directly  from  the  measurement. 

The  tangent  at  the  largest  Mohr  stress  circle  is  called  the  effective  yield  locus.  It 
encloses  the  a-axis  with  the  angle  cpe  (effective  angle  of  internal  friction).  Because  the 
largest  Mohr  stress  circle  indicates  a  steady  state  flow,  (pe  can  be  used  as  a  measure  of  the 
internal  friction  at  steady  state  flow. 
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yield  locus 


effective  yield  locus 


Figure  4.31  Yield  locus. 


The  angle  of  internal  friction  cpi  is  defined  as  the  local  slope  angle  of  the  yield 
locus.  Because  the  yield  locus  is  curved  instead  of  being  straight,  (p,  varies  along  the  yield 
locus.  It  is  sufficient  for  most  applications  to  state  one  value  of  (p,.  In  this  case,  the 
internal  friction  angle  is  defined  by  the  slope  (piin  of  the  linearized  yield  locus.  The 
linearized  yield  locus  results  from  the  tangent  common  to  both  Mohr  stress  circles  shown 
in  Figure  4.31  (Schulze,  1994). 


CHAPTER  5 
EXPERIMENTAL  RESULTS 

5.1  SERF  Experiments 

Eighteen  clay  samples  were  tested  in  total,  and  for  all  but  three  samples 
(#2,  #9  and  #15)  erosion  was  detected  by  the  load  cell.  The  experimental  conditions 
under  which  these  samples  (less  2,  #9  and  #15)  were  tested  are  listed  in  Table  5.1. 


Table  5.1  SERF  experimental  conditions 

Sample 

# 

Clay# 
(Bennet  Pottery  Supply) 

Pore  Water 
Composition 

Density, 

P 
(kg/m3) 

Range  of  Applied 

Shear  Stress 

(Pa) 

1 

10 

fresh 

1665 

2.9  -  5.8 

3 

10 

fresh 

1710 

2.8-5.9 

4 

10 

fresh 

1859 

4.2-8.1 

5 

10 

saline 

1928 

2.3  -  6.5 

6 

10 

fresh 

1928 

2.2  -  7.5 

7 

20 

fresh 

1435 

1.6-3.9 

8 

20 

fresh 

1537 

4.1-5.4 

10 

20 

fresh 

1721 

2.8  -  4.5 

11 

20 

saline 

1894 

2.9  -  5.8 

12 

20 

fresh 

1905 

4.1-6.5 

13 

75 

fresh 

1675 

1.1-4.5 

14 

75 

fresh 

1698 

2.0  -  4.5 

16 

75 

fresh 

1806 

2.9  -  5.7 

17 

75 

saline 

1940 

3.1-6.4 

18 

75 

fresh 

1963 

3.2-6.8 
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A  dimensionless  erosion  rate  as  a  function  of  dimensionless  shear  stress  plot  was 
obtained  as  may  be  seen  in  Figure  5.1.  For  that  purpose,  linear  trend-lines  were  fitted  to 
each  set  of  results  with  the  following  form: 

«-*.(*» -O  (5.1.1) 

where  e  is  the  erosion  rate  in  kg/m2-s,  ts  is  the  sample  shear  strength  in  Pa,  eN  is  the 
erosion  rate  constant  in  kg/N-s  and  xb  is  the  shear  stress  in  Pa.  A  summary  of  these 
results  may  be  seen  in  Table  5.2.  These  trend  lines  were  then  made  dimensionless  as 
follows: 
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rable  5.2  Density,  erosion  rate  constant,  anc 

shear  strength  values 

?or  15  clay  samples 

Sample  # 

Density,  p 
(kg/m3) 

(kg/m2-s) 

(Pa) 

1 

1665 

5.00  x  10"5 

3.14 

3 

1710 

1.91  x  10"5 

3.38 

4 

1859 

1.23  xlO"5 

4.49 

5 

1928 

1.73  x  10s 

0.60 

6 

1928 

1.20  xlO"5 

1.67 

7 

1435 

3.93  x  10"5 

1.93 

8 

1537 

1.19xl0"5 

2.81 

9 

1721 

7.72  x  10~5 

3.02 

10 

1894 

7.00  x  10"5 

1.58 

11 

1905 

7.10xl0"6 

2.00 

12 

1675 

4.80  x  10"5 

1.39 

13 

1698 

3.22  x  105 

1.73 

14 

1806 

1.14xl0"5 

1.76 

16 

1940 

6.89  xlO"5 

3.27 

17 

1963 

9.85  x  10"6 

1.88 

80 
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Figure  5.1  Dimensionless  erosion  rate  as  a  function  of  dimensionless  shear  stress  for  the 
clay  samples  run  in  SERF. 


It  should  be  pointed  out  that  the  selected  group  of  samples  did  not  erode  enough 
to  significantly  decrease  the  sample  radius,  and  only  one  cylinder  calibration  curve  was 
used  for  each  sample.  For  example,  sample  #6,  with  a  bulk  density  (wet)  of  1928  kg/m3, 
lost  a  total  of  3  g  during  the  entire  test.  This  corresponds  to  a  reduction  in  thickness  of 
0.007  cm,  or  0.18%  thickness  lost  (based  on  large  cylinder  volume). 


5.2  Flume  Experiments 

After  the  flume  was  calibrated,  it  was  run  with  the  clays  from  Bennett  Pottery 
Supply,  the  clays  obtained  from  the  EPK  Division  of  Feldspar  Corporation,  and  with  the 
five  angular  quartz  sands  from  the  Feldspar  Corporation.    A  total  of  44  samples  were 
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tested,  and  the  experimental  conditions  under  which  these  samples  were  run  are  listed  in 
Table  5.3. 


Table  5.3  Summary  of  flume  experimental 

conditions 

Sample 

# 

Sand 
Type 

Clay  Type 

Fine 
Fraction 

Pore 

Water 

Bulk 
Density, 
(kg/m5) 

Applied 

Shear  Stress 

(Pa) 

1 

N/A 

Bennett  Pottery 
#20 

1.00 

fresh 

1232 

0.13-  0.79 

2 

N/A 

Bennett  Pottery 
#75 

1.00 

fresh 

1259 

0.16-0.52 

3 

N/A 

Bennett  Pottery 
#10 

1.00 

fresh 

1240 

0.09  -  0.84 

4 

N/A 

Bennett  Pottery 
#10 

1.00 

fresh 

1288 

0.14-0.88 

5 

#30 

N/A 

0.00 

fresh 

1850 

0.10-0.59 

6 

#2 

N/A 

0.00 

fresh 

1850 

0.10-0.63 

7 

#20 

N/A 

0.00 

fresh 

1850 

0.20  -  0.93 

8 

F8-20 

N/A 

0.00 

fresh 

1850 

0.12-0.98 

9 

20-30 

N/A 

0.00 

fresh 

1850 

0.12  -  0.75 

10 

F8-20 

Kaolinite 

0.02 

fresh 

1832 

0.12-1.40 

11 

F8-20 

Kaolinite 

0.04 

fresh 

1886 

0.11-0.99 

12 

F8-20 

Kaolinite 

0.06 

fresh 

1895 

0.12-0.99 

13 

F8-20 

Kaolinite 

0.08 

fresh 

1911 

0.12-0.59 

14 

F8-20 

Kaolinite 

0.10 

fresh 

1974 

0.12-0.80 

15 

20-30 

Kaolinite 

0.02 

fresh 

1770 

0.12-0.56 

16 

20-30 

Kaolinite 

0.04 

fresh 

1851 

0.12-0.49 

17 

20-30 

Kaolinite 

0.06 

fresh 

1912 

0.12-0.28 

18 

20-30 

Kaolinite 

0.08 

fresh 

2013 

0.12-0.41 

19 

20-30 

Kaolinite 

0.10 

fresh 

2022 

0.12-0.42 

20 

20-30 

Kaolinite 

0.12 

fresh 

2086 

0.12-0.78 

21 

20-30 

Kaolinite 

0.14 

fresh 

2128 

0.12-1.33 

22 

20-30 

Mixture  #1 

0.02 

saline 

1863 

0.07  -  0.49 

23 

20-30 

Mixture  #1 

0.04 

saline 

1951 

0.12-0.70 

24 

20-30 

Mixture  #1 

0.06 

saline 

2169 

0.12-0.86 

25 

20-30 

Mixture  #1 

0.08 

saline 

2416 

0.12-1.81 
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Table  5.3  (continued) 

Sample 

# 

Sand 
Type 

Clay  Type 

Fine 
Fraction 

Pore 

Water 

Bulk 
Density, 
(kg/m*) 

Applied 

Shear  Stress 

(Pa) 

26 

20-30 

Mixture  #2 

0.02 

saline 

1832 

0.12-0.66 

27 

20-30 

Mixture  #2 

0.04 

saline 

1903 

0.12-0.65 

28 

20-30 

Mixture  #2 

0.06 

saline 

2102 

0.12-0.74 

29 

20-30 

Mixture  #2 

0.08 

saline 

2394 

0.12-1.07 

30 

20-30 

Mixture  #2 

0.10 

saline 

2539 

0.12-1.53 

31 

20-30 

Kaolinite 

0.02 

saline 

1853 

0.12-0.63 

32 

20-30 

Kaolinite 

0.04 

saline 

2050 

0.12-0.53 

33 

20-30 

Kaolinite 

0.06 

saline 

2140 

0.12-0.48 

34 

20-30 

Kaolinite 

0.08 

saline 

2307 

0.12-0.57 

35 

20-30 

Kaolinite 

0.10 

saline 

2417 

0.12-0.68 

36 

20-30 

Kaolinite 

0.12 

saline 

2550 

0.12-0.86 

37 

20-30 

Kaolinite 

0.14 

saline 

2574 

0.12-0.96 

38 

#2 

Kaolinite 

0.02 

fresh 

1847 

0.12-0.56 

39 

#2 

Kaolinite 

0.04 

fresh 

2057 

0.12-0.20 

40 

#2 

Kaolinite 

0.06 

fresh 

2215 

0.12-0.56 

41 

#2 

Kaolinite 

0.08 

fresh 

2301 

0.12-0.50 

42 

#2 

Kaolinite 

0.10 

fresh 

2434 

0.12-0.57 

43 

#2 

Kaolinite 

0.12 

fresh 

2554 

0.12-0.72 

44 

#2 

Kaolinite 

0.14 

fresh 

2571 

0.12-1.33 

The  flume  was  first  run  with  the  three  different  Bennett  Pottery  clay  beds.  The 
results  from  these  runs  may  be  seen  in  Figures  5.2  -  5.4.  These  results  were  used  to 
check  flume  calibration  as  discussed  in  Section  5.4. 
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Figure  5.2  Erosion  rate  plot  for  Bennett  Pottery  Supply  Clay  #20. 
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Figure  5.3  Erosion  rate  plot  for  Bennett  Pottery  Supply  Clay  #75. 
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Figure  5.4  Erosion  rate  plot  for  Bennett  Pottery  Supply  Clay  #2 
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The  flume  was  then  run  with  the  five  different  sand  beds,  each  with  a  different 
median  diameter  (d50),  and  the  critical  shear  stress,  as  defined  in  Section  3.2,  for  each  bed 
was  determined.  The  fluid  velocity,  and  consequently  the  bottom  shear  stress  was 
increased  incrementally  until  sand  grain  motion  was  observed  to  be  restricted  to  isolated 
spots  with  a  small  number  of  grains  moving  briefly  a  few  centimeters.  After  each 
increment  of  velocity  change,  the  water  surface  height  was  recorded  using  the 
manometers  to  determine  both  the  surface  slope  across  the  bed  and  the  discharge  through 
the  sluice  gate. 

The  critical  shear  stress,  the  Reynolds  number, Re,,  and  the  critical  Shields 
parameter,  6cr,  for  each  bed  may  be  seen  in  Table  5.4.  The  critical  Shields  parameter  for 
the  cohesionless  beds  was  plotted  as  a  function  of  Re,,  and  compared  to  values  obtained 
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by  Mehta  and  Partheniades  (1979),  Panagiatopolous  et.  al.  (1997),  and  Torfs  (1995). 
Figure  5.5  shows  the  results  from  the  flume  to  compare  well  with  the  results  of  these 
other  authors. 


Table  5.4  Critical  shear  stress,  Reynolds  number,  and  critical  Shields  parameter  for  sand 
beds 


Sand  Type 

^50 

(m) 

Critical 

Shear  Stress 

(Pa) 

«. 

(m/s) 

Re* 

Bcr 

#2 

0.00041 

0.352 

0.0187 

7 

0.0385 

#30 

0.00057 

0.360 

0.0190 

10 

0.0390 

20-30 

0.00083 

0.498 

0.0223 

17 

0.0371 

#20 

0.00094 

0.609 

0.0247 

21 

0.0400 

F8-20 

0.0014 

0.765 

0.0277 

35 

0.0338 

The  flume  was  then  used  to  test  Samples  11  through  46  (Table  5.3).  Results  from 
these  tests  are  seen  in  Figures  5.6  to  5.13.  Figure  5.12  is  a  comparison  plot  of  three 
different  sediment  beds,  each  with  fresh  kaolinite  as  the  fine  component  and  a  different 
d5o  for  sand,  to  illustrate  the  effect  of  grain  size  on  the  behavior  of  the  critical  shear  stress 
of  the  sand/kaolinite  mixture.  From  this  figure  it  can  be  seen  that  they^  changes  for 
each  bed. 

From  Figure  5.12  it  is  seen  that  y^  =  0.045  for  d50  =  0.41  mm,  vj/^  =  0.055  for 
d50  =  0.83  mm,  and  y^  =  0.55  for  d5o  =  1.40  mm,  which  indicates  that  the  critical  shear 
stress  is  not  affected  significantly  by  sand  size. 

Figure  5.13  is  a  comparison  plot  of  four  different  sediment  beds,  each  with  the 
same  d5o  (=  0.83  mm)  and  a  different  clay,  to  illustrate  the  effect  of  clay  on  the  critical 
shear  stress  for  sand.  It  is  seen  that  \|/rain  =  0.06  for  the  fresh  kaolinite,  y^  =  0.053  for 
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the  saline  kaolinite,  y^^  0.04  for  the  clay  mixture  #1,  and  v|/min=  0.03  for  the  clay 

mixture  #2,  indicating  the  critical  shear  stress  is  affected  by  the  cohesive  component  of 
sediment. 

Finally,  the  minimum  critical  shear  value,  Tcs,  for  each  series  of  beds  was  used 
with  Equation  3.3.34  to  determine  the  critical  Shields  parameter,  6cr,  for  each  series  and 
these  values  were  then  compared  with  those  calculated  for  Figure  5.5.  The  values  for  the 
parameters  in  Equation  3.3.34  were  the  same  as  those  given  in  Section  3.3  with  the 
exception  of  £>  =  65°  and  <J>eff  =  35°.  This  comparison  may  be  seen  in  Figure  5.14. 
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Figure  5.5  Comparison  of  critical  Shields  parameter  of  sand  beds  with  the  results  of 
other  researchers  and  with  the  Shields  curve. 


87 


0.00 


Kaolinite  and  F  8-20  Mixture 


I.OU 

ll0mt 

S.    1-40 

la 

♦ 

S    1.20 

sz 

CO 

«    1-00  < 

"■p 

♦ 

♦ 

♦ 

♦ 

g   0.80 
0.60 

0.05  0.10 

Fine  Fraction 


0.15 


♦     Kaolinite  Mixture 


Sand 


Figure  5.6  Plot  of  the  behavior  of  the  critical  shear  stress  of  the  Feldspar  F  8-20  and 
kaolinite  with  fresh  pore  water  series  of  beds. 
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Figure  5.7  Plot  of  the  behavior  of  the  critical  shear  stress  of  the  Feldspar  F  20-30  and 
kaolinite  with  fresh  pore  water  series  of  beds. 
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Figure  5.8  Plot  of  the  behavior  of  the  critical  shear  stress  of  the  Feldspar  F  20-30  and  clay 
mixture  #1  with  3  ppt.  saline  pore  water  series  of  beds. 
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Figure  5.9  Plot  of  the  behavior  of  the  critical  shear  stress  of  the  Feldspar  F  20-30  and  clay 
mixture  #2  with  3  ppt.  saline  pore  water  series  of  beds. 
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Figure  5.10  Plot  of  the  behavior  of  the  critical  shear  stress  of  the  Feldspar  F  20-30  and 
kaolinite  with  3  ppt.  saline  pore  water  series  of  beds. 
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Figure  5. 1 1  Plot  of  the  behavior  of  the  critical  shear  stress  of  the  Sand  #2  and  kaolinite 
with  fresh  pore  water  series  of  beds. 
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Variation  of  Normalized  Critical  Shear  with  Sediment  Size 
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Figure  5.12  Variation  of  the  normalized  critical  shear  stress  with  fine  weight  fraction  for 
three  different  sands. 
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Figure  5.13  Variation  of  the  critical  shear  stress  with  fine  weight  fraction  for  20-30  sand 
and  different  clays. 
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Critical  Shields  Parameter 
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+ 
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0.83mm  sand  &  clay  1 

o 
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Figure  5.14  Comparison  of  the  critical  Shields  parameter  of  sand  beds  with  kaolinite 
calculated  using  Equation  3.3.34  compared  with  the  results  of  Mehta  and  Partheniades 
(1979)  and  with  the  critical  Shields  parameter  of  sand  beds  without  kaolinite  calculated 
using  Equation  3.3.17. 


From  Figure  5.14,  it  is  seen  that  the  introduction  of  clay  particles  in  pore  water  to 
the  sand  grain  matrix  results  in  a  significant  decrease  in  the  particle  packing  angle  for  the 
sand  (from  65o  to  35o),  hinting  at  the  mechanism  causing  the  reduction  in  critical  shear 
stress  for  sand.  The  particle  packing  angle,  analogous  to  the  angle  of  repose  which  is 
often  referred  to  as  the  friction  angle,  by  definition  incorporates  the  effects  of 
interparticle  friction  between  the  sand  grains.  Consequently,  the  reduction  in  this  angle 
(from  the  introduction  of  interstitial  clay  particles)  represents  a  lubrication  between  the 
sand  grains. 
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5.3  Schulze  Shear  Tester  Experiments 

Sixty  samples  were  tested  in  total,  and  for  all  but  two  samples  whose  shear  stress 
values  were  less  than  the  unconfined  yield  strength  stress  circle,  reasonable  values  were 
detected  by  the  shear  cell.  The  experimental  conditions  under  which  these  samples  were 
tested  are  listed  in  Table  5.5. 


Ta 

3le  5.5  Schulze  Shear  Tester  experimenta 

conditions 

Sample 

# 

Sand 
Size 
(mm) 

Clay  Type 

Fine 
Fraction 

Pore 

Water 

Preshear 

Load 

GO 

Applied 

Shear  Stress 

(kPa) 

1 

0.41 

Kaolinite 

0.00 

fresh 

6000 

0.00  -  6.23 

2 

0.41 

Kaolinite 

0.00 

fresh 

4000 

0.00  -  4.50 

3 

0.41 

Kaolinite 

0.00 

fresh 

2000 

0.00  -  2.49 

4 

0.41 

Kaolinite 

0.00 

fresh 

1000 

0.00-1.81 

5 

0.41 

Kaolinite 

0.02 

fresh 

6000 

0.00  -  7.03 

6 

0.41 

Kaolinite 

0.02 

fresh 

4000 

0.00  -  4.90 

7 

0.41 

Kaolinite 

0.02 

fresh 

2000 

0.00  -  2.37 

8 

0.41 

Kaolinite 

0.02 

fresh 

1000 

0.00-1.93 

9 

0.41 

Kaolinite 

0.04 

fresh 

6000 

0.00  -  7.35 

10 

0.41 

Kaolinite 

0.04 

fresh 

4000 

0.00  -  5.28 

11 

0.41 

Kaolinite 

0.04 

fresh 

2000 

0.00  -  2.44 

12 

0.41 

Kaolinite 

0.04 

fresh 

1000 

0.00-1.88 

13 

0.41 

Kaolinite 

0.06 

fresh 

6000 

0.00-5.51 

14 

0.41 

Kaolinite 

0.06 

fresh 

4000 

0.00-4.51 

15 

0.41 

Kaolinite 

0.06 

fresh 

2000 

0.00-2.18 

16 

0.41 

Kaolinite 

0.06 

fresh 

1500 

0.00-1.67 

17 

0.41 

Kaolinite 

0.08 

fresh 

6000 

0.00  -  6.75 

18 

0.41 

Kaolinite 

0.08 

fresh 

4000 

0.00  -  4.29 

19 

0.41 

Kaolinite 

0.08 

fresh 

2000 

0.00-3.18 

20 

0.41 

Kaolinite 

0.08 

fresh 

1500 

0.00-1.57 

21 

0.41 

Kaolinite 

0.10 

fresh 

6000 

0.00-5.31 

22 

0.41 

Kaolinite 

0.10 

fresh 

4000 

0.00  -  3.95 

23 

0.41 

Kaolinite 

0.10 

fresh 

2000 

0.00  -  2.03 
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Table  5.5  (continued) 

Sample 

# 

Sand 
Size 

(mm) 

Clay  Type 

Fine 
Fraction 

Pore 

Water 

Preshear 

Load 

(K) 

Applied 

Shear  Stress 

(kPa) 

24 

0.41 

Kaolinite 

0.10 

fresh 

1500 

0.00-1.56 

25 

0.41 

Kaolinite 

0.12 

fresh 

6000 

0.00-5.12 

26 

0.41 

Kaolinite 

0.12 

fresh 

4000 

0.00  -  3.56 

27 

0.41 

Kaolinite 

0.12 

fresh 

2000 

0.00-1.89 

28 

0.41 

Kaolinite 

0.12 

fresh 

1500 

0.00-1.49 

29 

0.41 

Kaolinite 

0.02 

saline 

6000 

0.00  -  5.25 

30 

0.41 

Kaolinite 

0.02 

saline 

4000 

0.00  -  3.74 

31 

0.41 

Kaolinite 

0.02 

saline 

2000 

0.00-2.12 

32 

0.41 

Kaolinite 

0.02 

saline 

1500 

0.00-1.69 

33 

0.41 

Kaolinite 

0.04 

saline 

6000 

0.00-5.21 

34 

0.41 

Kaolinite 

0.04 

saline 

4000 

0.00  -  3.80 

35 

0.41 

Kaolinite 

0.04 

saline 

2000 

0.00-1.964 

36 

0.41 

Kaolinite 

0.04 

saline 

1500 

0.00-1.58 

37 

0.41 

Kaolinite 

0.06 

saline 

6000 

0.00  -  5.21 

38 

0.41 

Kaolinite 

0.06 

saline 

4000 

0.00  -  3.68 

39 

0.41 

Kaolinite 

0.06 

saline 

2000 

0.00-1.89 

40 

0.41 

Kaolinite 

0.06 

saline 

1500 

0.00-1.48 

41 

0.41 

Kaolinite 

0.08 

saline 

6000 

0.00  -  5.05 

42 

0.41 

Kaolinite 

0.08 

saline 

4000 

0.00  -  3.59 

43 

0.41 

Kaolinite 

0.08 

saline 

2000 

0.00-1.81 

44 

0.41 

Kaolinite 

0.08 

saline 

1500 

0.00- 1.42 

45 

0.41 

Kaolinite 

0.10 

saline 

6000 

0.00  -  4.86 

46 

0.41 

Kaolinite 

0.10 

saline 

4000 

0.00  -  3.23 

47 

0.41 

Kaolinite 

0.10 

saline 

2000 

0.00-1.77 

48 

0.41 

Kaolinite 

0.10 

saline 

1500 

0.00-1.27 

49 

0.41 

Kaolinite 

0.12 

saline 

6000 

0.00  -  4.64 

50 

0.41 

Kaolinite 

0.12 

saline 

4000 

0.00  -  3.36 

51 

0.41 

Kaolinite 

0.12 

saline 

2000 

0.00-1.72 

52 

0.41 

Kaolinite 

0.12 

saline 

1500 

0.00-1.26 

53 

0.41 

Kaolinite 

0.01 

fresh 

6000 

0.00  -  4.98 

54 

0.41 

Kaolinite 

0.01 

fresh 

4000 

0.00  -  3.74 

55 

0.41 

Kaolinite 

0.01 

fresh 

2000 

0.00  -  2.58 
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Table  5.5  (continued) 


56 

0.41 

Kaolinite 

0.01 

fresh 

1500 

0.00-1.68 

57 

0.41 

Kaolinite 

0.01 

saline 

6000 

0.00-5.13 

58 

0.41 

Kaolinite 

0.01 

saline 

4000 

0.00  -  4.24 

59 

0.41 

Kaolinite 

0.01 

saline 

2000 

0.00-2.12 

60 

0.41 

Kaolinite 

0.01 

saline 

1500 

0.00-1.57 

The  Schulze  shear  tester  was  first  run  with  the  series  of  beds  composed  of  Sand 
#2  and  kaolinite  with  fresh  pore  water.  The  results  from  these  runs  are  summarized  in 
Table  5.6  and  illustrated  in  Figure  5.15.  The  behavior  of  the  friction  factor  for  this  series 
with  change  in  the  fine  fraction  may  be  seen  in  Figure  5.16.  The  shear  tester  was  next 
run  with  the  series  of  beds  composed  of  Sand  #2  and  kaolinite  with  saline  pore  water. 
The  results  from  these  runs  are  first  summarized  in  Table  5.7  and  then  illustrated  in 
Figure  5.17.  The  behavior  of  the  friction  factor  for  this  series  with  change  in  the  fine 
fraction  may  be  seen  in  Figure  5.18. 

The  friction  factors  were  determined  by  extrapolating  the  critical  shear  stress  back 
to  a  consolidation  stress  of  1000  Pa.  The  unconfined  yield  strength  data  for  a  given  fine 
fraction  was  first  plotted  as  a  function  of  consolidation  stress,  calculated  using  the  applied 
normal  load  and  the  area  of  the  annular  disk.  This  provided,  for  each  fine  fraction,  a 
series  of  data  points  that  could  be  used  for  a  linear  regression.  The  line  was  extrapolated 
backwards  to  determine  the  unconfined  yield  strength  occurring  at  a  consolidation  stress 
of  1000  Pa,  which  was  a  consolidation  stress  much  closer  to  that  of  a  naturally  occurring 
bed  than  what  the  ring  shear  tester  could  reproduce.  The  angle  of  internal  friction  for  a 
given  fine  fraction  occurring  at  1000  Pa  was  determined  using  the  same  technique.  Once 
these  two  values  were  known,  the  shear  stress  at  incipient  flow  was  calculated  from 
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Tc=-occos<p( 


(5.1) 


where  ac  is  the  unconfined  yield  stress  and  <p,  is  the  angle  of  internal  friction.  From  the 
shear  stress  the  friction  factor  could  be  calculated  from  its  definition: 


N     1000 


(5.2) 


where  N  is  the  load.  This  process  was  repeated  for  each  fine  fraction,  resulting  in  a 
series  of  data  points  for  beds  of  Sand  #2  with  fresh  kaolinite  and  Sand  #2  with  saline 
kaolinite. 


Table  5.6  Schulze  Shear  Tester  results  for  fresh  kaolinite  series 


Fine 
Fraction 

Consolidation 

Stress 

(g) 

Unconfined 

Yield 
Strength,  gc 

(kPa) 

Angle  of 

Internal 

Friction,  <p, 

Shear  Stress 
(kPa) 

0.00 

2497 

1.791 

38.35 

0.702 

0.00 

4504 

2.198 

38.09 

0.865 

0.00 

6234 

2.193 

37.94 

0.865 

0.01 

1911 

0.477 

44.52 

0.170 

0.01 

3527 

0.581 

43.97 

0.209 

0.01 

5167 

0.675 

43.89 

0.243 

0.02 

2132 

0.551 

44.59 

0.196 

0.02 

3673 

0.689 

42.59 

0.254 

0.02 

5376 

0.826 

41.53 

0.309 

0.04 

2007 

0.584 

44.50 

0.208 

0.04 

3423 

0.713 

42.60 

0.262 

0.04 

5238 

0.858 

41.59 

0.321 

0.06 

1682 

0.586 

45.11 

0.207 

0.06 

2332 

0.634 

43.26 

0.231 

0.06 

3582 

0.776 

42.51 

0.286 

0.06 

5641 

0.938 

41.04 

0.354 

0.08 

1569 

0.894 

39.04 

0.347 

0.08 

4288 

1.932 

36.78 

0.774 

0.08 

6755 

3.053 

33.40 

1.274 

0.10 

1556 

0.856 

45.12 

0.302 

0.10 

2032 

0.962 

40.73 

0.364 
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Table  5.7  (continued) 


0.10 

3947 

1.428 

39.29 

0.553 

0.10 

5309 

1.399 

38.87 

0.545 

0.12 

1486 

0.875 

45.58 

0.306 

0.12 

1895 

0.982 

45.28 

0.345 

0.12 

3556 

1.322 

40.75 

0.501 

0.12 

5117 

1.376 

39.50 

0.531 

Table  5.8  Shulze  Shear  Tester  results  for  saline  kao 

inite 

Fine 
Fraction 

Consolidation 
Stress 

(g) 

Unconfined 

Yield 
Strength,  cc 

(kPa) 

Angle  of 

Internal 

Friction,  <p, 

Shear  Stress 
(kPa) 

0.00 

2497 

1.791 

38.35 

0.702 

0.00 

4504 

2.198 

38.09 

0.865 

0.00 

6234 

2.193 

37.94 

0.865 

0.01 

2497 

1.791 

38.35 

0.702 

0.01 

4504 

2.198 

38.09 

0.865 

0.01 

6234 

2.193 

37.94 

0.865 

0.02 

2123 

0.651 

44.39 

0.233 

0.02 

3738 

0.789 

42.39 

0.291 

0.02 

5254 

0.926 

41.03 

0.349 

0.04 

1580 

0.792 

42.69 

0.291 

0.04 

1964 

0.877 

44.27 

0.314 

0.04 

3801 

0.935 

41.61 

0.350 

0.04 

5214 

0.968 

41.26 

0.364 

0.06 

1482 

0.894 

42.78 

0.328 

0.06 

1893 

1.039 

46.10 

0.360 

0.06 

3677 

1.078 

39.76 

0.414 

0.06 

5209 

1.203 

39.73 

0.463 

0.08 

1424 

1.283 

44.88 

0.454 

0.08 

1806 

1.300 

40.76 

0.492 

0.08 

3591 

1.452 

39.42 

0.561 

0.08 

5055 

1.525 

40.20 

0.582 

0.10 

1766 

1.436 

45.99 

0.499 

0.10 

3323 

2.424 

40.06 

0.928 

0.10 

4865 

2.915 

36.50 

1.172 

0.12 

1264 

1.638 

38.63 

0.640 

0.12 

1717 

1.816 

34.90 

0.745 

0.12 

3362 

2.291 

38.20 

0.900 

0.12 

4640 

2.964 

40.98 

1.119 
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Sand  2  with  Fresh  Kaolinite 
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Figure  5.15  Comparison  plot  of  the  shear  stress  at  incipient  flow,  rc,  for  the  Sand  #2  with 
fresh  kaolinite  series  of  beds. 
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Figure  5.16  Plot  of  the  friction  factor  at  a  consolidation  stress  of  1000  Pa  as  a  function  of 
\\i  for  the  Sand  #2  with  saline  kaolinite  series  of  beds. 
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Sand  2  with  Saline  Kaolinite 
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Figure  5.17  Comparison  plot  of  the  shear  stress  at  incipient  flow,  zc,  for  the  Sand  #2  with 
saline  kaolinite  series  of  beds. 


Sand  #2  and  Saline  Kaolinite  and  a  Normal  Load  of  1  Kg 
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Figure  5.18  Plot  of  the  friction  factor  at  a  consolidation  stress  of  1000  Pa  as  a  function  of 
\|/  for  the  Sand  #2  with  saline  kaolinite  series  of  beds. 
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5.4  Comparison  of  SERF  and  Flume  Data 

The  shear  strengths  clay  samples  used  in  both  the  SERF  and  the  flume  were 
plotted  versus  the  dry  density,  and  a  fit  of  the  following  form  (Mehta  and  Parchure,  2000) 
was  applied  to  each  clay  type  (Bennett  Pottery  #10,  Bennett  Pottery  #20,  Bennett  Pottery 

#75): 

r,=C,(pd-p,)c'  (5.3) 

where  C\  and  C2  are  empirical  coefficients,  Pd  is  the  dry  density  in  kg/m3,  and  p;  is  64 
kg/m  (selected  to  represent  the  dry  density  of  fluid  mud  with  no  shear  strength).  For  the 
#10  and  #20  clays,  Cx  =  3.20  x  10"3,  and  C2  =1,  and  for  #75,  d  =  1.30  x  10"3,  and  C2  =1 
were  obtained  (Figures  5.19  through  5.21).  Figure  5.22  is  a  comparison  plot  of  the 
results  for  all  three  clay  types.  It  was  necessary  to  create  beds  of  much  lower  density  for 
the  flume  than  those  used  in  the  SERF.  Yet  the  flume  results  agree  with  those  from  the 
SERF.  Consequently,  it  was  concluded  that  the  flume  was  accurately  measuring  the 
shear  stress  on  the  beds  placed  in  it. 
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Combined  Data  for  #10  Clay 
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Figure  5.19  Comparison  of  values  for  shear  strength  for  Clay  #10. 
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Figure  5.20  Comparison  of  values  for  shear  strength  for  Clay  #20. 
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Figure  5.21  Comparison  of  values  for  shear  strength  for  Clay  #75. 
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Figure  5.22  Summary  plot  of  the  comparison  of  the  shear  strength  values. 


CHAPTER  1 
ANALYSIS  OF  RESULTS 

6.1  Introduction 

In  Equation  3.3.34  for  the  critical  shear  stress   fcr,  or  simply  xcr,  of  sand 

influenced  by  the  presence  of  clay,  rcr  has  been  shown  to  vary  with  the  weight  fraction  v|/ 
so  that  in  the  low  range  of  \\i ,  rcr  tends  to  decrease  from  xco  for  pure  sand  to  r„  at  some 
V/  =  Vmin-  rcr  then  increases  with  a  further  increase  iny,  passing  through  the 
point  where  it  once  again  equals  rco  at  yr .  This  is  trend  illustrated  in  Figure  6.1.  Table 
6.1  summarizes  the  results  from  several  experiments. 
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Figure  6.1  Schematic  drawing  of  the  effect  of  clay  fraction  on  the  critical  shear  stress  for 
sand  (Data  of  Sharif,  2002). 
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Table  6.1  Summarized  data  for  rco,  v|/min ,  rcs,  and  v|/r 


Investigator(s) 

^50 

(mm) 

(Pa) 

Clay  type 

T  min 

Tcs 

(Pa) 

N>r 

Torfs  (1995) 

0.23 

0.35 

Scheldt  Muda 

0.018 

0.31 

0.04 

Torfs  (1995) 

0.23 

0.35 

Kaolinite 

0.018 

0.31 

0.04 

Torfs  (1995) 

0.23 

0.35 

Montmorillonite 

0.028 

0.33 

0.07 

Huygens  and 
Verhofen(1996) 

0.17 

0.18 

Kaolinite 

0.036 

0.15 

0.05 

Panagiotopoulos  et.  al. 
(1997) 

0.15 

0.32 

Combwich  mudb 

0.022 

0.30 

0.11 

Panagiotopoulos  et.  al. 
(1997) 

0.22 

0.35 

Combwich  mud 

0.015 

0.32 

0.23 

Sharif (2002) 

0.20 

0.22 

Kaolinite 
(1.27  cm  bed  depth) 

0.040 

0.18 

0.08 

Sharif (2002) 

0.20 

0.22 

Kaolinite 
(2.29  cm  bed  depth) 

0.040 

0.17 

0.12 

Present 

1.40 

0.99 

Kaolinite 

0.045 

0.87 

0.11 

Present 

0.83 

0.70 

Kaolinite 

0.060 

0.35 

0.13 

Present 

0.83 

0.70 

Kaolinite 
(saline) 

0.050 

0.42 

0.11 

Present 

0.83 

0.70 

Clay  Mixture  #2C 

0.030 

0.47 

0.06 

Present 

0.83 

0.70 

Clay  Mixture  #ld 

0.020 

0.50 

0.04 

Present 

a  .,_  .....     „„„ 

0.41 

0.35 

Kaolinite 

0.036 

0.20 

0.07 

55%  illite,  30%  smectite,  15%  kaolinite 

56%  illite,  26%  expansive  clays  (undefined),  1 1%  kaolinite,  and  7%  chlorite 

50%  kaolinite,  35%  attapulgite,  and  15%  bentonite 

50%  kaolinite,  45%  attapulgite,  and  5%  bentonite 


104 

This  work  examines  a  physics-based  explanation  of  the  behavior  of  rcr  under  the 
premise  that  the  interstitial  clay  particles  alter  grain-grain  friction  (and  interlocking) 
between  sand  particles,  and  that  the  nature  of  lubrication  is  a  function  ofy.  This 

explanation  is  accomplished  through  the  use  of  two  different  hypotheses.  In  the 
following  sections,  these  hypotheses  are  presented  first,  and  then  the  two  are  compared  to 
determine  their  relative  merits. 

6.2  Quasi-Hydrod v  namic  Model  for  Sand  Lubrication 

In  an  effort  to  explain  the  interaction  between  the  sand  grains  and  the  fine 
fraction,  a  mechanism  of  clay  flow  that  possesses  some  of  the  features  of  hydrodynamic 
lubrication  is  considered.  The  basic  feature  of  this  quasi-hydrodynamic  flow  is  an  in  situ 
layered  flow  of  a  clay  film  between  two  sand  grains,  portrayed  schematically  in  Figure 
6.2. 


Intermediate 
Film 


Vi 


-*■  v, 


± 


Figure  6.2  Schematics  of  the  layered  flow  of  a  clay  film  between  the  surfaces  of  two  sand 
grains  where  S{  and  5i  are  the  roughness  of  the  two  surfaces  (Heshmat  and  Brewe,  1995). 
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The  behavior  of  two  surfaces  that  begin  to  slide  with  respect  to  each  other  is 
strongly  modified  with  the  introduction  of  a  lubricant  between  them.  For  example,  if  for 
two  such  surfaces  we  plot  the  coefficient  of  friction  /  against  the  ratio  fiU/P,  where  fi  is 
the  lubricant  viscosity,  U  is  the  relative  velocity  of  the  one  surface  moving  against  the 
other,  and  P  is  the  specific  load,  we  find  three  different  lubrication  regimes  (Szeri,  1998). 
This  is  illustrated  in  Figure  6.3. 

At  large  values  of  piU/Pf  is  low  and  is  proportional  to  fiU/P.  This  is  the  regime 
of  thick-film  lubrication.  It  occurs  when  the  film  thickness  exceeds  roughly  2.5  jxm  and 
is  approximately  described  by  Petroffs  law  (Shigley  and  Mischke,  1989).  Within  this 
regime,  /  is  dependent  on  fx.  With  decreasing  \M/P,  f  passes  through  a  minimum,  as 
indicated  schematically  in  Figure  6.3,  and  the  behavior  enters  the  regime  of  thin-film 
lubrication.  In  this  regime,  the  film  becomes  thinner  and  some  of  the  opposing  asperities 
(roughnesses)  touch.  /  now  depends  on  the  surface  roughness,  the  properties  of  the 
solids,  and  the  properties  of  the  lubricant  .  With  a  decrease  in  /nU/P,  f  passes  into  the 
boundary  lubrication  regime,  in  which  the  film  is  so  thin  that  its  properties  are  no  longer 
the  same  as  those  of  the  bulk.  If  U  or  \i  is  reduced  or  P  increased,  the  film  becomes 
thinner  than  the  height  of  some  of  the  asperities.  If  these  asperities  are  covered  by  a 
suitable  molecular  layer  of  lubricant,  they  will  not  weld  together  (Szeri,  1998). 

The  model  of  clay  lubrication  in  Figure  6.2  permits  an  evaluation  of  the 
performance  of  clay  lubrication  in  quantitative  terms.  In  an  incipient  sliding  contact 
between  two  sand  grains  with  clay  lubrication  showing  an  infinitesimal  section  of  it  in  the 
sliding  direction,  the  film  can  be  divided  into  regions  of  intermediate  films  and  clay 
lubricant  film.    Intermediate  films  are  of  the  order  of  surface  roughness  elements,  the 
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asperities,  of  the  sand  grains.  In  the  present  study  it  was  observed  visually  that  the 
surfaces  of  the  sand  particles  in  the  beds  tested  were  covered  with  thin  adhered  layers  of 
clay  almost  immediately  upon  mixing  the  two  sediment  types  together. 

As  may  be  seen  in  the  schematic  model  of  Figure  6.2,  due  to  clay-induced 
lubrication  an  intermediate  film  acting  as  a  boundary  layer  is  created  between  the 
interacting  sand  grain  surfaces  and  the  clay  lubricant  film.  A  transition  between  the 
velocities  of  the  two  sand  grains  (Vi  and  V2,  where  the  subscript  "1"  represents  the  top 
sand  grain  and  subscript  "2"  represents  the  bottom  sand  grain)  and  the  velocities  of  the 
edge  of  the  lubricant  film  {V\  and  Uz)  occurs  across  this  intermediate  film.  This 
intermediate  layer  has  the  flexibility  of  behaving  either  as  a  solid  or  as  a  semi-solid  in 
which  creep  occurs,  which  makes  the  transition  of  flow  velocities  from  Vi  to  U\  and  V2  to 
U2  possible,  as  shown  in  Figure  6.2.  Figure  6.4  illustrates  this  quasi -hydrodynamic  flow. 


Thick-film  lubrication 


^ PetrofTs  law 


tiWP 


Figure  6.3  Lubrication  regimes. 
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Figure  6.4  Quasi-hydrodynamic  model  of  sand  lubrication  film. 

Hydrodynamic  lubrication  is  frequently  explained  through  the  application  of 
Couette  flow  to  a  slider  bearing.  Figure  6.5  is  an  illustration  of  a  simplified  slider 
bearing.  The  same  technique  may  be  used  here  if  two  sand  grains  are  represented  in  the 
analogous  fashion.  As  illustrated  in  Figure  6.6,  the  upper  sand  grain  can  be  shown  to  be 
the  slider  block,  and  the  lower  to  be  the  reference  surface  (or  pad). 
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Figure  6.5  Slider  bearing  moving  to  right  at  speed  U  above  stationary  pad  (a),  and 
transformed  to  stationary  block  and  pad  moving  to  the  left  at  speed  U  (b). 
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Tmbulent  flow r- 

Bed  plane 


Figure  6.6  Flow  of  interstitial  clay  suspension  viewed  as  slider  bearing. 


With  reference  to  Figure  6.6,  we  will  consider  the  problem  in  the  x-y  plane, 
treating  the  grains  as  if  they  were  of  unit  thickness  in  the  third  dimension.  Consider  that 
the  upper  sand  grain  is  set  in  motion  at  a  velocity  U  by  a  turbulence-induced 
instantaneous  shear  stress  at  the  exposed  bed  plane. 

We  may  conveniently  consider  the  short  duration  over  which  this  shear  acts  to  be 
long  enough  for  the  forces  to  be  in  equilibrium.  Einstein  and  El  Samni  (1949)  found  that 
the  pressure  fluctuations  (used  to  calculate  the  average  lift  on  an  individual  sand  grain  as 
discussed  in  Section  2.2)  had  a  predominant  period  of  approximately  l/50th  of  a  second. 
These  researchers  also  reported  the  lift  pressure  as  a  function  of  duration,  which  can  be 
used  to  determine  a  reference  instantaneous  lift  pressure.  A  lift  pressure  (10.05  kg/m2) 
occurring  only  3%  of  the  time  is  chosen  so  as  to  obtain  a  maximum  value  for  the 
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instantaneous  lift  pressure.  Through  Equations  3.3.8  and  3.3.12,  it  is  seen  that  both  the 
shear  stress  and  the  lift  per  unit  area  are  proportional  to  the  same  flow  properties  and 
should,  as  a  result,  be  of  the  same  order  of  magnitude.  Consequently,  this  maximum 
instantaneous  lift  pressure  may  be  used  to  determine  a  maximum  value  for  the 
instantaneous  load  on  a  single  0.83  mm  sand  grain  (found  to  be  0.022  g)  that  is  also 
representative  of  the  maximum  value  for  the  instantaneous  shear  stress  on  an  individual 
grain  of  the  same  diameter.  The  above  load  value  was  determined  as  follows: 

Instantaneous  Lift  Pressure  =  (0.0100584  m  water)(1000  kg/m3)  =  10.05  kg/m2 
Instantaneous  Load  =  (10.05  kg/m2)(7t)(0.00083  m)2  =  0.02 17g 

Now,  due  to  the  non-parallel  orientation  of  the  upper  and  the  lower  grain  surfaces 
separated  by  a  viscous  clay  suspension  of  small  thickness,  i.e.,  h{x)IL«\  where  h  is 
O(10  |im)  and  L  is  O(103  |im),  the  resulting  pressure  gradient  over  distance  L  will  cause 
the  suspension  to  flow.  If  we  choose  the  upper  grain  to  be  stationary,  the  lower  one  to 
translate  at  a  velocity  -U,  and  apply  Couette  flow,  then  for  small  values  of  angle  6  a 
velocity  profile  in  the  suspension  may  be  derived. 

Starting  with  the  continuity  equation  (considering  only  the  x-y  plane) 

du      dv 

T*+B-y=°  (621) 

and  using  the  following  dimensionless  variables, 

u  x 

u,=—  x,=  — 

U  I 

U  L                                           (6.2.2) 

v  v 

v.  =  —  y,  =  — 

V  y      h 

where  V  is  the  maximum  velocity  in  the  v-direction,  the  continuity  equation  can  be 
expressed  in  the  following  form 
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dx. 


fL\d^-o 


yh 


dy. 


(6.2.3) 


It  is  assumed  0  is  very  small,  or  hIL  is  very  small  (h«L).    Consequently,  for  the  two 
terms  of  Equation  6.2.3  to  be  of  the  same  order,  v  «  u  since  Uh»  1.  The  momentum 


equations  are: 


du      dp 
dx      dx 


fd2u    aV 

+ 


dx2 


of) 


•  g- 


(6.2.4) 


(6.2.5) 


The  ^-direction  momentum  equation  can  be  made  non-dimensional  as  follows 


du, 
ox. 


dp, 


( 


dx. 


n 


pUL 


d2u. 


dxl      \pVh\h)    dy1. 


(6.2.6) 


where  p,  =  p/p0  and  p0  is  the  ambient  pressure  in  the  clay  suspension.  Substituting  the 
definition  for  the  Reynolds  number  (Re  =  pUL/n)  into  Eq.  6.2.6  gives  the  following 


expression 


du, 
dx. 


rPo)dp.     (l)d2u, 


dx, 


vRe, 


dx2. 


+ 


(  i  Yrfd2u 


■U 


[Re^h)    dy2, 


(6.2.7) 


Since  Uh»l,  and  since  this  is  laminar  flow  (Re  <  104),  it  can  be  seen  that  the  first  and 
third  terms  of  Equation  6.2.7  will  be  much  smaller  than  the  rest.  The  first  term  will  be 
0(1)  and  the  third  term  0(1/Re),  whereas  the  second  term  0(pJL)  and  the  fourth  term  is 
0(L  Ih  ).  Consequently,  the  first  and  third  terms  can  be  neglected.  This  leads  to  a  thin 
layer  flow  approximation  in  which  gradients  across  the  flow  in  the  narrow  direction  are 
much  greater  than  gradients  along  the  flow: 


Ill 


dp  _     d2u 


(6.2.8) 


which  is  very  similar  to  Couette  flow.  The  only  difference  is  that  now  d p/dx  depends  on 
x  [since  u  is  f(jc,y)]. 

The  boundary  conditions  for  the  fluid  velocities  are: 


u=-U     @     y=0 
M=0        @     y=h 


and  when  applied  to  Equation  6.2.8  result  in 


u  =  -U 


h2  dpy 


( 


2pidx  h 


I--?- 


v        "J 


From  Equation  6.2.10  the  (short-duration  steady)  discharge  Q  is  given  by: 

Q  =  \udy 


or 


giving  for  dp/dx 


Q_-Uh      h3   dp 
2        Ylpidx 


dp     -\2fx(Uh      A 
ir  +  Q 


(6.2.9) 


(6.2.10) 


(6.2.11) 


(6.2.12) 


(6.2.13) 


dx         2 
The  boundary  conditions  for  the  pressure  are 

p  =  p0      @     x  =  0,L 
where  p0  is  the  pressure  in  the  clay  suspension  at  both  x  =  0  and  x  =  L.  Applying  these 
boundary  conditions  for  the  pressure,  an  expression  relating  U  and  Q  through  geometry  is 
obtained 
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Q  =  - 


rdx 

2.1*1 

2  rdx 


(6.2.14) 


Using  geometry  to  express  h  as  a  function  of  x  gives 


h  =  hl+(h2  -h^ 

Substituting  Equation  6.2.15  into  6.2.14  gives 


r,     h.  -h7 
tan0=-! 2- 


Q  = 


UK 


(6.2.15) 


(6.2.16) 


where 


h  =  2hih2 


°     ht+h2 


so  that 


dp  _     6/iU 
~dx"     ~tF 


MJ 


Substituting  Equation  6.2.18  into  6.2.10,  we  obtain 
u  =  -U 


k     h)    2fi   K   hK      °\     h 


which  may  be  simplified  to 


u  =  -U 


fl-l). 


hj 


^(h-hjh-y) 


or 


(6.2.17) 


(6.2.18) 


(6.2.19) 


(6.2.20) 


u  =  -U 


ri-Z).3Uy 


\    " ) 


+  ^l(h>-h0h-hy  +  h0y) 


(6.2.21) 


or 
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u  = 


Uy  |  3Uy    3Uyh0     3Uy 
h        h         h2  h2 


h       h 


2  +3Uh0y2 


h3 


(6.2.22) 


Substituting  Equation  6.2.16  into  6.2.22  gives 


u=    Ui4Uyi6Qy    3Uy2     6Qy2 


Ir 


h- 


ft3 


(6.2.23) 


which  can  be  seen  to  be 


u(y)=  -U  + 


(MJ     6Q\     (3U     6Q} 


■  + 
h      h2 


v- 


iz+'h: 


(6.2.24) 


In  obtaining  Equation  6.2.24  it  is  assumed  that  the  grain  asperities  are  of 
negligible  height,  i.e.,  grains  have  smooth  surfaces,  so  that  they  do  not  significantly  effect 
the  flow.  The  shear  force  F  on  the  lower  grain  is  obtained  from 


r,     f    du 


/* 


dx  =  - 

(k-\)hL 


41n£-6 


k-\ 


k  +  \ 


UL  =  K.UL       (6.2.25) 


where  //  is  the  viscosity  of  the  suspension,  k  =  h2lh\  and 


K,  = 


V 


(k  - 1\ 


41nA:-6 


k-\ 
k  +  1 


Hence  the  coefficient  of  friction /is  obtained  as 


'-H? 


(6.2.26) 


which  is  Petroff  s  law  for  thick-film  lubrication  (Figure  6.3). 

As  juU/P,  a  measure  of  horizontal  viscous  force  relative  to  normal  load,  decreases 
/will  decrease  until  h  approaches  O(10"1-10°)  urn,  the  size  of  clay  particles.  Hence/ will 
increase  with  a  further  decrease  in  ixU/P  in  the  zone  of  thin-film  lubrication.  Finally, 
there  will  be  a  further  and  rapid  rise  in  /  (boundary  lubrication)  as  the  film  becomes  so 
thin  that  its  properties  are  no  longer  those  of  the  bulk  and  abrasion/interlocking  due  to 
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clay  particles  and  asperities,  O(K)"1)  urn,  takes  over.    In  accordance  with  this  quasi- 
hydrodynamic  behavior  Equation  6.2.26  is  extended  to  yield 


( 


!*L\kSW 


f  =  Kx   !=—  +K2e     I'l  (6.2.27) 

in  which  K2  and  K3  characterize  thin-film  and  boundary  layer  lubrication.  Eq.  6.2.26  is 
recovered  when  K2  =  0. 

From  a  reanalysis  of  the  data  of  Heshmat  and  Brewe  (1996)  using  dry 
Molybdenum  sulfide  powder  as  a  lubricant  for  journal  bearings,  mean  values  were 
obtained  for  K\  =  0.0431,  K2  =  1.00  and  K3  =  0.31  -  2.2,  as  shown  in  Figure  6.7. 

As  previously  discussed  in  Sections  4.2,  4.5,  and  5.5,  46  different  sediment  beds 
were  subjected  to  increasing  water  flow  until  incipient  motion  of  sand  grains  occurred. 
As  a  specific  example,  in  runs  17-23  and  32-39  (listed  in  Table  5.3),  840  \ua  quartz  sand 
mixed  with  1  urn  kaolin  was  tested.  In  runs  17-23,  the  pore  water  was  fresh  and  in  runs 
32-39  it  was  pre-equilibrated  with  3  ppt  solution  of  NaCl  in  water.  A  salinity  of  3  ppt 
was  specifically  chosen  for  the  kaolinite  flocculation  properties  discussed  subsequently. 
A  comparison  plot  of  the  results  of  both  of  these  series  of  runs  may  be  seen  in  Figure  6.8. 
All  of  the  results,  which  are  shown  in  Figures  5.9  to  5.16  and  6.8,  seemingly  mimic  the 
behavior  of  Equation  6.2.27  which  is  illustrated  in  Figure  6.3. 
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Reanalysis  of  data  from  Heshmat  and  Brewe  (1996) 
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10.000  12.000  14.000 


Figure  6.7  Plot  of  Equation  6.2.27  applied  to  the  data  of  Heshmat  and  Brewe  (1996). 


Comparison  of  Fresh  and  Saline  Kaolinite  for  20-30  Sand 
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Figure  6.8  Comparison  plot  of  both  series  of  kaolinite  and  20-30  sand  tests. 


Referring  to  Equation  6.2.27,  from  Engelund  and  Zhaohui  (1984)  the  viscosity  of 
hyper-concentrated  flows  of  kaolinite  is  given  by 

H  =  0.001  +  0.206C1 68  (6.2.28) 


El. 

mc 

c  = 

.  vc 

+  vw 

_vc 

rsc 

v, 
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where  C  is  the  volume  concentration  of  the  suspension  expressed  as  a  dimensionless 
fraction,  and  p.  is  in  Ns/m2.  The  pore  water  kaolinite  suspension  concentrations  in  the 
present  study  ranged  from  0.00  to  0.37.  These  ranges  were  obtained  by  using  the 
following  expression 


(6.2.29) 


where  pd  is  the  dry  density,  psc  is  the  granular  density  of  the  clay,  mc  is  the  mass  of  the 
clay,  Vc  is  the  volume  of  the  clay,  and  Vw  is  the  volume  of  the  water.  The  concentration, 
C,  is  related  to  v|/,  where  \|/  is  defined  as  follows 

y  = S =     :       N^T (6-2-30) 

mt  +mc  +mw      P,V,  +  pK(l-/ic)Vc  +  pVw 

where  /ns  is  the  mass  of  the  sand,  Vs  is  the  volume  of  sand,  ps  is  the  granular  density  of 
the  sand,  Vw  is  the  volume  of  the  pore  water,  and  nc  is  the  porosity.  Solving  Equation 
6.2.30  for  Vc  results  in  the  following  expression 


w(p  V  +  pV  ) 
C"T77^         Vi \  (6.2.31) 

Substituting  Equation  6.2.31  into  6.2.29  results  in 

c-v(p.v,/vt+pVw/v.)  (ft0^ 

~       ~^~7\ Vi S —  (6.2.32) 

However,  if  it  is  assumed  within  the  range  of  fine  fraction  used  for  this  study,  that  the 
contribution  of  the  volume  of  clay  to  the  total  volume,  V,,  is  negligible,  the  volume  ratios 
can  be  refined  as  follows: 
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V  V  V 

— = ! =  — 5 —  =  1-", 

v,     vs  +  vc  +  vw     V,  +  vw 

(6.2.33) 


=  nr 


v,     vs+vc  +  vw    vs+vw 

where  ns  is  the  porosity  of  the  sand  grain  matrix.  This  assumption  is  shown  to  be  valid  in 
subsequent  sample  calculations.  If  these  ratios  are  substituted  back  into  Equation  6.2.32, 
the  following  expression  results  relating  C  to  vy: 

C  ~  ~Z1\ Y^ r  (6.2.34) 

which  is  valid  for  the  range  of  fine  fraction  used  in  this  study.  Assuming  ps  =  psc  = 
2.65x1 0"6  kg/mm3,  p  =  l.OOxlO"6  kg/mm3,  nc  =  0.48  (this  value  is  discussed 
subsequently),  ns  =  0.50  (taken  from  Hough,  1957,  for  clean,  uniform  surface  sand),  and 
v|/  =  0.16  (the  upper  range  of  fine  fraction  used  in  this  study),  Equation  6.2.34  gives  C  = 
0.37  as  the  upper  range  of  concentration  for  this  work.  The  viscosity  consequently  varied 
as  follows: 

\i  =  0.001  +  0.206(0.0)' 68  =  0.001  , 

v     '  Ns/m2 

H  =  0.001  +  0.206(0.37)' 68  =  0.040 
The  flow-induced  drag  coefficient  is 

Cd=^y  (6.2.35) 

P"cr 

where  p  is  water  density  and  ucr  is  the  critical  velocity.  If  we  consider  ucr  to  be  invariant, 

then  since  /  can  be  considered  to  be  related  to  CD,  the  variation  of  /  can  in  turn  be 

considered  to  qualitatively  reflect  the  variation  of  xcr.     By  substituting  Equation  6.2.34 

into  Equation  6.2.29,  viscosity  n  is  seen  to  vary  with  \|/  as  follows 
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j«-  0.001  +  0.206{k(lrns)vP"s]T 


1.68 


Thus  the  behavior  of  \x,  and  consequently  //£//P,  can  be  considered  to  reflect  the  behavior 
of  \|/.  Therefore,  the  variation  of/ with  fiU/P  is  seen  to  be  qualitatively  analogous  to  the 
variation  of  xcr  with  y . 

6.3  Effect  of  Pore  Water  Salinity 

From  Figure  6.8  we  note  that  rco  =  0.70  Pa.  Further,  for  saline  water  xcs  =  0.42  Pa 
at  Ymin  ~  0.05,  and  for  fresh  water  xcs  =  0.35  Pa  also  at  v|/min  ~  0.06.  From  the 
corresponding  r„/rco  =  0.6  and  0.5,  we  may  infer  that  kaolin  in  fresh  water  provides  -17% 
greater  lubrication  than  in  saline  water,  and  that  \|/  =  0.06  marks  the  transition  between 
thin-film  and  thick-film  lubrication.  Both  of  these  trends,  the  location  of  i|/min  and  the 
reduction  in  tcs  with  fresh  pore  water,  are  explained  by  the  quasi-hydrodynamic  model  of 
clay  lubrication. 

The  lubricating  effect  of  kaolinite  in  pore  water  is  determined  by  the  orientation 
of  the  individual  plate-like  particles.  Kaolinite  is  characterized  by  particles  of  relatively 
large  size  and  a  low  cation  exchange  capacity  in  comparison  to  other  clay  minerals.  For 
example,  the  length  and  thickness  of  typical  kaolinite  particles  are  1  (am  and  0.1  um, 
respectively,  compared  with  equivalent  dimensions  of  0.1  and  0.01  um,  respectively,  for 
montmorillonite. 

Clay  particles  have  a  negative  charge  on  the  surface  and  a  positive  charge  on  the 
edges.  Hence,  certain  clays  can  be  flocculated  by  edge  to  face  mechanisms  (Figure  6.9) 
under  a  state  of  non-salt  flocculation.    The  addition  of  salt  in  the  suspending  medium 
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increases  the  concentration  of  negatively  charged  ions  which  gather  at  the  edges  of  the 
particles  and  thus  make  the  edges  less  effective  in  flocculation.  Consequently,  a  greater 
degree  of  parallelism  results  with  salt-flocculation  (Figure  6.10).  The  structure  of  clay 
with  edge-to-face  contacts  of  particles  is  more  open  and  has  lower  shear  resistance  than 
that  of  a  more  compact  structure  with  face-to-face  contact  of  particles  (Lambe  and 
Whitman,  1969). 


Figure  6.9  Freshwater  flocculation  in  pore  water  with  edge-to-face  contacts. 


Figure  6.10  Saltwater  flocculation  with  particle  orientation  approaching  parallelism. 
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Both  of  these  types  of  flocculation  are  observed  in  the  case  of  kaolinite  in  water  at 
different  salt  concentrations.  Kaolinite  in  fresh  water  flocculates  by  the  edge-to-face 
mechanism,  as  illustrated  in  Figure  6.9,  resulting  in  a  very  weak  particle  matrix.  The 
presence  of  concentrations  of  salt  in  the  suspending  fluid  (greater  than  approximately  1 
ppt),  like  the  concentration  (3  ppt)  used  in  this  study,  results  in  salt-flocculation  of  the 
type  illustrated  in  Figure  6.10.  In  this  case  a  matrix  of  lower  density  and  shear  strength 
than  the  low-salt  case  results,  but  still  having  greater  density  and  shear  strength  than  with 
fresh  water  flocculation  (Lambe  and  Whitman,  1969). 

The  addition  of  salt  at  3  ppt  concentration  resulted  in  a  matrix  that  was  both 
denser  and  stronger,  and  as  a  result  had  a  viscosity  higher  than  the  case  with  fresh  water 
pore  water.  From  Wasp  et  al.  (1977),  the  effective  viscosity  of  Bingham  plastic  fluid  is 
given  as 

where  T0  is  the  yield  stress,  Dp  is  the  pipe  diameter,  and  U,  is  the  transition  velocity.  This 
expression  may  be  used  to  determine  the  ratio  of  the  effective  viscosity  of  saline  kaolinite 
to  the  effective  viscosity  of  fresh  kaolinite 

r~es    _      os 


(6.3.2) 


where  the  subscript  "s"  represents  saline  kaolinite,  and  the  subscript "/'  represents  fresh 
kaolinite.  From  the  work  of  Einstein  (1941),  this  ratio  is  found  to  be  40:1.  Equation 
6.3.32  can  then  be  combined  with  Equation  6.2.28  to  determine  the  viscosity  of  saline 
kaolinite  as  follows 
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A*«  = 


(o.OOl  +  0.206C168)=  (40)[o.001  +  0.206(0.37)' 68]=  1.6  Ns/m2    (6.3.3) 


V 


(The  quasi-hydrodynamic  model  then  would  predict  this  bed  to  have  a  higher  friction 
factor,  since  in  the  thick-film  region  of  lubrication  Equation  6.2.27  is  essentially 
governed  by  Petroffs  law  where  the  friction  factor  is  directly  proportional  to  the 
viscosity.  With  an  increase  in  viscosity,  as  in  the  comparison  of  the  fresh  water  kaolinite 
to  saline  kaolinite  illustrated  in  Figure  6.8,  Equation  6.2.27  would  result  in  an  increase  in 
friction  factor. 


Figure  6. 1 1  Face-centered  cubic  unit  cell. 


If  we  assume  that  sand  grains  are  hexagonally  close-packed  in  the  bed,  we  can 
determine  both  the  number  of  layers  of  clay  particles  and  the  total  thickness  of  the  layers 
within  the  pore  water  of  the  resulting  face-centered  unit  cube  (Figure  6.11)  that  will 
provide  lubrication  between  sand  grains  at  \|/min.  As  may  be  seen  in  Figure  6.11,  in  this 
structure  the  grains  touch  along  the  diagonals  for  each  face  (and  not  along  the  edges  of 
the  cube).  The  length  of  the  diagonal  is  rg  +  2rg  +  rg,  or  4rg,  where  rg  is  the  sand  grain 
radius.  Hence  the  edge  length  is  rg  V8  .  Since  the  net  number  of  sand  grains  in  the  face- 
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centered  cubic  unit  cell  is  4,  the  total  volume  occupied  by  the  pore  fluid  is  the  total 
volume   of  the   cell  [rg  V8  J  minus   the   volume   occupied   by   the   four   sand   grains 


4(4/3^),  or  [(rgJ%J  - 4(4/3^ )' 


A  face-centered  cubic  structure  with  a  sand  grain  having  a  diameter  of  0.83  mm,  a 
clay  particle  size  of  lu,m,  and  a  fine  grain  fraction  of  0.066,  will  have  a  clay  layer 
approximately  20  p.m  thick  and  will  consequently  have  approximately  20  layers  of  clay 
particles.  This  calculation  is  given  below: 

Vtotal  =  (0.415V8)3  =  1.62mm3 
1A 


V      =4 

Ysand 


-  W(0.415)3  =  1.20mm3 


f4 
V      =4 

vclay 


-  W[(0.415  +  0.02 )3  -  (0.415)3]=  0.182mm3 


V      =V      -V      -V      =  0  24mm3 

vpore         v  total         vsand         vclay        U.^MIIIU 

= 2.65x10^(1 -0.48)Vclay 

V     2.65xl0"6 Vsand  +  2.65xl0"6(l  -  0.48)Vclay  +  1.00xlO-6(0.48)Vclay  +  l.OOxlO"6 V^  "  0"°66 


assuming  pchy  =  psand  =  2.65x1 0"6  kg/mm3,  Pw  =1.00xl0"6  kg/mm3,  the  porosity  of  the  clay 
=  0.48: 

n   _      PbV,o,  -me  (l850kg/m3)(0.005984m3)- 10.001kg 

C     P(VM-mJpJ     l(l000kg/m3)(0.005984m3  - (l0.001kg)/(2650kg/m3))J      ' 

(where  Vt0,  =  Vc  +  Vw),  and  the  clay  bonded  to  the  sand  grains  in  layers.  This  last 
assumption  is  based  on  the  observations  of  Torfs  (1995),  Pangiatopolous  et.  al.  (1997), 
and  visual  observation  made  during  the  laboratory  portion  of  the  present  study.  Torfs  and 
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Panagiatopolous  both  commented  on  the  homogeneity  of  the  sand  and  clay  mixtures  used 
in  their  work.  In  this  study  it  was  observed,  especially  in  the  sand  beds  with  very  low 
fine  fractions  (0.2  and  0.4),  that  the  beds  looked  to  be  made  up  of  "dirty"  or  discolored 
sand.  In  other  words,  the  clay  could  be  seen  with  the  naked  eye  to  bond  to  the  sand 
grains  giving  the  beds  the  appearance  that  they  were  made  up  of  a  different  colored  sand. 

It  is  interesting  to  note  that  the  thickness  of  the  clay  layer  (20  \xm)  at  the 
minimum  point  of  critical  shear  stress  (r„)  is  close  to  the  size  of  the  asperities  (~  20  u.m) 
for  a  "rough"  quartz  surface  (prepared  using  a  No.  220  grit  diamond  wheel)  as  reported 
by  Lambe  and  Whitman  (1969)  and  as  seen  in  Figure  4.21.  This  is  precisely  the  point  at 
which  the  transition  from  thin-film  to  boundary  layer  lubrication  occurs. 

A  face-centered  cubic  structure  with  a  sand  grain  having  a  diameter  of  0.83  mm,  a 
clay  particle  size  of  l^im,  and  a  fine  grain  fraction  of  0.13,  will  have  a  clay  layer 
approximately  40  |im  thick  and  will  consequently  have  approximately  94%  of  the 
interstitial  pore  space  filled  with  clay.  This  calculation  is  given  below: 

Vlotal  =  (o.415>/87  =  1.62mm3 


V      =4 

vsand 


-W(0.415)3=  1.20mm3 


-W  [(0.4 15  +  0.04)3  -  (0.415)3]=  0.38mm3 


V      =4 

vclay 


2.65x10-*  (l-0.48)Vclay 

V "  2.65x10"^^  +2.65xl0-6(l-0.48)Vclay  +1.00xlO-6(0.48)Vday  +1.00xlO-6Vpore 

V 

_*L  =  0.94 

V 

pore 
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In  this  case,  it  is  interesting  to  note  that  "space-filling"  is  occurring  at  a  fine  grain 
fraction  of  0.13.  This  is  an  important  point,  since  it  implies  that  v|/r  (~  0.13  for  fresh 
kaolinite  as  shown  in  Figure  6.8),  which  is  the  point  that  the  shear  stress  once  again 
equals  xco  and  the  sediment  starts  to  be  dominated  by  the  clay  component  (as  can  be  seen 
in  Figure  6.8),  appears  to  be  the  space-filling  fine  fraction. 

6.4  Shear  Resistance  Model  for  Clay  Lubrication 

Relative  sliding  between  particles,  which  is  the  main  mechanism  causing 
deformation  within  the  soil  mass  including  at  the  sand  bed  surface,  is  determined  largely 
by  shear  resistance  at  contacts  between  particles.  Considering  the  force  balance  of 
Chapter  3  (Equation  3.3.2),  this  model  seeks  to  explain  the  behavior  of  rcr  through  the 
change  in  shear  resistance  that  accompanies  a  change  in  y . 

For  inter-particle  shear  resistance,  which  is  frictional  in  nature,  there  are  two  basic 
laws  (Amonton's  laws):  (1)  shear  resistance  between  two  bodies  is  proportional  to  the 
normal  force  between  the  bodies,  and  (2)  shear  resistance  between  two  bodies  is 
independent  of  the  dimensions  of  the  two  bodies.  According  to  Lambe  and  Whitman 
(1969): 

(1)  On  a  submicroscopic  scale  most  surfaces  are  actually  rough,  and  consequently 
two  solids  will  be  in  contact  only  where  the  asperities  touch  one  another;  i.e.,  actual 
contact  is  a  very  small  fraction  of  the  apparent  contact  area  as  may  be  seen  in  Figure 
6.12. 
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(2)  Because  contact  occurs  at  discrete  sites,  the  normal  stresses  across  these 

contacts  are  high,  and  even  under  light  loading  reach  the  yield  strength  of  the  material  at 

these  sites.  Therefore,  the  actual  area  of  grain-grain  contact  Ag  is 

N 
Ag=—  (6.4.1) 

where  N ,  which  is  equivalent  to  P  of  the  quasi-hydrodynamic  model,  is  the  normal  load 
and  qu  is  the  normal  stress  required  to  cause  yield,  i.e.,  plastic  flow.  Since  qu  is  fixed  in 
magnitude,  an  increase  in  the  total  normal  load  between  the  two  surfaces  implies  a 
proportional  increase  in  the  area  of  actual  contact.  This  increase  in  area  is  a  result  of 
plastic  flow  of  the  asperities. 

(3)  The  high  contact  stresses  cause  the  two  surfaces  to  adhere  at  the  points  of 
actual  contact;  i.e.,  the  two  bodies  are  joined  by  chemical  bonds.  Shear  resistance  is 
provided  by  the  adhesive  strength  of  these  points.    The  maximum  possible  shear  force 

*max  IS 

r«  =  sAg  (6.4.2) 

where  s  is  the  shear  strength  of  the  adhered  junctions.  Combining  Equations  6.4.1  and 
6.4.2  leads  to 

Tmm=N-  (6.4.3) 

Since  s  and  qu  are  material  properties,  Tmax  is  proportional  to  N. 
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Figure  6.12  Microscopic  view  of  frictional  resistance. 


While  Equation  6.4.3  is  appropriate  for  sand  beds,  it  needs  to  be  modified  to 
include  the  effect  of  clay  in  pore  water.  Clay  surfaces  have  irregularities  not  in  the  form 
of  asperities,  but  in  the  form  of  step-like  mesas  and  plateaus,  and  this  fact  combined  with 
the  size  of  the  particles  results  in  a  different  shear  resistance  behavior  than  that  of  sand 
grains.  On  clay  surfaces,  the  steps  are  on  the  order  of  10  -  100  A  compared  to  an  asperity 
size  of  about  1000  -  10,000  A  for  sand  grains.  Even  when  clay  surfaces  come  into  close 
proximity  over  almost  their  entire  area,  they  may  not  come  into  direct  contact.  Foreign 
molecules  on  the  surfaces,  including  sorbed  water,  are  not  squeezed  out  from  between  the 
surfaces  unless  the  normal  stress  exceeds  approximately  5.50  x  108  Pa.    Rather,  these 
foreign  molecules  participate  in  the  transmission  of  the  normal  stress.  Hence  with  clays, 
the  normal  stress  is  carried  over  much  greater  areas  and  is  transmitted  throughout  the  clay 
matrix  rather  than  over  a  few  asperities,  as  in  the  case  with  sand  grains  resulting  in  plastic 
deformation.    This  is  the  case  even  with  fresh  water  flocculation  (Figure  6.9),  which 
results  in  an  edge-to-face  orientation  that  is  more  nearly  similar  to  the  asperity  contacts 
discussed  for  granular  particles,  except  that  in  the  case  of  clays  each  contact  probably 
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contains  only  one  asperity.  However,  because  of  the  small  size  of  the  clay  particles,  the 
normal  load  is  spread  across  several-fold  more  asperities  than  with  the  case  of  sand 
grains.  Consequently,  the  load  on  each  asperity  is  orders  of  magnitude  smaller  than  with 
the  case  of  sand  grains.  This  is  best  illustrated  with  an  example. 

A  sand  grain  contact  is  illustrated  in  Figure  6.13(a).  The  grains  are  assumed  to  be 
spheres  and  the  diameter  of  the  contact  area  between  the  two  grains  is  given  by  (Lambe 
and  Whitman,  1969) 


d„  = 


12(l-v2) 

E  /?,  +  fl2 


lN-Afr 


(6.4.4) 


which  is  simply  an  equation  for  the  diameter  of  the  contact  area  between  two  elastic 
spheres  as  derived  by  Lambe  and  Whitman  (1969). 


0.02  mm 


(a) 


(b) 


Figure  6. 13  Sand  grain  arrangement  for  calculation  of  shear  resistance  example. 
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In  the  above  equation,  the  Poisson's  ratio  v  =  0.31  for  quartz,  Young's  modulus  E  =  7.58 
x  1010  Pa  for  quartz,  and  Ri  and  R2  are  the  grain  radii  (which  may  be  assumed  equal). 
This  equation  is  combined  with  Equation  6.4.1  to  determine  the  average  contact  stress 

(qavg  =  N/A)  between  two  quartz  spheres: 

1 

N        4AT3 

<U=7  = r  (6-4-5) 

nC]Rl 

where  C\  =  0.0004  (1/Pa)1/3.  If  we  assume  the  grains  are  oriented  in  a  cubic  array  such  as 
in  Figure  6.13(b),  we  need  only  consider  a  horizontal  plane  through  the  unit  cell  itself. 
From  Figure  6.13  we  note  that  there  will  be  four  contacts  in  a  unit  cell,  and  consequently 
the  average  contact  stress  will  be  one  fourth  of  the  average  contact  stress  for  a  single 
contact 


1 

N  N> 


*       A  - 

""     tcC?R~3 


(6.4.6) 


where  Atol  is  the  total  contact  area  for  the  unit  cell  equal  to  four  times  the  area  of  a  single 
contact.  Given  that  plastic  yield  for  quartz  occurs  at  a  normal  stress  of  1.03xl010  Pa 
(Lambe  and  Whitman,  1969),  Equation  6.4.6  can  be  used  to  determine  the  normal  load 
and  from  it  the  confining  stress  needed  to  cause  plastic  yield: 


(  *V    r.  .       .  ,2^3 


N  = 


9«*C*R> 


v 


(l.03xl010  Pa^C^l.OxlO"5  mf 


=  721g         (6.4.7) 


q  =  JL=  — 721g        =1.8xl0'2Pa  (6.4.8) 

A*u      (2.0xl0"5  mf 

where  q  is  the  over-all  confining  stress  applied  over  the  area  of  the  unit  cell,  Act\\. 
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If  we  now  consider  that  each  sphere  is  not  perfectly  smooth  but  contacts  its 
neighbors  on  one  asperity  with  a  tip  radius  of  1000  A,  Equation  6.4.6  can  again  be  used 
to  determine  the  normal  load  and  the  confining  stress  that  must  be  applied  to  cause  plastic 
yield: 


yv  = 


»Y 


qavsc2xw 


<7  = 


N 


2 

(l  .03x10'°  Pa)rC2  (l.OxlO-7  mfi 

33.2g 


4*u      (2.0x1 0~5  mf 


=  8.3xl0luPa 


=  33.2g         (6.4.9) 


(6.4.10) 


If  we  consider  the  dimensions  of  kaolinite,  l^im  long  x  0.1  |im  wide  x  0.01  ^m 
thick,  and  assume  the  porosity  (nc)  to  be  0.48,  we  can  determine  the  number  of  particles 
per  unit  volume 

Vsolid  =(!-">,=  0.52m3  (6.4.11) 


#  particles  = 


solid 


0.52  m3 


Vparticle      1.0xl0-21m3 


=  5.2x10 


20 


(6.4.12) 


The  number  of  particles  contacting  a  unit  area  may  be  assumed  equal  to  (Lambe  and 
Whitman,  1969) 


Number  Particles 

Number  Particles > 

2 
3 

Unit  Area 

,      Unit  Volume     i 

Number  Particles 
Unit  Area 

(5 

.2xl020^ 
lm3 

2 
.1 

=  6.47x 

1013 

(6.4.13) 


(6.4.14) 


If  we  consider  that  each  particle  makes  two  stress-carrying  contacts  on  a 
horizontal  plane,  we  find  that  the  number  of  contacts  per  unit  area  is  twice  the  number  of 
particles  per  unit  area,  or  1.29  x  1014. 
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If  we  further  consider  each  particle-particle  contact  with  the  edge  to  have  a  radius 
of  0.01  urn  and  that  plastic  flow  starts  at  6.9x1 09  Pa  (Lambe  and  Whitman,  1969),  these 
values  can  be  used  with  Equation  6.4.6,  to  determine  the  confining  stress  that  must  be 
applied  to  cause  plastic  yield 


N  = 


(6.9xl09Pa)rC12(l.29xl014Xl.0xl0-8mJ 


3 

20 


=  1.25xl0^g        (6.4.15) 


,.f.,""B''.ija,lft>  (6.4.16) 

A*n      (2.0xl0"5  mjr 

From  this  example  it  is  seen  that  the  normal  load  at  a  discrete  contact,  N,  required 
to  cause  plastic  yield  increases  significantly  from  the  case  where  the  overall  confining 
stress  is  carried  by  a  sand  grain  matrix  in  which  the  grains  are  modeled  as  spheres  with 
1000  A  asperities  (N  =  33.2  g  as  shown  in  Equation  6.4.9)  to  the  case  where  the  stress  is 
carried  by  a  kaolinite  suspension  with  a  porosity  of  0.48  (N  =  1.25  x  1020  g  as  shown  in 
Equation  6.4.15). 

Since  the  friction  factor  (f=  TmJN)  is  inversely  proportional  to  the  normal  load, 
this  helps  to  explain  the  behavior  illustrated  in  Figure  6.14.  When  v|/=  0.0,  the  confining 
stress,  due  to  a  normal  load  is  carried  completely  by  the  sand  matrix.  The  normal  load  at 
each  discrete  contact,  as  a  result,  is  low  (~  33  g)  and  the  friction  factor  is  high  (f=  0.68). 

When  v)/>\j/min ,  the  confining  stress  is  carried  completely  by  the  clay  matrix.  The 
normal  load  at  each  discrete  contact  is  much  higher  (-1.25  x  1020  g)  and  the  friction 
factor,  consequently,  much  lower  (f~  0.3). 

At  low  values  of  v(/(<\|/min ),  when  assuming  a  face-centered  cubic  structure  for 
the  sand  matrix  it  was  shown  in  Section  6.3  that  the  clay  layer  is  not  thick  enough  to 
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cover  the  asperities  of  the  sand.  As  mentioned  previously,  a  sand  grain  having  a  diameter 
of  0.83  mm,  a  clay  particle  size  of  lu,m,  and  a  fine  grain  fraction  of  0.066,  will  have  a 
clay  layer  approximately  20  ^m  thick  which  is  approximately  the  height  of  the  asperities 
on  rough  quartz.  Consequently,  when  v|/<ymin  the  over-all  confining  stress  is  carried  by 
both  the  sand  matrix  and  the  clay  suspension  in  between  the  sand  grains. 

As  f  increases  from  zero,  more  clay  particles  are  added  to  the  sand  matrix  and  the 
clay  increasingly  shares  the  load  carried  by  the  sand  particles,  resulting  in  a  decrease  in 
the  value  of  N  for  the  sand  grains.  Considering  Equations  6.4.1  and  6.4.2,  it  can  be  seen 
that  this  decrease  in  the  normal  load  carried  by  the  sand  matrix  will  result  in  a  decrease  in 
the  contact  area  between  the  sand  grains,  Ag,  and  a  corresponding  decrease  in  the 
maximum  possible  shear  force,  T^,  for  the  sand,  which  will  in  turn  result  in  a  decrease 
in  the  friction  factor  for  the  sand  grains  (f=  0.15). 

With  fresh  water  flocculation  for  kaolinite  there  is  an  increased  number  of  edge- 
to-face  contacts  than  with  salt  water  flocculation,  which  results  in  a  lower  friction  factor 
for  the  sediment  beds  with  fresh  water  kaolinite.  From  Equation  6.4.15  it  can  be  seen 
that  the  normal  load  is  directly  proportional  to  the  number  of  contacts.  Consequently,  as 
shown  in  Figure  6.14,  for  a  given  value  of  v|/  the  friction  factor  for  the  fresh  kaolinite  will 
be  less  than  the  friction  factor  for  the  saline  kaolinite.  Also,  this  should  lead  to  a  greater 
difference  between  the  two  beds  as  \|/  increases,  because  the  clay  suspension  will  have 
more  of  an  influence  on  the  behavior  of  the  friction  factor  as  it  will  carry  more  of  the 
over-all  confining  stress. 
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♦  Saline  Kaolinite  ■  Fresh  Kaolirite 


Figure  6.14  Plot  comparing  the  results  from  the  Schulze  cell  for  the  saline  kaolinite  and 
fresh  kaolinite  series  of  sand  beds. 


6.5  Comparison  of  Model  Concepts 

The  two  models  are  compared  in  Table  6.2,  their  strengths  and  weaknesses  listed, 
to  determine  their  relative  merits. 


Model 


Quasi- 

Hydrodynamic 

Model 


Table  6.2  Comparison  of  mode 


Strengths 


Effective  in 

quantitatively  predicting  the 
behavior  over  the  entire 
range  of  lubricating 
thickness. 

Effective  in  explaining 
the  behavior  of  the  mixture 
once  the  sand  grains  are 
separated  enough  to  have  a 
clay  layer  between  them. 


concepts 


Weaknesses 


Focus  is  on  the  thin  film 
and  thick  film  regions  of 
lubrication. 

Does  not  provide  a 
detailed  explanation  for  the 
observed  behavior  in  the 
boundary  layer  region  of 
lubrication. 
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Table  6.2  (Continued) 

Model 

Strengths 

Weaknesses 

•           Effective  in  explaining 

•           Focus       is       on      the 

the  behavior  of  the  mixture 

boundary    layer    region    of 

when  the  sand  grains  are  in 

lubrication. 

Shear  Resistance 

contact  or  slightly  separated. 

•           Does     not     provide     a 

Model 

detailed  explanation  for  the 
observed  behavior  in  the  thin 
film  and  thick  film  regions  of 
lubrication. 

From  Figure  6.8  it  is  seen  that  the  quasi-hydrodynamic  model  predicts  a  value  for 
\|/min  for  kaolinite  that  is  between  0.04  and  0.06.  This  value  matches  the  y^  values 
obtained  from  the  data  for  kaolinite  beds  of  other  researchers;  the  data  from  Torfs  (1995) 
results  in  y^  =0.02,  the  data  from  Huygens  and  Verhofen  (1996)  results  in 
^min=  0.04 ,  and  the  data  from  Sharif  (2002)  results  in  \|/min  =  0.04 . 

It  is  also  noted  that  assuming  a  face-centered  cubic  structure  for  the  sand  grain 
matrix,  \^/min  is  the  fine  weight  fraction  at  which  the  thickness  of  the  clay  layer  is  close  to 
the  size  of  the  asperities  of  the  sand  grains  (~  20  urn)  at  which  the  transition  from  thin- 
film  to  boundary  layer  lubrication  occurs. 

From  Figure  6.14  it  is  seen  that  the  shear  resistance  model  predicts  y^  =  0.01  for 
kaolinite.  This  value  does  not  match  well  with  the  v|/min  values  obtained  from  the  data  for 
kaolinite  beds  of  other  researchers  0.02-0.04.  In  addition,  this  value  for  y^  results  in  a 
clay  thickness  much  less  than  the  size  of  the  asperities  of  the  sand  grains  (~  20  nm),  as  a 
sand  grain  having  a  diameter  of  0.83  mm,  a  clay  particle  size  of  1  ^im,  and  a  fine  grain 
fraction  of  0.01,  will  have  a  clay  layer  approximately  3  ^m  thick  assuming  the  sand 
grains  are  organized  in  a  face-centered  cubic  structure.  It  is  also  noted  that  the  minimum 
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friction  factor  (~  0.15)  occurs  when  the  thickness  of  the  clay  layer  is  approximately  3  |im, 
from  which  it  may  be  concluded  that  the  strongest  lubricating  effect  occurs  when  the 
over-all  confining  stress  is  carried  by  both  the  sand  matrix  and  a  thin  clay  suspension  in 
between  the  sand  grains. 

Finally,  it  should  be  noted  that  v|/min  was  not  found  to  represent  the  space-filling 
weight  fraction,  contrary  to  the  suggestion  by  Torfs  et  al.  (2000).  Assuming  that  sand 
grains  are  organized  in  a  face-centered  cubic  structure,  the  sand  grain  diameter  is  0.83 
mm,  and  the  clay  particle  size  is  1  |im,  it  is  not  until  the  fine  grain  fraction  reaches  a 
value  greater  than  0.13  that  the  pore  volume  not  occupied  by  clay  becomes  negligible  (< 
1%  of  the  total  volume).  In  other  words,  it  is  not  until  a  fine  fraction  greater  than  0.13  is 
reached  that  the  clay  suspension  fills  the  entire  volume  of  pore  space  within  the  sand  unit 
cell.  This  is  an  important  point,  since  it  implies  that  \|/r  (~  0.13  for  fresh  kaolinite  as 
shown  in  Figure  6.8),  which  is  the  point  that  the  shear  stress  once  again  equals  tco  and  the 
sediment  starts  to  be  dominated  by  the  cohesive  component  (as  can  be  seen  in  Figure 
6.8),  appears  to  be  the  space-filling  fine  weight  fraction. 


CHAPTER  7 
SUMMARY  AND  CONCLUSIONS 

7.1  Summary 

It  was  the  goal  of  this  study  to  develop  an  explanation  for  the  effect  of  clay 
particles  on  the  critical  shear  stress  tcr  for  sand.  To  meet  this  goal  the  following  tasks 
were  carried  out: 

1 .  An  expression,  based  on  a  force-balance  for  a  sand  grain  at  the  surface  of  a  bed,  was 
developed  for  the  critical  shear  stress  of  sand  beds  in  the  presence  of  interstitial  clay 
particles,  rcr,  in  the  range  of  fine  fraction,  \j/,  from  0.0  to  0.15. 

2.  A  narrow,  unidirectional  flow  flume  with  a  false  bottom  designed  for  emplacing 
sediments  beds,  was  used  to  determine  the  relationship  of  the  critical  shear  stress  of  the 
sand  beds  as  a  function  of  clay  weight  fraction. 

3.  The  expression  for  the  critical  shear  stress  of  sand  beds  in  the  presence  of  interstitial 
clay  particles  was  applied  to  the  results  from  the  flume  tests  to  determine  the  effect  of  the 
fine  grain  fraction  on  the  effective  particle  packing  angle,  Oeff- 

4.  A  quasi-hydrodynamic  model  of  clay  lubrication,  a  mechanism  of  interstitial  clay  flow 
which  possesses  some  features  of  hydrodynamic  lubrication,  was  developed  to  explain 
the  interaction  between  sand  grains  and  the  fine  fraction. 

5.  A  Schulze  ring  shear  tester  was  modified  for  use  with  clays  and  then  used  to  explain 
the  behavior  of  xcr  of  sand  at  low  values  of  \\f. 

6.  A  shear  resistance  model,  which  relates  to  the  nature  of  shear  resistance  between 
particles,  was  developed  to  explain  the  dependence  of  rcr  on  the  fine  grain  fraction. 
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7.2  Conclusions 

The  following  are  the  main  conclusions  of  this  study: 

1.  From  the  application  of  the  expression  for  the  critical  shear  stress  rcr,  of  sand 
influenced  by  the  presence  of  clay,  xcr  was  shown  to  vary  with  the  weight  fraction  y  so 
that  in  the  low  range  of  y ,  zcr,  tended  to  decrease  from  xco  for  pure  sand  to  tcs  at  some 
y-Ymin-  Tcr  then  increased  with  a  further  increase  iny,  passing  through  the 
point  where  it  once  again  equals  rco  at  yr . 

2.  It  was  determined  that  y^  is  associated  with  the  presence  of  clay  particles.  From 
tests  on  four  different  sediment  beds,  each  with  the  same  dso  (=  0.83  mm)  and  a  different 
clay,  v|/min  was  found  to  vary  from  0.020  to  0.060. 

3.  It  was  determined  that  y^  is  less  significantly  affected  by  the  sand  size.  From  tests 
on  four  sediment  beds,  each  with  the  same  clay  (kaolinite)  and  a  different  d50,  y,™  was 
found  to  only  vary  from  0.045  to  0.060. 

4.  The  quasi-hydrodynamic  model  for  clay  lubrication  was  found  to  explain  the  observed 
behavior  of  zcr  withy.  This  model  predicted  a  value  for  y^  for  kaolinite  (0.04)  that 
matched  the  v|/min  values  obtained  from  the  data  for  kaolinite  beds  of  other  researchers; 
the  data  from  Torfs  (1995)  resulted  in  y^  =  0.02 ,  the  data  from  Huygens  and  Verhofen 
(1996)  resulted  in  y^  =  0.04 ,  and  the  data  from  Sharif  (2002)  resulted  in  y^  =  0.04 . 
This  result  was  explained  through  ymin  being  the  fine  weight  fraction  at  which  the 
thickness  of  the  clay  layer  is  close  to  the  size  of  the  asperities  of  the  sand  grains  (~  20 
urn)  at  which  the  transition  from  thin-film  to  boundary  layer  lubrication  occurs. 
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5.  The  shear  resistance  model  for  clay  lubrication  was  found  to  offer  an  explanation  for 
the  observed  behavior  of  rcr  with  vj/  in  the  boundary  layer  region  of  lubrication.  This 
model  predicted  a  value  for  y^  (0.01)  that  did  not  match  well  with  the  y^n  values 
obtained  from  the  data  for  kaolinite  beds  of  other  researchers  (0.02-0.04).  This  model 
failed  because  it  predicted  a  value  for  v|/min  that  resulted  in  a  clay  layer  approximately  3 
u,m  thick  which  is  much  smaller  than  the  size  of  the  asperities  of  the  sand  grains  (~  20 
u,m),  and  consequently  in  the  boundary  region  of  lubrication. 

6.  It  was  determined  that  \|/r  is  the  space-filling  weight  fraction  and  the  point  at  which  the 
behavior  of  xcr  of  the  sand  grain  matrix  starts  to  be  dominated  by  the  clay  suspension. 
From  tests  on  two  sediment  beds,  each  with  the  same  clay  (kaolinite)  and  the  same  d50  (= 
0.83  mm),  \|/rwas  found  to  be  the  point  at  which  the  pore  volume  not  occupied  by  clay 
becomes  negligible  (<  1%  of  the  total  volume). 

7.  Salinity  was  found  to  effect  sand  lubrication  in  the  boundary  region  through  changes  in 
flocculation  and  in  the  thick-film  region  through  changes  in  viscosity.  Viscosity  was 
found  to  increase  from  0.040  Ns/m2  to  1.60  Ns/m2  as  a  result  of  addition  of  salt  to  pore 
water,  and  from  0.001  Ns/m2  to  0.040  Ns/m2  as  a  result  of  increase  in  the  fine  fraction. 

8.  It  was  determined  that  the  shear  resistance  in  the  boundary  region  of  lubrication  was 
dramatically  decreased  with  the  presence  of  interstitial  clay  particles  due  to  an  increased 
number  of  load-carrying  contacts.  The  normal  load  at  a  discrete  contact,  N,  required  to 
cause  plastic  yield  was  found  to  increase  significantly  from  the  case  where  the  overall 
confining  stress  was  carried  by  a  sand  grain  matrix  in  which  the  grains  were  modeled  as 
spheres  with  1000  A  asperities  (N  =  33.2  g)  to  the  case  where  the  stress  is  carried  by  a 
kaolinite  suspension  with  a  porosity  of  0.48  (N=  1.25  x  1020  g). 
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7.3  Recommendations  for  Future  Work 

The  following  recommendations  are  made  for  further  work: 

1.  Consolidation  of  sediment  beds  composed  of  differing  sand  sizes  and  clay  types  needs 
to  be  conducted  to  determine  the  effect  of  consolidation  on  vj/^  .  Consolidation  has  been 
reported  in  several  previous  studies  (e.g.,  Mehta  and  Partheniades,  1982)  to  dramatically 
effect  the  behavior  of  clay  beds,  and  consequently  has  the  potential  to  effect  the  behavior 
of  Vmin  through  changes  in  the  clay  fraction  of  sediment  mixtures.  All  of  the  beds  in  this 
study  were  consolidated  for  a  period  of  72  hours  before  testing.  In  future  work,  beds  of 
similar  composition  (sand  grain  size,  clay  type,  and  fine  fraction)  should  be  consolidated 
for  varying  time  periods  to  determine  if  consolidation  has  an  effect  on  the  behavior  of 

T  min  ' 

2.  The  two  models  for  clay  lubrication  presented  in  this  study  give  differing  values  for 
\|/min,  and  this  discrepancy  warrants  further  study.  The  quasi-hydrodynamic  model  for 
clay  lubrication  predicted  a  value  for  y^  for  kaolinite  of  0.04,  and  the  shear  resistance 
model  for  clay  lubrication  predicted  a  value  for  \|/min  for  kaolinite  of  0.01.  A  potential 
cause  for  this  discrepancy  is  that  while  the  dimensional  analysis  that  leads  to  the  Stribeck 
curve  is  reasonable  for  the  "thick-film"  lubrication  regime  where  the  grains  are  not  in 
contact,  it  is  quite  possible  that,  in  the  lower  region  of  the  curve  (the  boundary  layer 
lubrication  regime)  where  the  grain  asperities  become  significant,  other  variables  such  as 
grain  asperity  size  might  become  contributing  factors.  For  example,  if  a  modulus  of 
elasticity,  E,  is  used  to  represent  the  grain  material  strength  (with  units  of  stress)  with  the 
length  scale  of  the  grain  asperities,  d,  then  a  potential  additional  dimensionless  parameter 
incorporating  the  normal  load  per  unit  width,  P,  would  be 
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(P/d) 

which  is  a  ratio  of  material  strength  to  applied  stress.  In  this  ratio,  an  increase  in  the 
normal  load  per  unit  width  would  be  equivalent  to  a  decrease  in  the  length  of  the 
asperities  (both  would  increase  the  stress)  or  a  decrease  in  the  grain  material  strength  for 
the  same  load.  It  would  consequently  be  expected  that  lower  values  of  this  parameter 
would  result  in  a  more  rapid  transition  to  the  "thick-film"  regime  of  lubrication  where  the 
grain  asperities  play  no  role.  If  this  were  the  case,  a  family  of  curves  would  be  expected 
that  should  all  converge  to  the  thick-film  curve  for  larger  values  of  dimensionless 
viscosity  (Figure  7.1),  analogous  to  smooth"  turbulent  flow  where  the  resistance  is 
governed  only  by  viscosity  and  the  roughness  plays  no  role. 
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Figure  7.  1  Stribeck  Curve  with  new  dimensionless  parameter. 
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Consequently,  the  application  of  various  normal  loads  per  unit  width  over  beds 
with  differing  grain  asperities  needs  to  be  conducted  to  determine  the  effect  of  normal 
load  per  unit  width  on  \\fmin . 
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